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PREFACE 

A number of original works on the application of 
mathematics to biology have been published during 
the past few years, but there is still no work of a 
synthetic and didactic type, bringing together all the 
research on this subject. M. Kostitzin's book fills 
this gap. 

M. Kostitzin is a very distinguished mathematician. 
He is thoroughly familiar with all the resources of 
analysis and its most recent applications. He has 
shown his ability and skill by solving problems that 
had defeated the efforts of his predecessors. 

But to deal with questions of an entirely novel kind, 
such as those with which M. Kostitzin is concerned, it 
is not enough to be a good geometrician ; other gifts 
and other abilities — found in the highest degree in this 
book — are necessary, as well as extensive knowledge 
beyond the regions of pure mathematics. 

M. Kostitzin himself clearly expresses this idea in 
the first chapter of this work, in which there will be 
found many general considerations on mathematical 
reasoning, on fundamental hypotheses in general 
science, on statistical methods and on the many sources 
of error in passing from experiment to analysis and in 
trying to arrive at general laws. This philosophical 
viewpoint at once gives the reader a favourable im- 
pression of the spirit in which the work is conceived. 

The book begins with general mathematical ideas 
on the frequency curves, the correlation of phenomena 
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and the differential equations that are encountered in 
biology. It next goes on to the general laws governing 
the circulation of organic matter, of carbon, oxygen, 
and nitrogen, which give a synthetic idea to life and 
its environment. 

One chapter is devoted to the Logistic Law, first 
stated by Verhulst and so often applied and verified 
by Pearl and other statisticians. 

The consequences of these problems of multiplica- 
tion, taking the most varied factors into account, form 
a new chapter. In particular a study of poisoning due 
to the products of metabolism will be found; the 
experimental work of Regnier and Lambin is familiar ; 
I had taken up the question from a mathematical 
point of view, obtaining an integro-differential equation 
which I had simplified and reduced to a differential 
equation. 

But it is to M. Kostitzin that the integration of this 
equation is due, an integration accomplished in a very 
ingenious manner, as will be seen in this book. This 
result has both a mathematical and a biological interest, 
and offers a fine example of the analytical audacity 
of M. Kostitzin. 

I have one comment to make on this point. In this 
particular problem it has been possible, by an artifice, 
to avoid integration of the integro-differential equation 
and to reduce it to a differential equation, but it is 
clear that integro-differential equations cannot always 
be avoided, and that they must play a part in biological 
problems. 

I shall not go into details of a host of examples in 
which the author studies population, the fluctuations and 
the growth of a population composed of several species, 
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taking into account food obtained directly from the 
environment or by the species eating each other. 

From this group of ideas we pass easily to the study 
of symbiosis and parasitism, and so to the growth 
of organisms and of tissues considered on the lines of 
a population. But the author carefully avoids push- 
ing this comparison too far; he shows the essential 
difference which must exist between a case in which 
the components of a group are independent and one in 
which they are not. In an organism the functions of the 
different elements, and their ways of nourishing and re- 
producing themselves, must be taken into consideration. 

In order to omit nothing that can be examined from 
a geometrical point of view, the author also deals with 
the form of living objects and the work ends with a 
general discussion of evolution and selection, in which 
mathematics plays a very important part. 

In the treatment of all these questions the author 
shows perfect familiarity with the many researches 
carried out by naturalists, especially recently, and his 
work successfully combines mathematical and experi- 
mental studies. To give an idea of the superposition 
and of the necessary fusion of fact and theory, the 
author employs a very striking image, drawn from 
romantic literature. 

It is unnecessary to add that the book will be of 
the greatest value to naturalists, and that it will give 
mathematicians a host of problems to investigate. 
From the point of view of the progress of mathematical 
biology, a new science that is developing so rapidly 
and so successfully, the publication of M. ICostitzin's 
book marks a very important date, 

Vrro Volterra. 



CHAPTER I 

INTRODUCTION 

i. The Purpose or this Book 

This book is intended for biologists who are to some 
extent familiar with mathematical analysis, and for 
those mathematicians who, being naturalists, do not 
believe that the sole aim of mathematics is the enun- 
ciation of formulae. I have tried to discover the rela- 
tions between my methods and their results, but in 
this region of learning, where so much is still to be 
done, the disparity between its methods and reality 
is often evident. From this there follows the appearance 
of duality which this book presents, as if it were com- 
posed of two different works whose pages had been 
intermingled. Literature includes a well-known exam- 
ple of this in the Kater Murr of the famous German 
writer, E. T. A. Hoffmann. In that book the matter 
of fact philosophy of a cat is contrasted with the shafts 
of wit, the dreams and the romance of Kreisler, the 
musician. Anyone who has read the Kater Murr 
realises that in spite of its rambling appearance, the 
book possesses a deep inner unity. There is romanticism 
and realism in my book, but there is unity as well. 
If I have not always succeeded in making it manifest, 
it is sometimes because the problem does not lend 
itself to this treatment and sometimes because the 
biological material was lacking. 

This little book differs from other works published 
under similar titles, which are rather text books of 
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mathematics for biologists than a mathematical study 
of biological problems. I do not deny the value of 
such works for information, for quick reference, etc., 
but I do not believe that such a selection of certain 
sections in one science for the use of workers in another 
can be successful. Each science has its own language 
and its own logic, and it is only by keeping to its own 
principles that it can retain its full value in applications. 
It is impossible to say what portions of a science are 
necessary and what are superfluous in any given 
application. It would be interesting to review all the 
mathematical text books published during the last 
fifty years intended for physicists and naturalists, 
beginning with that of H. A. Lorentz (1882). One 
would see how certain sections, long neglected by 
physicists and yet very important from the mathe- 
matical point of view, suddenly appeared and took a 
prominent place in these works. A biologist who has 
need of chemistry ought thoroughly to study that 
science, and not a few sections specially selected and 
adapted. In the same way, to be able to apply mathe- 
matical methods it is necessary to study what consti- 
tutes the very essence of these methods — their ideas — 
and not a few applications of the calculus. It is of 
greater value to consult a logically ordered text book 
than a compilation including scraps of analysis and 
odds and ends of the calculus of probabilities, decked 
out with biological examples that have often been 
dragged in by the scruff of the neck. 

In the following pages I have given a collection of 
biological problems capable of mathematical treatment. 
These are the problems of the multiplication of groups 
of organisms and of the organism itself considered as 
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a complex of cells or of organs, or as a thermal machine, 
or as a major unit composed of, but not reducible to, 
a number of minor units. From the mathematical 
point of view, all the problems considered present a 
certain unity. The fundamental ideas, the variables 
that are used, the co-efficients of the equations, are 
the result of long statistical elaboration, but the logical 
apparatus is purely analytical. The equations nearly 
all belong to the type of differential equations or 
integro-differential equations of the first order. That 
rich mine, opened independently by A. J. Lotka and 
V. Volterra is far from being exhausted. Numerical 
results are not always reached, but qualitative results 
are obtained, the only ones capable of verification, 
because the numerical biological material with which 
they deal is too uncertain or too heterogeneous for 
accurate corroboration. 

I should like to thank Madame J. Kostitzin for her 
invaluable help. I thank sincerely M. E. Rabaud and 
M. P. Montel, who agreed to the inclusion of this work 
in the Collection Armani Colin. I am most grateful 
to M. M. Prenant, who has done me the kindness of 
reading the manuscript, and to M. G. Teissier, whose 
conversation has given me many interesting sugges- 
tions, and to M. R6gnier, who sent me his experimental 
results. I am particularly grateful to M. Vito Volterra 
for the encouragement that he has never ceased to 
afford me. 



2. Mathematics and Natural Science 

The history and the logic of Science teach us that the 
application of mathematics to science is an inevitable 
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stage and a necessary condition of its development, 
Mathematics itself conceals under its array of logic an 
empirical foundation, sometimes so hidden by super- 
structures that mathematicians forget its very exis- 
tence. It is this basis which comes into play every 
time one appeals to common sense, it is this that is 
discussed in the study of the logical and psychological 
foundations of infinitesimal analysis. 

Everyone is agreed that figures have the right of 
entry into Natural Science, but when it becomes a 
matter of reasoning about these figures, or of dealing 
mathematically with them, there is encountered a resis- 
tance, a repugnance. Why is this ? It must be admitted 
that this resistance does not always come from narrow- 
mindedness or from conservatism. Reasoning in 
general does not frighten a naturalist, but mathe- 
matical reasoning startles him, because he is in 
the habit of verifying each step by experiment. In 
reasoning, experimental results are subjected to a 
series of logical operations. The accuracy of the 
ultimate result depends on that of the initial data, and 
also on the number and nature of the logical operations 
caixied out between the premise and the conclusion. 
In ordinary reasoning, this number is not very great, 
and successive stages are always verifiable. In mathe- 
matical reasoning, the steps are taken too quickly, and 
. a result is reached which may appear to be, and very 
often is, arbitrary or untrue. The reason for this is 
that in establishing a biological equation, the problem 
is simplified by sacrificing a number of factors or a 
number of details, and these sacrifices distort the 
results that emerge from the logical method. Inexact 
yet acceptable initial hypotheses are transformed into 



INTRODUCTION 17 

gross errors by a few turns of the crank of the logical 
apparatus. Simplification gives rise to paradox. This 
is a fact that cannot be denied, and one that is common 
to all experimental sciences. In purely mathematical 
problems, the mathematician is certain of his premise 
and the number and nature of the operations is of 
interest only from the aesthetic point of view ; but in 
dealing with applications of mathematics, the position 
is quite different. A long period of development has 
been necessary to secure unquestioning acceptance of 
the agreement between reasoning and experiment in 
mechanics, physics and astronomy. Perhaps in biology 
this development will be more rapid by virtue of an 
acquired impetus, In any case, long and continuous 
collaboration between mathematicians and biologists 
is necessary before reaching the same assurance as in 
the physical sciences, 

3. The Statistical and the Analytical Methods 

Mathematics gained entry into Natural Science by 
the statistical door, but this stage passes to the stage 
of analysis as in all quantitative sciences. The function 
of the statistical method is to clear the ground, to 
establish a certain number of empirical laws, and so 
to facilitate the transition from statistics to analysis. 
The labour of doing this is considerable, and important, 
but when it is done the first duty passes to mathe- 
matical analysis, which, at this stage in the growth of 
a quantitative science, is alone able to trace out the 
causes of phenomena and to deduce from them all 
the logical consequences. At this stage even approxi- 
mate hypotheses as to the true nature of the phenomena 
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are often more useful than empirical laws, calculated 
lege arlis, with all essential accessories. Thus Ptolemy, 
in his Optics, studies the problem of the refraction of 
light. He gives several results of experiments on the 
passage of light from air to water and from air to glass. 
His figures are quite sufficient to confirm Descartes' 
law. With these figures a modern investigator could 
construct a polynomial, reproducing them with exact 
precision, and could calculate their probable errors, 
but such a formula taken as a basis of research would 
be not only useless but misleading, and not even 
Descartes, had he been only a statistician, could ever 
have discovered the law of refraction. This example is 
not intended in any way to belittle the statistical 
method, but merely to deny it the character of a 
universal master-key, which is often attributed to it. 



4. Example : The Growth of a Population 

To illustrate more clearly the difference between the 
two methods, let us consider a biological problem, that 
of the growth of a homogeneous population. Let p 
be the number of inhabitants of a given area at a time t, 
and let the composition of the population be compared 
at times t and t -\-h. The difference between p(t -f- h) 
and p(t) depends on: 

1". The individuals born in that period (t, t + h). 

2. The individuals dying in the same period. 

3. The individuals entering the area (immigrants). 

4. The individuals leaving the area (emigrants). 

If the interval h between the two times is not great, 
it may he assumed (and statistics confirm this), that 
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all these increases and decreases are proportional to h. 
On the other hand, the same common sense, supported 
by statistics, suggests that the number of births and 
deaths are proportional to p. Thus we may write : 

No. of births in {t, t + h) ~ Nph ' 
No. of deaths in (t, t -f h) ~ Mfth 

where N and M are uncorrected co-efficients of birth- 
rate and mortality. As to emigration and immigration, 
these are also proportional to h, 'but their dependence 
on p is more complicated. If Ik and Ek be these incre- 
ments, we get : 

ft(t + h) - ft(t) = Nfth — Mph + lit — EA 

or, dividing by h 

fi(t + h)~ ft(t) 



h 



= Np - Up + I - E (r) 



i.e., the average rate of increase of the population. 
Proceeding to the limit as h approaches zero, we get 

ft'lt) = N£(*) - Mp{t) + I - E (2) 

It may be pointed out that this step is not altogether 
justifiable, since p is essentially discontinuous. This is 
a condition of great significance, the neglect of which 
may lead to fallacious results. It does not do so 
when p is sufficiently large, but when p is small, or 
when the density of the population is small, special 
precautions must be taken in order that the calculated 
results may agree with reality. 

The co-efficients N, M, I, and E are not constant but 
depend on both ft and t. In fact, the uncorrected 
fertility N decreases and the uncorrected mortality M 
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increases with growth of the population. Of all the 
hypotheses on the nature of this dependence the linear 
form 

N = n — vp, M ~ m + fip (3) 

is the simplest and is the most in accordance with the 
statistical data. Putting the above expressions (3) into 
equation (2) we get 

p'(t)^ E p(t)-hp*(t) + I~E (4) 

where 

e = n—m and h ~ v + fi {5) 

We shall call s the co-efficient of increase, h the 
limiting co-efficient, n the intrinsic co-efficient of fer- 
tility, v the limiting co-efficient of fertility, m the 
intrinsic co-efficient of mortality, p the complementary 
co-efficient of mortality, and all the co-efficients in 
biological equations of the type of equation (4) the 
vital co-efficients./ We shall see later that very 
simple hypotheses on emigration and immigration 
enable us to replace E and I by appropriate functions 
oip. 

The vital co-efficients depend on all the factors that 
can influence a population, climate, area colonized, 
food, social organization, etc. Further, it often happens 
that the simple relations of equation (3) are not always 
verified, yet at the same time equation (4) may be 
accurate while equation (3) is not. In other cases more 
complex hypotheses are necessary, and equation (2) 
assumes a more complicated form. It will be shown 
later that the simple form of equation (4), with all 
the consequences that follow from it, conveniently 
expresses a very considerable number of biological 
facts. 
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5. Historical Remarks 

We conclude this chapter with a few historical 
remarks. Natural science has used the statistical 
method for a long time, but it is to the great English 
scholar Karl Pearson that we owe a complete revision 
of the method and the foundation of a new science- 
biometry, 

Equation (4) was first given by Verhulst in 1845 in 
discussing an isolated population. Since the works of 
Verhulst special biological problems have from time 
to time been investigated by the analytical method, 
. but it is to Lotka, and independently to Volterra, that 
the honour belongs of having systematically applied 
the method to more general problems. The work of 
Lotka in physical biology has inaugurated a bio- 
mathematical movement in Anglo-Saxon countries. 
The brilliant work of Volterra, so fertile in mathematical 
suggestion, has started a similar movement on the 
continent of Europe. Pearl's half theoretical, half 
experimental researches on the laws of population 
growth have provided the basis of modern demography. 
J. B. S. Haldane's very important contribution to the 
mathematical theory of Natural Selection should also 
receive mention. 



CHAPTER II 
SOME MATHEMATICAL IDEAS 

In this chapter we shall briefly consider some mathe- 
matical ideas that will be useful later. 

i. Frequency Curves 

Suppose that we have to study statistical data 
referring to any biological character, such as size in 
a group of animals. We subdivide the 'range be- 
tween the smallest and 
the largest into a number 
of equal intervals and 
count the number of 
- animals whose size falls 
within each of these 
divisions. These figures 
can be represented by a 
polygon of frequency in 
'" which the abscissae are 
the sizes and the ordi- 
nates are the numbers of 
animals included in the 
corresponding size- 
interval. This gives a 
polygon which for the 
best-known infusorians, 
■p, Paramecium, has the 

M I.-FREQ^CY POLYGON :J araCteristiC i0 ™ ° f 

Size o! Paramecium (after Jennings). Fig. I. 




SOME MATHEMATICAL IDEAS 23 

This polygon shows that the smallest and the largest 
individuals are rare, that the average size is the com- 
monest, and that the polygon is to all intents and 
purposes symmetrical about a median line. This 
distribution can be represented by many curves, but 
for many reasons the Gaussian curve is chosen. Let 
I be the size, l m the average size, N the total number of 
individuals and Aw the number whose size falls within 
the interval (/, I -f- A/). Then 

A W = ^*-W- '•>'** (i) 

h 

A frequency curve of this kind is completely deter- 
mined by three parameters, of which the first, l m , gives 
the position of the peak, the second, N, the total num- 
ber of individuals, and the third, h, the distribution 
about the average. Law (i) is that of the distribution 
of chance errors about a mean. By supposing that the 
distribution of size verifies this law we implicitly admit 
that this biological character only obeys the law of 
chance, that is to say that the distribution of size 
about the average l m is governed by no special factor 
and that the statistical material is homogeneous. 

In. the special case of Paramecium these conclusions 
are not true. The material is heterogeneous, for 
Jennings' statistics do not distinguish between indivi- 
duals of different ages, and £.ge is a factor which cer- 
tainly influences the parameters of curve (i). In fact, 
curve (i) more or less resembles the polygon of Fig. i, 
but does not represent it satisfactorily, and one gets 
the impression that a curve of the type 

An = A x e~ **'<'- <■>' Al + A 2 e- *'""" ( »>' Al + A a e~ *•'('- «• M 
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would be preferable. Consequently it is always neces- 
sary to standardize the material, that is to say to study 
the distribution of characters among individuals of the 

same age, the same race, 
etc. In many cases the 
heterogeneous nature of 
the material is expressed 
in the frequency polygon 
itself, in the form of 
several peaks. This is true 
of a mixed population 
composed of two forms 
of Paramecium (Fig. z). 
Peak a corresponds to 
Paramecium aurelia, and 
peak A to Paratnecium 
the distribution can be 
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Fig. 2. — Frequency Polygon : 

Body -length of Paramecium 

(after Jennings). 

caudatum. In this case 
represented by the curve 

A« = A x e" v "~ '•'' M + A/- *•*<'- «' M (2) 

Frequency curves are often found with two or three 
peaks, but it is not always easy to ascertain the reason 
for the complexity. Sometimes it is a difference in 
age, as in birds with a seasonal diet ; the age-groups 
are separated by a considerable interval and the average 
size of the young is not the same as that of the adults. 
Sometimes it is sex dimorphism. In some cases the 
activity of predatory foes selectively destroys the 
individuals of particular size- or age-groups and so 
creates an apparent heterogeneity in a homogeneous 
population. Again, features of the environment, such 
as temperature or humidity, act upon biological char- 
acters and if the statistical material reflects the whole 
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of the results on a homogeneous population the fre- 
quency curve may have several peaks. 

From a mathematical point of view, the operation 
of dissecting a frequency curve with several maxima, 
that is to say the determination of all the parameters 
in an expression of type (2) presents considerable 
difficulty and uncertainty even in the simplest cases. 
This uncertainty is further increased by the fact that 
the frequency curve of Gauss is not the only one: it 
corresponds to the simplest schemata of probability 
only. Such a schema of probability is a vessel filled 
with balls of two or of different colours. One ball is 
taken from the vessel and is then returned to it so 
that the conditions in subsequent tests shall be un- 
changed. In such a series of experiments the distribu- 
tion of deviations from the mean corresponds to law (1) . 
Hence Gauss's law is valid only when the tests are 
completely independent of one another. 

Other methods of conducting the experiment can 
be suggested, the balls not being returned to the vessel, 
or being returned according to certain more or less 
complex rules which would produce a dependence 
between successive tests, or replacement may be de- 
layed, conditioned, etc. 

In this way frequency curves differing from 
(1) are obtained. Pearson, in studying the distri- 
bution of biological characters, realized the impossi- 
bility of representing them all by curve (1) or its 
variations. He inti*oduced a new group of frequency 
curves by generalizing the differential equation of 
the curve 

4>=*Ae- kH *-'« ) ' (1') 
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Taking the derived logarithm of (i') we get the 
equation 

£ = *»{*»- x)- 3 ^ (3) 

This equation shows that the relative growth of <■/> is 
proportional to , the difference (x m —x). Pearson as- 
sumed that the co-efficient of proportionality <x is not 
constant, and substituted for equation (3) 

f = x m ~ x 

</> oc -j- a x x ~\- a, z x 2 

There are several arguments in favour of this general- 
ization; it is sufficient to say that urn schemata can 
be constructed to reproduce the majority of Pearson's 
frequency curves. The curves resulting from the 
solution of equation (4) have very different appearances 
and can account for several distributions where curve 
(1) fails to do so. Further, the many factors deter- 
mining size and other biological characters are related 
to, and dependent on, each other. Hence all the curves 
that analysis produces must be admitted on equal 
footing. For each curve schemata of probability can 
be found, justifying its use -for the representation of 
distribution of characters. For a mathematical natur- 
alist this implies freedom of choice and the end of the 
monopoly of type (1). 

2. Evolution of a Frequency Curve 

A frequency polygon may vary with time — for 
example, it may vary with age if one is studying the 
distribution of a character in a homogeneous popula- 
tion. This polygon may be represented by a frequency 
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curve of a type determined by co-efficients varying 
with age. Suppose, for instance, that the character 
studied is human stature, and that we have frequency 
polygons for each age of human life. We can represent 
these polygons by Gaussian curves and calculate for 
each age the average stature, the number of individuals 
of this age and the dispersion. All these parameters 
have a real biological meaning and the study of their 
variation with age gives valuable results. 

However, this procedure will not serve in every case. 
It often happens that the frequency curve of a certain 
type can represent only a part of the development of 
a population, and it is necessary to look for another 
expression. In other cases the existence of a simple 
causality determines the distribution and makes statis- 
tical representations valueless. The following is a 
biological example in which this condition obtains. 

The females of a species of moth Lithocolletis populi- 
foliella Tr. deposit their eggs on the leaves of poplars. 
Polejaeff has studied the relation between the number 
of eggs laid on a leaf and the area of the leaf, and also 
the relation between the number and the mortality 
of the caterpillars on a leaf. Fig. 3 gives the three 
successive stages of contamination of the leaves by 
these moths. Let / be the area of a leaf and fi the 
number of eggs laid on a leaf. Fig. 3 shows that on 
three successive ovipositions — June 13, June 17, June 
2i — the number of eggs per leaf is proportional to the 
area 

i»-A/, (5) 

but A, the co-efficient of utilization, varies with the 
date — > 
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Date 


Junk 13 


June 17 


June 21 


A 


•42 


■9.1 


1-28 



There is a certain slowing down in the rate of utilization 
of the leaves. 
How are these facts to be interpreted? The first 
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Fig. 3. — Infection of Poplar Leaves by Moths 

explanation, suggested by equation (5) is the existence 
of a relation between the number of eggs laid and the 
path taken over the leaf by the insect during oviposi- 
tion. This obvious explanation is not in accordance 
with the fact that the moth does not lay more than 
two eggs a day, while on the first day of laying the 
number of eggs varies from 4 to 13. Hence the eggs 
on any one leaf have not all been laid by a single 
female. Moreover, Fig. 3 shows that the number of 
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small leaves used increases between the first and fifth 
days of oviposition. The ratio between the maximum 
and minimum area per egg decreases, as the following 
figures show — ■ 



Date 


June 13 


June 17 


June 21 


Max. area 
Min. area 
Max. area 
Min. area 


6'8o 

•40 

17 


2-30 
•20 

"•5 


I-IO 

•25 
4'4 



On the ninth day of egg-laying the utilization of the 
leaves is greater than on the first day, but is not 
complete. PolSjaeff sees in these figures a sign of 
preference for large leaves. In my opinion the moth 
tries above all to ensure a sufficient food supply for its 
caterpillars', on the first day of oviposition it visits 
large leaves, avoiding, as far as possible, leaves that 
are already occupied. On the following days eggs are 
laid wherever there is space available, with no prefer- 
ence for large or small leaves. The increasing difficulty 
in the search for an unoccupied area is expressed in the 
retardation of increase of the co-efficient of titilization A. 
From the insect's point of view this method of egg-lay- 
ing is evidence of very bad adaptation. The presence 
of many caterpillars on one leaf facilitates epidemics, 
attracts parasites and predators, and creates many incon- 
veniences for the caterpillars themselves. Chance plays 
a great part in the wanderings of these insects, but the 
final result is subordinated to a very simple causality. 



3. Correlation and Causality 

Search for a causal bond between two phenomena 
is not easy. Galton has suggested a practical method 
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which is of service in certain cases. Suppose that 
observation gives two series of numbers 

yv y* ■ ■ -,y n \ 

between which a relation is suspected of the nature 

F(*. y) = o (7) 

Suppose that the function F(x, y) be expanded accord- 
ing to Taylor's theorem — 

F(x, y) = a + a-iX + a 2 y + . . . (8) 

and that terms above the first order are negligible. 
Equation (7) can then be re-written 

«o + «i# -I- «ay = (9) 

If this expression of the relation between x and y is 
permissible, the points 

(%i,yi), K.y 2 ). • • ■ . (**.y«) 

cannot lie far from the straight line (9) . 

Imagine such a straight line and estimate the sum 
of the distances between it and the points [x k , y k ) 

f - V (a ° + UiX * + °*y*)* / Tn » 



For simplicity, put Xs=—yx k Y = - ^> y k 

k-x k- 1 

A -:2-' B -=2*' c -;5>- 

k *• x A - 1 A _ j 
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cos Q — 



VV -f a£ Va* + 



a 



sin .0 — — 



Equation (10) then becomes 

E = \ {d — x h cos — y k sin 0) 3 

or E = n(d* — zKd cos — 2 Yd sin + A cos 2 

-I- 2C cos sin 6 + B sin 2 0) (10') 

Seeking the values of d and of 6 which make this 
expression a minimum, we get by equating to zero the 
derivatives of E in respect of d and d 

d — X cos — Y sin = o (11) 

A — B 

dX sin 8— dY cos 0+C cos 20 sin 20 — o (12) 

Equation (11) shows that the point (X, Y), the centre 
of gravity, of the points [x k , y k ) lies on the line (8) . 
Equation (12) gives the value of in terms of these 
points: a(C -XY) 

tan 2 ° = A-X»-B + Y» " (I3) 

Further, the minimum value of E satisfies the 
equation 
E 2 -E'»(A-|-B-X a - Y a ) 

+ <n 2 [(A - X») (B - Y») - (C - XY) a ] == o (14) 
It can easily be shown that 

A > X a ; B > Y 2 ; (A- X 2 ) (B - Y 2 ) > (C - XY) a 
and therefore we may write 

r*(A - X 2 ) (B - Y 2 ) - (C - XY)« (15) 
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in which the co-efficient r" is equal to or less than 
unity. Equation (14) becomes 

E a ~EnfA + B-X 2 - Y a ) 

+ »»(i - r 1 ) (A - X 2 ) (B - Y a ) = (16) 

This equation has two positive roots, of which the 
smaller gives the required minimum. Suppose that all 
the points (x kl y k ) lie on the line (9). In this ideal case 
E = 0, but equation (16) shows that E can only dis- 
appear when a — t 

Further, in the case of symmetrical distribution of 
the points (x k , y k ) about the point XY we get 

In all other cases r 2 is less than unity, and it is 
natural to use this co-efficient to measure the possi- 
bility of representing a group of points by a straight line. 
The co-efficient r is called the co-efficient of correlation. 

The co-efficient of correlation is zero when it is 
impossible to represent the observations (6) by a 
straight line, but this does not mean that the magni- 
tudes of * and y are independent : it simply means 
that the function F(x, y) cannot be even approximately 
replaced by the linear expression (9) . This important fact 
is often forgotten and it is always worth remembering. 

4. Two Connected Series of which the 
Correlation is Zero 

We give an example of two connected series of which 
the correlation is zero although the variables are not 
independent. Suppose that two biological variables 
are functions of the time t. Let a> be an interval of 
time during which the variables x and y have been 
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observed. The sums X, Y, A, B, and C may be replaced 
by the integrals 



A = i 

co 





* <o "I" <0 




* f, + 01 


I 




I 




x = - 


x.dt, 


Y = - 


y.*. 


ft) 




co 




1 


'o 


* 


'« 


* 'n + <" _ 


* *o -1- « 




I 




I 


xKdt, B==- 


y*.dt, 


c = - 




<o 


CJ 



xy.dt. 



Suppose next that x and y are periodic functions of 
time for the interval co and that 

C = XY (17) 

In this case the co-efficient of correlation is zero, 
however closely x and y may be united by a causal 
bond. This is what occurs in the case of a mixed 
population composed of two species one of which feeds 
on the other and where the densities of the populations 
of the two species are sufficiently small for the limiting 
terms in each species to be negligible. Let x be the 
number of the species preyed upon and y the number 
of predators. The variables x and y satisfy a system 
of differential equations 

%' = e x x — h 12 xy 

y' = - e a y -I- h 2l xy. 
The term sjX implies that in the absence of predators 
the first species will increase according to the Malthusian 
law ( 

the term — e 2 v implies that in the absence of food the 
second species will be driven to extinction according 
to the exponential law 

the term — h n %y implies that each encounter between 
a predator and its prey is followed by destruction of 



(18) 
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the latter ; the term h 21 xy expresses the benefit which 
the predators gain from this. It will later be seen that 
a simple relation exists between x and y, 

X HyHe. - *«*- h# — h, (19) 

which shows that % and y are periodic functions of t. 
Let a> be this period. Integration of equation (18) gives 

X = X Y = il, XY-C (20) 

"21 "12 

or, in other words, the relation of equation (17) is 
verified and the co-efficient of correlation is zero, and 
nevertheless x andy may be connected by equation (19). 
The case represented by equation (18) is not a special 
one, nor is it purely theoretical. Periodic or periodo- 
morphic* oscillations are in fact observed in mixed 
populations and in this case it would be actually 
misleading to apply the method of correlation, with 
the formidable calculus that it requires, to the search 
for a dependence between biological variables that 
cannot be detected in this way. This does not mean 
that the method of correlation is wholly useless, but 
it greatly restricts its application and provides a 
further argument for the analytical method. 

5. Systems of First Order Differential 

Equations 

The following equation has been obtained for a 

homogeneous population : 

p' = ep - hp* + I - E. (21) 

Suppose that the terms I and E, which express immi- 

* By this term I describe oscillations that are not strictly periodic. 
I wish to avoid the terms pseudo-psriodio and quasi-periodic, as 
these have been used for some time in a specialized sense. 
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gration and emigration, depend only on the function 
p(t). Equation (21I may then be written 

£'(<)= W)] _ (22) 

Similaily, in the case of a mixed population composed 
of different groups p v p%, . . . , p n , these variables 
satisfy differential equations of the first order. 

Pk = HPv P2. • ■ • . Pn) (A == 1, a, . . . , «.) (23) 
In most cases the function <f> k is a polynomial of the 
second order in relation to the functions p x , p 2 , . . . , p„. 
Even with this simple hypothesis, equation (23) can 
be successfully resolved in a few simple cases only. 
Yet in the absence of a numerical solution, we must 
rest content with a qualitative study which can none 
the less give very interesting biological results. 
Consider the values 

Pi = a v A = H> • ■ ■ > P« = «» (24) 

which solve the algebraic equation 

<h = o, f /> 2 = o, . . . , <j> n — (25) 

If pi, p z , . . . , p n are given initial values equal to (24), 
these functions remain constant and the mixed popu- 
lation does not change; system (24) then represents 
a stationary state or a state of equilibrium. Several 
equilibrium states may exist and it is essential to dis- 
tinguish those that are stable from those that are not. 
Stability may be variously defined. We may speak 
of stability when small variations in the co-efficients 
are followed by only small variations in the results. The 
study of this stability will be kept for a later chapter. 
We may speak of stability when in the neighbourhood 
of a stationary state (24) small variations in initial 
values lead to solutions tending, when t-^-oo, towards 
the values (24). This stability will be considered here. 
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Suppose p v p % , . . . , p„ to be the initial values in 
the state of equilibrium (21) and let 

Pi — «i + h> Ih = a 2 + 4. - • • » fin = «» + L 
where £ v £ 2 , . . . , £ „ are of very small magnitude, so 
that their squares, products, etc., are negligible. Sub- 
stitute these values in equations (23). Since the 
functions p v </> 2 , . . . , 4> n are polynomials in relation to 
Pi, fiv ■ ■ • ,p« we ma Y write, taking equation (25) into 
account „ 



and equations (23) become 



,pn} ~Z,3aJ n 



£' 
C fc 



Tt ( - <*-■■ "'» 



(26) 



These equations can be resolved into exponential 
functions. Suppose that 





+ K/' nt 


(ft - 1, 2, 


• • 


• > ») (27) 


The exponents {.i v fi z , . . 
equation of the nth order. 


. , /.<„ satisfy the algebraic 


3« x 


~t*> 


3rt 2 


3«n 


* 




3% 




3^2 


_ 3 £a 
3«„ 




■ « (28) 


3« n 




30 2 

3« 2 


' ' 3«„ 


/* 





This equation, known as the characteristic or secular 
equation, has w roots, ^ l( fi a , . . . , /i„. 
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These roots may be real or imaginary. If they are 
all negative or with negative real parts, the equilibrium 
is stable, for the fluctuations (27) tend towards zero. 
If only one of the exponents is positive or zero, or 
with a real part positive or zero, the equilibrium is 
unstable. 

Consider, as an example, the case in which n = 2. 
Equation (28) becomes 

or fi'—pC + D^o (29) 

C 3<$j 3c/> a t-\ tyl <V2 d<f>l 3^2 / \ 

3% 3a 2 3fl x 3« 2 3a 2 3« x 

When C < 0, D > 0, E = C 2 — 4D > o, there is a 
node, that is to say the two solutions of the character- 
istic equation are negative and the point (p v p^ 
approaches [a v fl 2 ), following an arc (Pig. 4A). 

When C < 0, D > 0, E = C 2 — 4D < 0, there is a 
focus, that is to say the roots of equation (29) are partly 
real and negative, and the point [p lt p^j tends towards 
(a v # 2 ), describing a spiroidal curve (Fig. 4B). 

When C = 0, D>o,E<o there is a cycle, that is 

to say the two solutions of equation (29) are purely 

imaginary and the point (p v p^ describes a closed 

curve round the point (a v # 2 ), (Fig. 4c). This is what 

■_ happens in the case considered in the preceding section. 

Let us return to the general case. Equations (25) 
have several sets of solutions of type (24) and each of 
these corresponds to a stationary state. Some of them 
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are composed of negative or of complex expressions 
and hence have no interest here; the others may be 

stable or unstable. 
In the case in which 
there are several com- 
patible stable states 
the final end of the 
biological system de- 
pends on the initial 
values p^o), fafp), 
■ . ■ , p„ip) in the 
sense that each final 
state is surrounded 
by a region under its 
influence and if the 
initial state is in one 
of these regions, the 
corresponding final 
state is the limit. 

The variables p k 
are generally posi- 
tive. What happens 
if one of them van- 
ishes? This implies 
that one of the groups 

_ _ composing the mix- 

FiG. 4.— Curves 03? Relative Growth p a nnnnlitinn i= »v 
fob Two Groups, p l and p 3 ea P°P"iatlon IS ex- 

tinct. For its re- 
appearance new circumstances are essential, for ex- 
ample, an immigration or the hatching of some eggs 
which might have remained hidden. If there is no 
reason for the reappearance of a lost group, equations 
(33) must be revised to determine the effect of the 




SOME MATHEMATICAL IDEAS 39 

extinction of the group. It is not sufficient to replace 
by zero the variable that corresponds to the lost 
group for it may have an indirect influence on the 
vital co-efficients of the other groups in the population. 

If the numbers of one group approximate to zero, 
equations (23) must be handled cautiously, for bio- 
logical variables are essentially discontinuous, and the 
reduction of a group to a few individuals may be equiva- 
lent to its disappearance, which equations (23) do not 
provide for. On the contrary, biological equations, in 
order to be logical, must be concerned with large 
numbers, and for small numbers the form of these 
equations would be different. 

Other comments will be made in the examination of 
concrete biological problems. 

6. Biological Differential Equations 

We conclude this chapter by some consideration of 
biological differential equations. In almost every case 
that will be studied these equations belong to the type 



n n 

P'n : =e<-ic+ / ends— Pk / h>P> 



(k = i, 1, . . . , n) (31) 

The term a fc expresses the constant sum of emigra- 
tion and immigration. The co-efficients s kk almost 
always represent the intrinsic multiplication of the 
group ft k . The term B klt p e , in which h ¥> s, corresponds 
to possible transfers from one group to another, as for 
example in the case of a species divided into age-groups. 
The term h kk p z k represents the limiting action on the 
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group p k of all limiting factors save the action of other 
groups in the same population. The term h k ,p k p a repre- 
sents the interaction between the groups. It will be 
seen that other interpretations of these groups are 
possible. 

These equations produce 2" separate points, which 
represent stationary states. We are here concerned 
with a fairly common special case, that of an isolated 
system composed of a certain number of groups, without 
the possibility of passage from one to the other. In 
this case we have the equations of Volterra : 

p\ ~p!c ( s k ~~ /> Kp, J (A = i, a, - • • , «) (32) 

\ s •== I ' 

The 2" equilibrium states always remain theoretically 
possible. It is easy to enumerate them all. First we 
have the state 

1\ = 0, P„ = 0, . . . , P n = (33) 

which is only stable if all the co-efficients of increase 
are negative. 

e fc <o (ft = I, 3, . . . ,«). (34) 

Then there are n stationary states of the type 

n kk 

. . • , P. » (35) 
These equations imply that all groups save one dis- 
appear. Some of these states may be negative and of 
no interest to us. A state of type (35) is stable if the 
inequalities 

e fc > o, h kk > o, 
eJhk<h, k s k , (s = i, 2, . . . , It— x, k + 1, . . . , 11) (36) 
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are satisfied. These n states may be compatible with 
each other. Consider, for example, two states of the 
order k and I. We may have simultaneously 

£k > 0. hh > 0, eA-fc < h fK £ h (s ¥> k) 
e, > 0, h n > 0, e„.h u < htft {s' ¥- 1) 

It follows that 

eJh-k < fh k e k 

and therefore that 

hdhi < hnJhi 

These inequalities are not contradictory. On the other 

hand, state (35) is incompatible with state (33). 

n(n — 1) . , 

1 hen there are — - stationary states of the type 

p _ fc l aa E %h\ p _ £ a^ii ~ £ i%2 / , 

''ll"22 — "12"21 "U"22 — "12"21 

P 3 = 0, . . ,,P, = o. 

We may continue in this way and obtain all the 
equilibrium states of the system. Without pausing 
over the stability and the compatibility of these states, 
we shall only mention the equilibrium state in which 
all the groups of the system share, This state is defined 
by the linear equations 

n 

H = y A*.P. (A = I. 2, . . . , n) (38) 

JL i 

This state is stable if the secular equation (28) has 
all its roots ( ; « x , jj, z , , . . , ft, n ) negative or with nega- 
tive real parts. 



CHAPTER III 

LIFE AND ENVIRONMENT 

i. Circulation of Organic Mattku 

Let us consider the life process as a whole, and in its 
relations with inorganic matter. Snch matter exists 
in several forms, but the atmosphere is the chief source 
from which living organisms absorb the elements 
necessary for the formation of complex compounds 
and of living matter. The process of the assimilation 
of the absorbed matter consists of several stages. 
Plants absorb atmospheric carbon, converting the 
carbon dioxide into starch and liberating the oxygen, 
and comparable to this are the organisms living in the 
sea which absorb carbon dioxide and convert it into 
chalk. Nitrogen is a more inert component of the 
atmosphere than oxygen or carbon dioxide ; it is fixed 
by soil bacteria living symbiotically with certain plants. 
Thus nitrogen becomes available for further syntheses, 
ending in the formation of living matter. Other micro- 
organisms may be mentioned, which fix other elements 
necessary for the production of organic matter, as 
well as natural factors, such as thunderstorms, 
which facilitate the absorption of nitrogen by living 
organisms. 

In spite of their variety! all these processes can be 
summarized in the form of simple differential equations, 
capable of supplying information as to the nature of the 
life process in its general relations with the inorganic 
world. 

43 
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We may rapidly review the principal sources and 
ultimate destinations of each substance taking part in 
these reactions. 

1. Atmospheric Oxygen. Owing to its chemical affin- 
ities, the oxygen of the air cannot be considered as one 
of the original constituents of the earth's atmosphere ; 
it is rather the result of the activity of living organisms. 
It is evolved by plants in the process of the assimilation 
of atmospheric carbon, and absorbed by all living 
objects in respiration as well as in human industries. 
It is impossible to point to any source of this gas other 
than the chemical actions occurring in plants. 

2. Atmospheric Carbon. This exists chiefly in the 
form of carbon dioxide. Its sources are volcanic action, 
respiration in animals and plants, decomposition of 
organic matter and inorganic oxidation of the earth's 
crust accompanied by the evolution of this gas. The 
catises of its removal are the absorption of carbon 
dioxide by plants and by chalk-producing organisms. 
An essential difference may be noted between the 
organisms that directly absorb carbon dioxide and 
those that only make use of carbon already assimilated. 

3. Atmospheric Nitrogen. This is assimilated by 
organisms with the help of microbes in the soil and, 
perhaps, liberated in certain decomposition processes. 
It cannot be asserted that nitrogen is not evolved in 
some chemical reactions in the earth's crust, but there 
is no direct proof of this. 

4. Atmospheric Hydrogen. The lower layers of the 
atmosphere contain only traces of this gas, which takes 
part only in infinitely small quantities in the life 
processes. 

5. Water, Despite the importance of water in organic 
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life, the study of its circulation does not yield such 
useful information as does that of the atmospheric 
gases. 

6. Light Elements of the Earth's Crust. These ele- 
ments enter into the composition of living matter to a 
certain extent, but it would be premature to consider 
their circulation. The role of water as solvent, as 
carrier and as distributor is of fundamental importance 
in this, but the chief methods of circulation are far 
from obvious. Hence we are here confined to the 
mathematical study of the circulation of the three 
principal atmospheric gases. 

2. The Carbon-Oxygen Cycle 

Consider first the particular problem of the carbon- 
oxygen cycle. Let 

x = the total weight of free atmospheric oxygen. 
y 1 — that of atmospheric carbon dioxide. 
y 2 — that of carbon dioxide dissolved in the sea. 
s — that of these substances in chalk-producing or- 
ganisms, 
w — that of these substances in animals. 
v = that of these substances in plants. 
s = that of these substances dispersed through the 
earth's crust. 

The relations between these variables can be sum- 
marized by a simple diagram (Fig. 5), or by the follow- 
ing differential equations : 

*' = — «»« + Ki~ ««)w — % (1) 

y\ - <x 32 « — (au— a 4a )w + e'gt (2') 
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y\ = e"62 - y*y» 
«'=—(— a 13 + a 32 + <%)« -h fiuv 

v' — (a 14 — a 41 + a 24 — « 43 — - a iS )v — /?«v 

s' = a 35 u + k^u + a 65 ^— e' 5 a— £ "s2 + e is 

z' — — a^ + r^a 

The structure of these equations is very simple. 
Each term corresponds to one of the arrows in Fig. 5. 
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(2") 
(3) 
(4) 
(5) 
(6) 




Fig. 5. Circulation of Oxygen and Carbon 

Thus equation (1) shows that oxygen is consumed hy 
the respiration of animals (— a iz u), and of plants 
(— oi. i4 v), by inorganic oxidations of the earth's crust 
(~~ %), and that it is liberated in the process of assimi- 
lation by plants. Equation (2') shows that the atmo- 
sphere receives the carbon dioxide evolved by animals 
(a 32 w) and plants (oc^m) in the processes of respiration 
and decomposition of living matter, and that plants 
assimilate carbon dioxide (— a Zi v) : the term e' 62 repre- 
sents what the atmosphere receives as liberated carbon 
dioxide, from volcanic eruptions and other chemical 
reactions in the earth's crust. Human industrial 
activity, which is rapidly destroying mineral sources 
of energy, is wholly neglected: in a few centuries 
these resources will be brought to nothing, and mankind 
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will have to seek other supplies of energy. The carbon 
dioxide set free during this interval will be distributed 
between the atmosphere and the sea water, and hence 
our equations can be applied without the introduction 
of supplementary terms; it is enough to replace the 
numerical co-efficients by others more closely corre- 
sponding to the partial pressure of atmospheric carbon 
dioxide. 

Equation (2") refers to the hydrosphere. The total 
weight of carbon dioxide dissolved in the ocean is 
thirty to forty times greater than that in the atmo- 
sphere. To simplify these equations it is assumed that 
this proportion does not vary and that the ocean 
receives, through the intervention of the atmosphere, 
a constant amount, s" B2 , of carbon dioxide supplied 
by the earth's crust; the term (— yzy 2 ) represents 
dissolved carbon dioxide absorbed by marine chalk- 
producing organisms. 

In equation (3) the terms (« 13 «, — a 3a w) represent 
the amounts that aniznals add to or extract from the 
atmosphere; the term puv shows that animals live at 
the expense of plants and that process is regulated 
by the so-called principle of encounters between pre- 
dators and victims ; the term (— a 35 u) represents the 
fertilization of the soil by the products of animal 
metabolism. 

Equations (4) to (6) do not differ in their mode of 
construction from the above. 

The most difficult problem, and one that is not 
exactly represented by the equations (2') and (2") is 
that of the relation between atmospheric and oceanic 
carbon dioxide. According to Schloesing, the masses 
of oceanic water control the carbon dioxide in the sense 
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that the water loses carbon dioxide when the pressure 
of that gas in the atmosphere decreases, and absorb 
it in the reverse circumstances. This action ensures a 
degree of constancy in the co-efficients oc 32 , c/. Si , ct 42 , and 
consequently in the co-efficients « 13 , a 14 , a 41 , but makes 
uncertain the values of e' 5a and e" E2 . As for oxygen, 
these co-efficients are not invariable, and depend on 
the pressure of this gas in the atmosphere, but the 
weight of atmospheric oxygen is sufficiently great to 
ensure relative constancy of oc 32 , a 24 , a. l2 , a 13 , a M , and a tt . 
Thus the assumption of these constant co-efficients does 
not distort the truth very much, and the conclusions 
drawn from equations (1) to (6) remain valid, 

It should be noticed that each term appears twice 
in these equations, once with a positive and once with 
a negative sign. Hence, on adding these equations we 
get: 

*' + y\ + y\ + «' + «' + s' + z' = (7) 

and therefore 

X _f_ y l +y 2 + u + v + S + Ze -M (8) 

which expresses the commonplace of the constancy of 
the weight of the globe. 

Equations (1) to (6) may be arranged in three 
groups : 

1. The variables x, y, and s do not occur on the right- 
hand side. Thus equations (1), (3'), and (5) enable us 
to calculate these variables when all the others are 
known. 

2. Equations (2") and (6) contain only two variables, 
y z and z. 

3. Similarly, equations (3) and (4) contain only the 
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two variables u and v. Starting with equations (3) 
and (4) , we may write : 

u' — «(— X -f flv) and v' = v(fi — j3u) (9) 

where X = a 32 -f- a 35 — a 13 

and fj. = oc 14 — a 41 -|- a 24 — a. 12 — a d6 

The parameter X is positive, for in the absence of 
plants animals cannot exist ; the pai-ameter /.t is posi- 
tive, for in the absence of animals plants develop almost 
without hindrance. 

Equations (9) are only special cases of equation (18) 
of Chapter II, and equation (19) of Chapter II then 
becomes 

«V = Hfi"" + ^ (10) 

where H is a constant depending on the original values 
of w and v , 

H = « ' , V«~ /to, '" / *'• 

Time does not enter explicitly into equation (10). 

To each positive value of H not exceeding a certain 

limit there corresponds in the graph of the variables 

(m, v) a closed curve (Fig. 4, (hi) ) which is not cut by 

the analogous integral curves. This implies that u and 

v are periodic functions of the time t. Let a> represent 

the period, and let 

u X 

u = r ~ and v a = - a 

P P 

be the co-ordinates of C, the common centre of all the 
curves (10). This period depends not only on the 
co-efficients X, j,i, and /?, but also on the constant H ; 
its value is zero when « = u and v a = v a , and infinite 
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when H — 0. It has already been seen that the mean 
values of u and v are equal to u$ and v . 



CO 



u.dt = u ~L and — 
p a> 



v.di = v = - (11) 



This is a purely imaginary periodicity, lor equations 
(9) are not true to a first approximation, and the vital 
co-efficients are not constant, but vary slowly. Actually 
this imaginary periodicity corresponds to a real pheno- 
menon, to the rather irregular but periodomorphic 
fluctuations through which the system (u, v) passes 
and which might be confounded with periodic oscilla- 
tions when the observations are not spread over a 
sufficiently great interval of time. The curve repre- 
senting (u, v) may have cusps, which correspond to 
critical occasions, and limits, which correspond to 
relatively stationary conditions. In any case, it is 
seen that, apart from external factors, which may 
disturb the equilibrium (volcanic activity, glaciation, 
etc.) the actual structure of the organic world must 
produce a longer or shorter series of oscillations. .The 
disappearance of these oscillations, owing to the adapt- 
ation of the organic world to a new set of circumstances, 
leads to a new stationary state, to a new stability which 
is never permanent. 

Consider now equations (2") and (6), which may be 
written more simply, omitting the indices : 

y' = e — yzy, z' = — a.z + yzy (12) 

The action proceeds towards a stable stationary con- 
dition 

«■ e . .• 

y ^~ *„=.- 13 

Y a 
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The word 'stable' must be interpreted in the sense 
that the state (13) can be maintained for a longer or 
shorter time, but it is certain that the co-efficients of 
equation (12) vary and that the state (13) varies with 
them. 

Suppose that until a moment t a the volcanic effect 
is weak, so that* 

e/? < 4a 2 

and that the stationary condition of equation (13) is 
almost reached. Suppose further that from this moment 
until a moment t t a period of volcanic activity occurs 
with an average effect E far greater than e, 

E >e. 

Equilibrium is disturbed and the system (y, z) tends 
towards a new equilibrium state 

Y Q == - = y Z a ~ ~ > z 
y a 

Suppose that this state is almost reached at the time 
t v It is apparent that the new ultimate mass of oceanic 
carbon dioxide, Y , is equal to the original mass y , 
but the total mass of marine chalk-producing organisms 
has increased. 

Assume that after this recrudescence of volcanic 
activity, normal relations are re-established, and that 
for t > t x the effect of carbon dioxide is again equal to 
e. The state (y, z) tends to the original condition of 
equilibrium (13), describing decreasing spirals. This 
process is at first almost periodic, and the masses of 
carbon dioxide contained in the ocean and the atmo- 
sphere pass through a series of maxima and minima, 

* -^ «o is much less than, ^- = is much more than. 
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The same is true of the total mass of marine chalk- 
producing organisms, and in consequence terrestrial 
organisms are experiencing a critical occasion. 

The co-efficients of equation (12) can be approxi- 
mately calculated. Volcanoes and other internal 
sources evolve on the average io 8 metric tons of carbon 
dioxide in a year. The total mass of living marine 
organisms is of the order of xo u metric tons and it 
may be assumed that the mass of chalk-makers is of 
the order of io 13 metric tons. On the other hand, the 
mass of carbon dioxide dissolved in the ocean is io 14 
metric tons. Thus we get : 

e =. io 8 , y =. io u , z — xo ls , «■ = *o~ 5 , y — xo~ 19 

These values practically verify the inequality 

ejS < 4« 2 

necessary for the establishment of periodomorphic 
oscillations after a disturbance of equilibrium. 

Are any traces of such oscillations to he found? 
Perhaps the series of chalk deposits may be explained 
by this mechanism. It seems that the marine deposits 
of the Black Sea show a periodicity of the order of two 
to three thousand years. I have tried to relate these 
oscillations to the quaternary glaciations, but another 
periodic mechanism affords a better explanation of all 
the phenomena accompanying glaciation, as well as of 
the glaciations themselves. 

We pass on to the consideration of equations (1), 
(2'), and (5). Suppose that the time is one of relative 
stability and that the stationary condition (13) is 
almost reached as far as pelagic organisms and dissolved 
carbon dioxide are concerned. Integrate equations 
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(1), (2'), and (5) in respect of t, from t to t -f- o>, putting 
throughout 

u, = ^55 ^ _ £ja 

We find that the average increments of x, y, and s are 
equal to : 

x[t -f- to) — xU) 11 x . .A , . 

= ~ «U a + (««— «1«) « - %B (14) 

= «3s| + (««a- a at) « + E'52 ( I 5) 

The ultimate result, or at least the actual tendency 
of this process depends on the vital co-efficients. If 
these co-efficients are in agreement with certain in- 
equalities that are easy to establish and which lend 
themselves to simple biological interpretations, the 
atmospheric gases will be in a. state of rhythmical 
increase or decrease, according to the periodic oscilla- 
tions. Suppose, for example, that 

«14— **! < o 

which means that plants evolve more oxygen than 
they use in respiration and also that a 13 and % are 
small enough to ensure a positive sign for equation 
(14). In this case the quantity of free oxygen will be 
continuously increasing. 
Suppose, again, that 

a 42—«24<° 

that is to say that plants absorb more carbon dioxide 
than they evolve in respiration, and also that « 32 and 
s',52 are small enough to ensure a negative sign for 
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equation (15). In this case atmospheric carbon dioxide 
will be continuously decreasing, save for the regulating 
action of the oceans, but this reserve itself will not be 
inexhaustible. 

Other consequences which follow from equations (14) 
to (16) will not be considered. 

3. The Nitrogen Cycle 

The above results may be completed by a study of 
the circulation of nitrogen. To simplify matters, sup- 
pose that the micro-organisms able to fix atmospheric 
nitrogen live symbiotically with plants : this will not 
appreciably affect the results and will diminish the 
number of variables and of equations. Let 

x — the total weight of free atmospheric oxygen. 
y — the total weight of carbon dioxide in the atmo- 
sphere and the ocean. 

z ~ the total weight of free nitrogen. 

u = the total weight of these substances in animals. 

v — the total weight of these substances in plants. 

s = the total weight of these substances in the earth's 
crust. 

The relation between these variables may be repre- 
sented by the accompanying diagram (Pig. 6) or by the 
differential equations : 



*' = --»l4« + («6t~«lli)* 






(17) 


y' — a 42 M — (a 25 — a^)v 






(18) 


z' = — « aB u — e 






(19) 


u' = — (a 4a + a 4a — « M )w 4- jiwo 






(20) 


V' = (<x 15 — a fil + « 36 — a B2 + a 36 - 


aee) 1 " - 


- fiuv 


(21) 


s' = a^M H- («, 5 — a 6B )v + £ 






(22) 
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the formation of which is similar to that of the equa- 
tions in the preceding section. Suppose, as before, that 
in the absence of plants the balance of animal life will 
be diminishing, giving, 

Similarly, 



'•i* ■ 



P = «15 ~ «51 ' 



mosphirlc i 
lia'nqan j 




Fig. 6. — Atmospheric Circulation 

In these circumstances, equations (17) to (22) give 
results that do not differ from those of the previous 
paragraph. The variable z is always decreasing. Equa- 
tions (20) and (21) have again periodic solutions. 
Adding the equations of the atmospheric gases, we get ; 
x' + y' 4- z' 
= «(a 42 — a 14 ) — (a 1E — a 81 + a 26 ~ a 62 + a 8B )u — e 

Integrating this equation in respect of t from t to 
t + w, we get the average rate of disappearance of the 
atmosphere 

x[t + 0) +y{t + m) + z(t + o>)~ x{t) -y(t) - z{t) 



CO 

Aa 56 + /ia M 



£ < O 
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Hence the atmosphere, and the nitrogen in it, are 
continuously disappearing. The maximum speed of 
this disappearance can be estimated by simple con- 
siderations. The total mass of atmospheric nitrogen 
is of the order of 4 x lo ls metric tons. The earth's 
crust, at least in its upper layers, contains an amount 
of the same order of magnitude. Suppose that this 
combined nitrogen be of atmospheric origin, and that 
the process has so far lasted for 2 X 10 s years, then 
the annual loss of atmospheric nitrogen does not exceed 
2 X 10 8 metric tons, and by this computation the 
atmospheric nitrogen will not be exhausted before the 
passing of two thousand million years. 



CHAPTER IV 

GROWTH OF A POPULATION : THE 
LOGISTIC LAW 

In this chapter the growth of a.n isolated population 
will be considered. Such a population closely follows 
the law of Verhulst, a fact that is confirmed both by 
observation and experiment. 



i. Differential Equation for an Isolated 
Population 

In the first chapter we arrived at the differential 
equation 

p>(t)^ep(t)~hfi*(t)+I-E (I) 

where e is the co-efficient of multiplication, h the limit- 
ing co-efficient, I the co-efficient of immigration and 
E the co-efficient of emigration. The co-efficient of 
multiplication s is the difference between fertility n 
and mortality m. 

e = n — m 

For an isolated population with neither emigration 
nor immigration, equation (i) becomes 

p'{i)^sP(t)~hp*(t) (2) 

and in this form is known as the logistic, equation. 
The co-efficient h is generally small compared to e, 
and when p(t) is small enough the second term of the 

56 
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right-hand side may be neglected, which gives the 
Malthusian form of the population law 

f(t) = ep(t) ( 3 ) 

whence p{t) = p[o)e^ (4) 

One could give many biological examples of exponen- 
tial growth according to this law ; it happens whenever 
the necessities of life are sufficiently plentiful to abolish 
all competition between individuals. This occurs in 
the colonization of virgin territory so large that for a 
long time each colonist cannot 'see his neighbour's 
smoke.' As the population multiplies a sense of being 
cramped arises and the limiting term must be taken 
into consideration. 

Example. The following is an example which shows 
that in optimum conditions fertility may remain con- 
stant for a long time. It deals with a small rodent, 
Microkis armlis Pall, which multiplies with great 
rapidity. The following table gives the increase of the 
descendants of a single pair, mortality being neglected. 



Months 





2 


+ 


6 


8 


10 


12 


14 


p observed . 
p calculated 


2 

2 


5 
4-5 


16 
10 


20 
22 


40 
50 


109 
112 


200 
252 


283 

282 



Calculated values are based on the formula p — 2e ,4 * 
The agreement is as close as possible, for pregnancy in 
these animals lasts three weeks, and the difference in 
time between the taking of a census and the birth of 
the young is often considerable. 

Mote on the Limiting Co-efficient. Another interpreta- 
tion of the limiting term, differing from that given in 
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Chapter I, is possible. Numerical growth of a popu- 
lation may be impeded by competition between 
individuals, and the effect of this may be taken 
to be proportional to the number of encounters per 
unit of time. This number is proportional to p^it), 
from which is derived the form of the limiting term 
— hp 2 (t) in equation (2) . It may be noted that this form 
of the limiting term is too simple and takes no account 
of external limiting factors which have nothing to do 
with either encounters or competition. It is the whole 
of these factors that give its form to this term, and 
experience amply confirms the legitimacy of using it 
in the logistic law. 

Chemical Interpretation. Equation {2) may be given 
a slightly different form — 

dp e 

-~ = hp (cu — p) where ca = j {5) 

In this form the logistic equation expresses the progress 
of a reaction of which the velocity is proportional to 
the mass p of the products and to the mass (<u — p) 
of the reactants. This interpretation is meaningless 
as far as a population is concerned but it may be 
applied to the reactions in an organism. 

2. The Logistic Law 
Write equation {5) in the form 

codp It 1 \ , 

p[a>~p) \p to-pj * 
Integrating, we get, when p < a 

^logfc^ (6) 
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and when p > w 

d =»log$fc=$ (60 

(p~a>)p 

These formulae give the logistic law of population- 
growth : 

p ~ hTi^W^ (7) 

The limit of p as if -*■ 30 is equal to w whatever the 
initial value of p ; when p < co, p increases to a>, 
and when p > co, p decreases to co. 

Wlien p < — > curve (7) passes through a point of 
inflexion at 



e p, 







where the value of p is — - and when — <p a the maxi- 

3 2 

rauxn value does not exist. 

The fundamental conclusion is therefore : an isolated 
population tends towards a limit which only depends on 
the vital co-efficients; it is independent of the initial value 
ofp . This conclusion is only valid if the whole of the 
factors determining the vital co-efficients remain un- 
changed. These co-efficients depend on the neighbour- 
hood occupied, on its resources, on its stage of mechan- 
ical development, and on the social organization of the 
population. 

The United States. It is interesting to find that the 
logistic law 

, 197.373,000 (8) 
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closely expresses the population of the United States. 
This is shown by Fig, 7 and the following table : 



Year 


p Observed 


p Calculated 


Yeah 
1870 


^I OnSBRVSD 


/> Calculated 


*79° 


3'9 X 10' 


3-9 x 10" 


387 X JO* 


39-4 x ro» 


18110 


5-3 


5'3 


18U0 


SO- 3 


5o-» 


1810 


7-1 


7.3 


1890 


63-i 


6i-8 


X820 


9-0 


g-8 


1900 


76-r 


76-9 


1830 


is-9 


zyx 


S9» 


93-3 


92-0 


1840 


I7-I 


i?-5 


1920 


106-5 


109-4 


1850 


23-3 


23-i 


*930 


133-3 


133-9 


i860 


31-5 


311-4. 









The agreement is almost baffling. Since 1790 the 
United States have waged five major wars, conquered 
new territory, colonized vast spaces, developed # a 



zoo 




1790 1800 



1850 



Fig, 7. — Example oj? a Logistic Population 

Growth of Population in the United States. 

formidable industrial organization, and have received 
and absorbed a mass of emigrants from all parts of 
the world, yet, in spite of all this, the logistic law closely 
expresses the growth of the population. It is seen that 
the point of inflexion at which accelerating growth 
gave place to diminishing growth was only passed in 
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1914 and that a state of saturation is a long way ahead. 
At the moment the population density is only 20 persons 
per square kilometre, which is still far below the 
European saturation-point. 

Prance and Belgium. As far as European countries 
are concerned, Verhulst in 1845 calculated the condi- 
tions of saturation for France and Belgium. For 
Belgium he prophesied a maximum population of 
6,600,000, and for France about forty millions. The 
present population of Belgium is (1930) 8,092,000, 
corresponding to 266 persons per square kilometre, 
where the forecasts of Verhulst showed 224. The 
difference can be explained by the astonishing rise of 
industry in Belgium, and by the acquisition of the 
Congo which secured for the country sufficient addi- 
tional wealth to support the extra population. 

The present population of France, exclusive of 
foreigners, is in remarkable agreement with Verhulst's 
forecast. The theoretical maximum is almost reached, 
which explains the low birth-rate and other social 
phenomena which perturb public opinion. This state 
of affairs is only a temporary disturbance which will 
quickly pass with the growing value of the colonies 
and the decentralization of the metropolis. It is enough, 
further, to compare the rapid growth of the French 
population of Canada with the very slow growth of 
population in France to understand the biological and 
economic nature of these phenomena. 

Let us return to the population of the United States 
and compare the concentration in the New World 
(20 persons per square kilometre) with that already 
reached in Belgium (266 persons per square kilometre). 
These figures show that in America the economic and 
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psychological factors characteristic of the early days 
of cultivation still predominate, while in Belgium every 
corner of soil is utilized. It is to this state of things 
and to this spirit that the co-efficients 

e = -03134, co — 197,273,000, h ~ i-6 x 10" IO 

correspond in equation (8). 

A time will come when the more complete utilization 
of the soil will be as inevitable in America as in Europe, 
and then the modified vital co-efficients will give an 
extension of the population curve incompatible with 
the arc of 1790 to 1930 determined by equation (8). 

It should be noted that nothing is as deceptive as 
the beautiful agreement between the observed and 
calculated values. A series of observed values is, in 
fact, equivalent to a narrow band rather than to a 
curve, and in this band can be traced a number of 
curves, corresponding very well with the conditions of 
the problem within the period of the observations, but 
showing real divergence beyond these limits. 

Growth of a Bacterial Colony. After this sociological 
example a biological one taken from Lotka may be 
given. It deals with the growth of a colony of bacteria 
[B. dendroides) observed by G. Thornton. In this 
example we have : 



Aon, Days 





1 

2-Oz 


2 


3 


+ 


5 


Area 1 observed 
occupied ! 
in sq. cms. \ calculated 


•1$ 


I3'53 
I3'08 


3«'30 
37-05 


47'50 
47'39 


49-40 

4?'0* 

-. .. ... 



The calculated values were obtained from the formula 



P = 



■2524 



- 3'I28< 



+ -005X25 



It is seen that the agreement is as close as possible. 
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Bacillus coli. Regnier has investigated a number of 
cultures of Bacillus coli varying at the original inocula- 
tion from 500 to 500,000 bacteria per c.c. In these 
experiments the medium was not renewed and the 
toxic action of accumulating metabolic products ulti- 
mately modified the growth-curve, but its beginning 
agreed closely in all experiments with the Malthusian 
law (4), with e — i«i2, whatever might be the original 
number of germs. We shall return to these experi- 
ments in discussing the toxic action of the products 
of metabolism. 

3. Experimental Verification of the 
Logistic Law 

The example of the United States population shows 
a remarkable agreement between the observations and 
the logistic law. This agreement is not, however, suffi- 
cient, for it is always possible to express a series of 
numbers that increase and asymptotically approach 
a limit by a curve having the same properties. The 
logistic equation may be replaced by another more 
general equation — 

*' = *&) ="I>{P)~VU>) (9) 

where ^(p) is the increase of population, strictly speak- 
ing, and f(p) is the total of limiting factors. It is 
admitted that the action of these factors is scarcely 
perceptible at the beginning of the process, and that 
when p is sufficiently great the order of growth of 
f(f) is higher than that of <j>{p). We then have the 
following inequalities : 

<!>{£) > viP) (0 < P < ">) 

4{P) < y(P) (<» < P) 
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and therefore 

p'(a>) — (j>(co) — f(w) ~ 'F{co) ~ o 
F'(«>) = f (a>) - y/(a) < 0. 

Let F(j>) be expanded by Taylor's theorem to the 
neighbourhood of the point w. li (<D~p) is small we 
set 

V(-P) («-£)F» 

p' = -(a>-p)F'(w) 

and therefore 

p~co~~ Ae x * Mt (io) 

In the same conditions the logistic equation gives 

la> — p n )co . , v 

p~a~ J — ~~ e~ '' (n) 

Po 

Suppose, on the other hand, that (f>(p) is expanded 
in the same way and p is very small 

where <£ — ° an d ( iV > ° an( i f(P) disappears when 
p — o. In these conditions, when p is small, equation 
(9) becomes 

P' = MP 
and we again get the Malthusian law 

P = P^'\ (12) 

Thus in these broad hypotheses on the otherwise 
arbitrary functions <J> and y, the initial and final parts 
of the growth of a population follow exponential laws. 
Junction between these phases is very easy, and in 
consequence the S-form of the logistic curve may 
belong to a number of theoretical curves, all satisfying 
equation (9). 
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Admittedly, in the case of the logistic law the two 
exponents are identical : 

f(o)=~F'(a>), 

but this criterion is not certain, for the calculation of 
exponents is often difficult. 

Suppose now that by some means we succeed in 
modifying the co-efficient of increase, for example, 
by regularly removing a part of the population pro- 
portional to p. Equation (9) becomes 

P' = MP) — m P — w(P) 

This complementary mortality modifies co, but it is 
not possible to say in what direction, for when the 
functions <j> and y are unknown, the solution II of the 
equation 

4>iP) — m P ~ wiP) ~ ° 

"is not related to that of the equation 

4>(P) -y>(p) = o 

by any simple expression. On the contrary, in the case 
of the logistic equation we have 

p' = (e — m)p — hp % 

and the modified maximum of the population becomes 

_ b — m E — m , . 

n = __— « co (13 

h e 

In consequence if at the end of this experiment the 
new maximum II is found to be connected with the 
former maximum co by expression (13), there is a very 
cogent argument in favour of the logistic law. 

This argument occurs in the numerical results of 
the experiments of Gause on cultures of Paramecium. 
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By taking special precautionary measures to eliminate 
metabolic products, and by regularly renewing the 
medium and the food supply, Gause obtained popula- 
tions which reached the stationary state and maintained 
it for a relatively long time. For the purpose of 
numerical control Gause made periodic abstractions of 
the populated medium. In three series of experiments 
the fraction removed was -x, -2, and -3 of the medium 
and hence of the population.. If the logistic law is 
exact, the maxima in these series should be equal to 

ly* = e — -i ; l\h — e — -2 ; TL 3 h = « — -3. (14) 

Eliminating h and e from these three equations we 
get the following relation between the three maxima 

sn 2 - n a + n 3 . (15) 

In the numerical results obtained by Gause there 
are the following average values : 





n, 


n. 


n. 


Paramecium aurelia . 
Paramecium bursaria 


475 
340 


375 
240 


305 

190 



which closely verify equation (15). This confirmation 
of the logistic law is much more significant than the 
close agreement between the theoretical curve (2) and 
the empirical curve (8) , The limiting term in equation 
(2) is as well represented as the term of multiplication ep. 



CHAPTER V 
POPULATION AND ENVIRONMENT 

In this chapter the growth of a population will be 
studied while very varied factors will be taken into 
account: food, available area, immigration, emigra- 
tion, epidemics, complexity, residual actions, etc. 

i. Non-isolated Population 
Consider the equation 

f{t) = ep(t)-~hp'(t) + I~-E (I) 

Several causes can be foreseen that attract or repel 
immigrants. So far as human population is concerned, 
there are classes of immigrants who are attracted by 
the rise of a country ; there are on the contrary others 
whose wish to expatriate themselves increases with the 
growth of the population and who seek a virgin land. 
In an animal population, great migrations are almost 
always provoked by lack of food. Predators especially 
are mobile and seek the places where food is abundant. 
These two streams may alternate, and create local 
fluctuations that may be periodic or periodomorphic. 

Express the total of these two streams by a poly- 
nomial 

I _ E = „ + Vl p + V 2 (2) 

in which the co-efficients v , v x , and v 2 are positive or 
negative constants. Equation (i) then becomes 

p' - v + (e + vjp -{h- v$>* (3) 
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The terms of which expression (2) is composed have 
not been, invented ad hoc. The term v expresses the 
constant fraction of the stream across the frontier, the 
term v x p the positive or negative attraction of a popu- 
lation or the regular -withdrawals from bacterial cul- 
tures. The term a,^ 2 refers to sudden migrations of 
predators, etc, A few special cases will be considered. 
Unfavourable Case, It may happen that 
v Q < 0, • £ -j- v 1 <o, h~v s >o ; 

p'(t) is then always negative and p(t) diminishes and 
approaches zero. 
Favourable Case. Suppose that 

v > 0, h — v a > 0. 

The polynomial v + (e + v x )p — (h — v^p & has a nega- 
tive value with which we are not concerned, and a 
positive value which may be denoted by co; then, 
when p a < co, the derivative ■p'if) is always pbsitive, 
the population p[t) increases and tends towards the 
limit co; when p > to, p'(t) is always negative and 
p(t) decreases towards the same limit 10. 
Unlimited Immigration. Let 

v > 0, e + v t > 0, h — v 2 <o; 

in this case the population theoretically increases con- 
tinually and tends towards infinity. In fact the term 
Vgp* only appears in exceptional circumstances and 
immigration never lasts long. 

Opposing Factors. Suppose that the constant stream 
v Q is negative, but the attraction of the neighbourhood 
is considerable, so that 

e + v x > ; 

in this case two results may appear ; 
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(«) [e + v 1 ) z + 4.v (h — z/ a ) < ; 

^'(£) is always negative and p(t) tends to zero. 

the polynomial v + (e -|- v t )p — hfi 2 has two positive 
limits, co and co 2 > m x ; when j!> < oi v p'(t) is always 
negative and pit) tends to zero; when co x < p < co 2 , 
p'(t) is always positive and pit) increases towards co 2 ; 
when p > co 2 , p'(t) is always negative and p(t) dimin- 
ishes towards co 3 . Hence in this case the final end of 
the population depends on the original value p Q , while 
the two ultimate limiting values and a z do not depend 
on p . As for w a , it is an unstable stationary state. 

We shall not consider other possibilities which do 
not differ essentially from the four cases studied. 

2. Periodic Fluctuations of a Homogeneous 
Population 

A homogeneous population under the influence of 
intermittent factors may pass through periodic fluctua- 
tions which the classical form of the logistic equation 
does not provide for. A population that has reached a 
certain density may become the victim of an epidemic 
which only stops when the density has fallen to a limit 
which renders the propagation of the epidemic im- 
possible. Other causes may produce the same effect, 
for instance, alternation of the attraction of great cities 
with a return to the country, 

We shall consider the evolution of a population of 
this kind, assuming that it is originally governed by 
the logistic law 

f =^ B p — hp* (4) 
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and that the epidemic appears as soon as ft reaches a 
certain value Q 

At this stage the vital co-efficients become 

E < e, li >h 

and equation (4) is replaced by 

p'^Ep-Hp*. (5) 

E 
Suppose that Q > — :• The population decreases, 

A point is reached when the density falls below a 
certain limit which puts an end to the spread of the 
epidemic, Let q be the corresponding value of p : 

E 
a > — 

At this moment the epidemic ceases, and the popu- 
lation again begins to grow, again following equation 
(4) until the incidence of a fresh epidemic. In this way 
there are periodic oscillations of p between q and Q. 
The period Q, of these oscillations may be calculated. 
In the first part of the cycle, where p = q, we get, 
using formula 6, from Chapter IV : 

e q{e~hp) 
Hence the time taken by the first part of the cycle is 

T ~t ~t - I i oc Q (8 ~ ;t g), (6\ 

Ll - h h ~ s l% qi^-lQ) (6) 
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Similarly we get for the second part of the cycle, 
E being positive, 

which gives for the final point 

T^^-^glog^g^, (7) 

and the whole period is equal to 

a S (e-hQ)q^E g Q(H<?-E) U 

This representation of a two-phased biological process 
is of course not strictly accurate and needs some revision, 
but it presents a reasonably correct picture of the 
mechanism of the formation of periods of relaxation. 
This type of periodicity is often, though incorrectly, 
contrasted with classical oscillations, for example with 
elastic oscillations. The latter are, in fact, only a 
special case of the former. Oscillations of relaxation 
are much more frequent in nature than classical oscil- 
lations and to them must be largely attributed the 
maintenance of order and the regulation of states of 
chaos. 

The sudden appearance of hordes of small rodents 
and their no less rapid disappearance may be explained 
in this way. This phenomenon shows an irregular 
periodicity with a period in some cases of 9 to n years, 
in others of 3-5 to 4 years. Some authors attribute this 
periodicity to the action of sunspots, but there is no 
need for this, for purely terrestrial and even local factors 
are sufficient to explain it. Consider these factors. 
First, the observations of Kalabuchow, Sewertzoff, 
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etc., have demonstrated the importance of the den- 
sity of an animal population in the appearance, 
spreading, and disappearance of plagues. Further, a 
local increase of density attracts predators in large 
numbers. These two factors will in two or three weeks 
succeed in destroying 97-98 per cent, of a population 
of small rodents, and the density falls to the level at 
which the disease cannot spread. The population, 
reduced to 2 per cent, of its maximum, finds it refuges 
sufficient and its food abundant; the new conditions 
favour a slow increase in numbers, which in a few years 
re-establishes a state of affairs favourable to another 
wave of disease. Of course these conditions are never 
exactly reproduced and the fluctuations are not periodic 
but periodomorphic. 

The speed of reproduction in mice is so great that 
in the first part of the cycle limiting factors may be 
neglected, putting h — 0. In the second part of the 
cycle, on the contrary, E is very small in comparison 
with H, and may be neglected. In this way formulae 
6-8 may be given a simpler form : 



T Z 1 Q 
1^-log- 

T - £=1 

o^i og Q + y? 



(9) 



These equations enable us to calculate the vital 
co-efficients 



1 1 Q 

£ = — log -, 

h 9. 



H=a 



Q?T 2 



(10) 
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when Q, p, T 1( and T a are known. In the case of mice 

9,0 
T x ~ 4 years, s „ 5Q> i og - ^ 4 

? ? 

whence e ~ 1. 

This value of s corresponds very well with the rate 
of multiplication of mice in natural circumstances, 



3. Isolated Population Poisoned by its own 
Metabolic Products : Residual Action 

It often happens that the existing state of a system 
depends not only on the initial values of the variables, 
but on its whole history. In plant associations one 
species often "tills the soil " for another by the accumu- 
lation of nitrogen compounds. From the mathematical 
point of view this fertilization is expressed by comple- 
mentary integral terms in the biological equations. In 
an animal population the accumulation of metabolic 
products may seriously inconvenience the whole popu- 
lation or one or other of its members, and this is ex- 
pressed by a fall in birth-rate and an increase in 
mortality. If it be assumed, as Vdlterra has done, that 
the total toxic action on birth-rate and death-rate is 
expressible by an integral term in the logistic equation, 
we get the following integro-differential equation : 



fi'(t} = ap(t)-hp*(t)~cp(f) 



K{t~T)p[r)dr, (10) 



The function of residual action K{t — t) is a diminish- 
ing one. Generally we may write 

K(o) » 1. 
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It may happen that at the end of a certain time T the 
residual action ceases to exert an effect, which gives 
K(t) = o t > T 

P'(t) = ep[i) - hp*{t) - cp[t) f Kit - r)p[r)d r . (il) 

J (- T 

On the other hand, it may happen that residual 
action lessens but persists indefinitely, which gives 
lim K{1) — a < i. 
t-> a 

Simplification. Equation (10) cannot be integrated 
in its general form ; to simplify it put 

K(*-t) = i (12) 

when it becomes 

P'(t}~ep(f)-hj>*(t)-cp®l>(f> (13) 

where P(/) is the integral 

V(t) ~ j' o P(r)dr (I 4 ) 

Hence p{t) = P(<) £'(() = V"{t) 

and equation (13) becomes 

P" = eP'~AP' il -cPP'. . (15) 

Integration. This equation can be integrated. Multi- 
plying it by e hV we may write 

it l J dtl \h^ h* hi 

which after being integrated gives 
s c 



A 



h A a 



= J* 



r h h WJ 



(16) 
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From this we get the solution of equation (13) in 
parametric form : 

/ e c \ ,-t, e c c_ 

'-U> (I8) 

These equations give the population ft and the time 
t as functions of a parameter P. 

Discussion o£ the Solution. The population p cannot 
be negative. Hence the right-hand side of equation (17) 
vanishes for a finite value of P which we will call P . 

CO 

This value is reached when t — o> , which implies that 
the population disappears at the end of a very long time. 
Moreover, supposing that 

h<- h = a> 

we find that p increases with sufficiently small values 
of t. Hence there exists a maximum value for the popu- 
lation which is -reached in a finite time. Let p m be this 
maximum and t vl the time. It can easily be shown that 
e c , c -|- eh — pjfi s 

*- = »-*-» log — « <* (I9) 

Substituting this value for P in equation (18) we get 
i m . After having passed this maximum, ft tends to- 
wards zero. There is always one point of inflexion in 
the interval (t m , :»), In the case where 



3 , c x / a , a c , c 

4 4/i 2 4^ A 2 A 4 
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there is still one in the interval (o, t m ). Finally, when 
p Q > (a, there is no maximum, p'{t) is always negative 
and the population tends towards zero. 

Hence, in every case the accumulation of metabolic 
products is fatal to a population. 

Approximate Formulae. Expressions (17) and (18) 
are not in a convenient form, and may therefore be 
replaced by approximate formulae which are not 
precisely accurate, but which adequately represent 
the empirical curves. Suppose that c is small in 
relation to //. Then, when t < t m we may calculate 
p by the logistic formula, replacing co by p m . This 
gives 

when t > t m , equation (18) may be written thus : 

fV 



ds 



eh + c— chs— ce"~'' ls ~ p " 



(2!) 



using equation (20). We can replace e~ h[i '~" p " i) by its 
expansion 



fl ~/.(s-P m ) „_ 


i-ft( S _P ro )+: 


and, moreover, 


C 


Equation (21) now becomes 


t~t m ~- 


Pm ds 

y 2 -« a 



2 



where 2p m = cy % . 
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On integration we get 



t~t m = - log ■' 



whence 



c °y-P + P„ 



V-¥ m = yth C -tlM (22) 

# = -rn (<>o (23) 



It may be noted that when A = o this formula 
exactly represents the growth of a population p from 
the first moment. 

R^gnier's Experimental Populations. These results 
may be applied to J. Regnier's observations, to which 
reference was made in the last chapter. He was working 
with cultures of Bacillus coli and Staphylococcus, in 
a nutrient medium that was not renewed. The popu- 
lation grew, passed a maximum and finally disappeared. 
To explain the decrease and final disappearance there 
is a choice of three hypotheses : (i) lack of food, (z) lack 
of space, (3) poisoning of the medium by metabolic 
products. In different ways these three factors lead 
to the same result — total disappearance of the popula- 
tion. The third hypothesis seems more probable than 
the others, and we are here concerned with its verifica- 
tion. Unfortunately the counting of bacteria is an 
operation involving such gross errors that the results 
can never be very certain. 

The co-efficient for B. coli has already been deter- 
mined — e — 1-14: the others must now be found. 
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The following table gives the results of calculation and 
experiment. 



Age or 




P 


Population 


w Millions 


in Houks 


Observed 


Calculated 





•5 


■5 


1 


7'3 


4'9 


4 


39-8 


47'3 


6 


380 


417 


8 


80?. 


2126 


24 


2178 


39So 


3" 


3992 


3992 


48 


3574 


3«3 


96 


2050 


23.54 


192 


354 


3^7 



We gat p m = 3992, * TO = 30 hours. 

Equation (23) gives, when i = 192 



Ch 8lC = JW = 3-36 



whence 



334 



•023 



Equation (19) then gives //■ < -0003. 

The calculated values of ft agree very well with the 
observed values at the beginning and at the end of 
the process but not very closely between 8 and 24 hours. 

The same thing is found for Staphylococcus : 



Ago (liours) 





2 4 


P. ote. (millions) 


■5 


G'9 :)=■+ 


I', calo. (millions) 


'5 


4-1 33'« 



fi 8 34 30 4.8 96 19a 
332 628 3049 536O 4029 3385 390 
260 1569 5300 336C 5071 '4610 290 



The second line gives the average of four experiments 
carried out in identical conditions. The divergence 
about the averages is very great. The third line gives 
results calculated with the following parameters 

Pm — 5366, t m = 30, £ sa 1-05, C = -0265, h < -0002. 
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In spite of the imperfect nature of these results, 
some interesting facts may be stated. First, the almost 
equal values of the co-efficients of intoxication c for 
the two species shows that a chemical factor is involved 
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Fig. 8. — Bacterial Population 
'Staphylococcus) 

to which the two species react almost identically. The 
limiting co-efficients are very small in comparison 
with the co-efficients of intoxication, c, and h for B. coli 
is higher than for Staphylococcus. The co-efficients of 
multiplication £ are almost identical for the two species. 
Finally for both species the interval 8 to 24 hours gives 
a great deviation from the logistic law, of the same sign 
in both cases. This deviation must be ascribed to the 
effect of the two other factors mentioned. 



CHAPTER VI 

POPULATION AND ENVIRONMENT 

{continued) 

Differential equations of the logistic type or of the 
more general type (Equation 9, Chap. IV) have one 
serious defect : the growth of a population is considered 
without taking into account the sources of food-supply 
or the area occupied. These factors have a varied and 
often unexpected action. We shall consider some 
important cases. 

1. One Population Feeding on Another 

Imagine a species of animal feeding on plants or on 
another animal. Let u be the total number or biomass 
of the predators and v the biomass of the prey in such a 
population. This problem has already been considered 
for populations that are not dense. In the more general 
case the relations between prey and predator can be 
represented by the system : 

v' = SjV — h n v 2 — h^uv 1 , . 

u' = — e a « -j- h n uv ■— h^u 2 J * ' 

A full consideration of these equations must be kept 
for Chapter IX ; we are here confined to a few remarks 
about these equations and their solutions, especially 
from the point of view of the predators. 

First, the form of equations (1) implies a sufficiency 
of food for the prey, and when this is not so a third 

80 
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equation must be introduced (see Chapter IX). 
Secondly, the heterogeneous structure of two popula- 
tions, as for example in the case of an animal that 
undergoes metamorphosis, makes considerable caution 
necessary in drawing conclusions. In some cases it is 
enough to modify the vital co-efficients ; in others it 
may happen that the mode of feeding is different for 
larvae and adults, and that a scarcity of food that is 
fatal to adults may be of no consequence to the larvae. 
The species has then to be divided into two or more 
groups. 

Two types of solution are possible for equations (i), 
one of which is favourable to the predators while the 
other is not. In the first case we have 

e z h n <« 1 & ai ; 

an inequality which may be interpreted in several 
ways. In any case it is favourable to the predators 
and is of no disadvantage to the prey ; the possibility 
of a final equilibrium results from it. If the limiting 
co-efficients h n and h w are small enough, the system 
tends towards a limiting state, passing through 
decreasing fluctuations which sometimes bear the 
appearance of periodicity. 
In the second case the alternative inequality 

is satisfied and is wholly unfavourable to the predators ; 
at the limit 

Km u — o Hm v = -~^ 

t -9- DO * -*■ DO 711 

A peculiar paradox may be noticed in this case. 



82 MATHEMATICAL BIOLOGY 

It may happen that at the beginning the food-supply 
v is present in such quantity that the inequality 

«'<) = «o(— £ 2 + ^21 U — ^22 M o) > ° 

is satisfied; « then begins to increase, passes a maxi- 
mum and falls to zero. To an observer who sees only 
the predators without noticing the prey, such a case 
of a population passing a maximum to total extinction 
in the limit would be very difficult to interpret. 

2. A Population Feeding on Abundant Food 
Imagine that the food-supply v is so abundant that 

the loss of h 12 uv units in unit time is unappreciable ; 

the first of equations (i) then takes the logistic form 

v' = v(e x ~ h n v) ; 
whence 

«W=— r -r^-r^Tt where "i-p- (2) 

If the increase of the species u begins at the moment 
at which that of v is almost finished, v may be replaced 
by «! in the second of equations (i) and a logistic 
equation is again obtained for the predators : 

m' = w(E a — h^u), 

where the co-efficient of multiplication E 2 is equal to 

E 2 = h 21 co 1 — e 2 . 

If E 2 is positive the increase of the predators takes 
place normally. This example shows the complicated 
structure of the co-efficient of multiplication of a species 
and its dependence on external circumstance. 

If, on the other hand, the increase of v is not 
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complete expression (2) must be introduced into the 
second of equations (1). The solution of this equation 



L fii 



e~ 



».* 



I -f- A 22 M 



J L Ei 



(3) 



e~ e ' s ds 



is more comphcated than the logistic law, but when 

^2 11 "^ ^1 21 
it has the same character and u tends to the limit 

toi» = ■■■- * 
< -> x "•u'hz 

In the opposite case u tends to zero. 

3. A Population Feeding on a Limited Diet 
not Renewed 

In this case, which is very common in nature and 
in human society, system (1) takes the form 

u' = u{— s 2 -f h n v - h n u) \ , . 



v' ~~ k n uv 

These equations can easily be integrated. Put 

-t 
U = u{s)ds, 
J 

The second of equations (4) gives 
v = v a e~ !,,,u 
when the first can be written as 

u' = «(— f 2 + h 2t v e- h " v — h 2z u), 
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or 

(«' + A 23 « 2 )e"" u + ea«c*" u = A 21 u ms ( "»- Mu , 

which can be integrated giving 



u — F(u) = 



. £ 2 ^O 

% + t; 



n 22 — "is w aa 

* = f*- (6) 

Thus the problem is solved. The expression F(s) is 
positive when s = o and negative when s = do. It 
therefore vanishes for a value of s = U . This means 
that the predators finally tend to zero. Tf at the 
beginning w'(o) = o, that is to say if 

h 21 v > e z + h 2i u 

the number of predators increases at first, passes a 
maximum u m , and then falls to zero. It is easy to 
calculate a m , U m , and t m . A very interesting biological 
application of these formulae will appear later, 



4. Relation between Density of Population 
and Density of Food-supply 

In certain cases it is possible to study the relative 
density of a population in relation to its food-supply. 
It has several times been suggested that the study of 
density in relation to area should be replaced by that 
of density in relation to food, and it is undeniable that 
the suggestion sometimes gives a satisfactory result. 
We shall consider a biological example taken from 



POPULATION AND ENVIRONMENT 85 

Gause. This author studied the development of two 
uniform populations, one of which consisted of Para- 
mecium caudatum and the other of Paramecium aurelia. 
These protozoa are fed on bacteria in a saline solution 
regularly renewed and unfavourable to the multiplica- 
tion of bacteria. In some experiments the Paramecium 
only received a half-ration per day. In every case 
one-tenth of the population was periodically removed 
for checking purposes. If these circumstances are taken 
into account, the vital co-efficients are found to have 
the values shown in the following table. 





e 


Jl 


P. caudatum, ordinary ration 
P. caudatum, half-ration. 
P. aurelia, ordinary ration . 
P. aurelia, half-ration . 


I'2Q 

1-06 

1-28 
1-28 


■0089 

•0155 

■0265 
•0480 



The co-efftcient of multiplication is relatively stable, 
while the limiting co-efficient is inversely proportional to 
the daily ration. 

To explain this result, suppose that the density 

p — - of a population p in relation to a food-supply 

a satisfies the logistic equation 

P=p(e-h P ) {7) 

where the co-efficients e and h axe constant. When the 
food supply is maintained at a constant rate as in 
Gause' s experiments, equation (7) becomes 



p'~p{s~-p) 



(8) 



in perfect agreement with experiment. This result 
may be taken as yet another argument in favour of the 
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logistic law as the basis of the mathematical theory of 
population growth. 

It may be noted that there is no need for equation 
(7) in order to justify equation (8). The mutual pres- 
sure between the elements of a population becomes 
greater as food becomes scarcer. Hence it is quite 
natural that if a maintains a constant value the limiting 
co-efficient should be inversely proportional to a, while 
e does not vary. If a varies, equation (7) no longer 
fulfils the conditions of the problem and is inconsistent 
with equation (1). 



5. Population and Area: First Case 

The area occupied by a population is one of the most 
important limiting factors. In the case of experimental 
populations (Drosophila, Protozoa) the influence of 
the neighbourhood on birth-rate and mortality quickly 
makes itself felt. Each individual makes its immediate 
neighbourhood inconvenient for the others. After 
its death, until its complete disintegration, it still 
occupies a certain area, in this way hampering the 
increase of the survivors, even without taking into 
consideration its defiling of the medium. Two conse- 
quences may follow. In the first, this hampering may 
persist indefinitely. In the second, a dead individual 
is definitely disintegrated at the end of a certain time 
and the space it occupied becomes available once again. 
These two hypotheses finally lead to very different 
results. 

Let p be the density of a population, s the room 
occupied per unit of volume by living or dead indivi- 
duals, a the volume of a live individual and /J that of a 
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dead one. Suppose, first, that the birth-rate is propor- 
tional to p and to the room available, 1 — s ; secondly, 
that the death-rate is proportional to p and s ; thirdly, 
that the action of other limiting factors is not negligible 
and above all affects the birth-rate; fourthly, that /? 
is constant. These suppositions allow us at once to 
write the equations : 

p' = pn(i — s) — mps — hp 2, (9) 

** 
s = pa. -f- mp p(u)s(u)du. (10) 

J 

In equation (9) the term pn(i — s) corresponds to 
the birth-rate, the term mps to the death-rate, and the 
term hp z to all the unspecified limiting factors. In 
equation (10) pat. denotes the room occupied by the 
living and the integral term that occupied by corpses. 

Transformation o£ Equations. Differentiation of 
equation (10) gives 

s' — tup' + mfips 

or, including equation (9) 

s' — p[cm — s(a.n -\- xm — /Sot) — - h«.p] (n) 

Equations (9) and (11) form an integrable system. 
Eliminating t from these equations 

dp n— (m -f- n)s — hp 



ds an — hap — s{a.n + am — fim) 

and let z be a new variable defined as the quotient 

n— hp 

ks 

dp dz 

Then -f- = — z — s — . 

ds as 



(12) 
(13) 
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The numerator of equation (12) becomes 
n — (m + n)s — hp = s(hz — m~- n) 
while the denominator may be written 
an — hap — s(an -\- am — fim) = s(aJtz — an — can -\- ffm) 
Put 
F(*) = ahz 2 -\- z(h — aw — am -f- fttn) — (m -\- n) (14) 
Equation (12) becomes integrable 

ds ahz — an — am 4- Bin , , . 

7 = YW~ dz fl5) 

Expression (14) disappears for two values of z 
z-i < o < # 2 . 

The original values of p and s are connected by the 
evident relation 

S = ap 

It follows that in every case s and p satisfy the 
inequalities 

5 — ap > 0, s < 1 

Integration. Integration of equation (15) gives s as 
a function of z 

s = s e-* w (16) 

where 

8 ah£ — cm — aw 4- /?w 



#*) = 
We then get 



F(f) 



# (17) 



^ = --^8-*W (18) 

Equations (16) and (18) give the population p and 
the volume occupied s in the form of a parameter, 
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taking no account of time. To introduce this term, the 
equation 

must be integrated, which gives the time t as a function 
of the same parameter z 

C*. du 

(19) 



F(«) 



j — us a e~ *f» 



Equations (16), (18), and (19) completely solve the 
problem. Since the polynomial F appears in the de- 
nominator in the integral in (17) we must make sure 
that the critical values z ± and z 2 do not occur in the 
interval of integration, that is to say that the expression 
F always has the same sign and that z varies in the 
same direction. 

Effective Integration. Integrate the following. Put 
ahg — an — can -f- fim _ A x A 2 



and we get the following formulae 

ahz 1 — an — <xm -\- fim 1 + <xz 2 



A l = 



z<tkz 1 -\- h — <m — vm -\- flm a.(z 2 — gj) 



afe 2 — aw — aw + [5m I -f- a&i 



2afe 2 -+■ 7i — aw — am + /Sot «(% — z 3 ) 

where A x + A 2 = 1, 

which shows that A x and A 2 are positive, less than 1, 
and that the sign of the difference A B — A x is the 
opposite of that of the expression h + «w -|- am — pm. 
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First Case: z 9 > z 2 . In this case the expression F(z) 
is always positive, z is a diminishing function of the 
time and formulae (16), (18), and (19) become 



°^ z-z x f \z-z. 






■x z . 
du 



«/»(«— 2J [u — z^j 



.h °\« — z,J \u — z, 



n\ 



(20) 



It is to be seen that the space occupied increases 
more and more, the population reaches a maximum and 
then diminishes. Two results are possible. It may be 
that p disappears for a value of z for which s is less 
than 1, The third of formulae (20) shows that this 
occurs in a finite time t v What happens afterwards? 
When the value s = 1 is reached s no longer varies 
and equations (10) and (11) may be replaced by 

s = z, [t > t^ 

Equation (9) becomes 

p> = - m p - hp l 
which gives 

mp x e~ •'<'- '<> 



m -f hp x — hp x e~ m ~ (l) 
When t -^ 30, p tends towards and hence 
lim p — o, lira s = x 



(21) 
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Second Case; z t < z < 2 a . The expression ~F(z) is 
always negative, z is an increasing function of time 
and integration gives 



s = s. 



%b %\ \ ' I %v %i 



z»—z 



n z n — z 



A-i / p. *■ \ A a 



p = - - S Q Z 

h \z— z x J \ z* — 2 

1 f du 



<*& J (^2— «) («— %)^M 



(22) 



In this case the density of the population falls; the 
space occupied may reach a minimum but ultimately 
increases. The final result is always expressed by the 
formula 

lim p = 0, lim s < i (23) 

t~>y> <-»■:» 

Consequently the process always ends in the dis- 
appearance of the population. 



6. Population and Area: Second Case 

Another case, very different from the first may now 
be considered. Suppose that the dead individuals are 
definitely disintegrated by the end of a certain time 
and that the space they occupied becomes available 
once more. To fix our ideas put : 

= 06*-* 

where t — r is the time that has passed since the death 
of the individual. The first three hypotheses of the 
last section remain unaltered. In these conditions the 
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variation of p is again expressed by equation (9) ; 
equation (ro) may be replaced by 

s = a.p + m<x,\ p(u)s{u)e° {H ~ l) du, (24) 

J 

and equation (11) by 

s' — — as -\- a.p{a -\-n-~ns-~ hp) (25) 

The inequalities 

1 > s > wp 

are always satisfied. The difference between this case 
and the last consists in the existence of stable stationary 
state p , s , satisfying the equations 

n-{m + n)s —hp = o) .„ 

— as + «p[o + n — ns — hp) — o J ^ ' 

In the former case the population p finally disap- 
peared at the end of a finite time ; here an indefinite 
existence is possible. 

These results hold good when /S is expressed as other 
functions of the time. If p(t) has a lower limit 

fi(t) > B > 

the population p always disappears at the limit. If 
on the contrary p(t) tends towards zero, a limiting 
state that is stationary and stable becomes possible. 
These conclusions have been reached on condition that 
the metabolic products have no toxic action. If such 
action appears the disappearance of the population 
always occurs at the end of a longer or shorter time. 



CHAPTER VII 

FURTHER MATHEMATICAL TREATMENT 
OF POPULATION-GROWTH 

i. Population Composed of Different 
Age Groups 

Vital co-efficients always depend on age. Infant 
mortality in man considerably exceeds the average 
mortality that appears in biological equations. More- 
over, in civilized society limiting factors act indirectly 
on children, through their parents. 

"With animals that exhibit metamorphosis, as for 
example, with insects, it is always necessary to con- 
sider the larvae as well as the imagos, and the larvae 
often live in a different land of environment. It is 
therefore quite impossible to make use of the logistic 
equation to represent all these cases. Each separate 
population must be treated individually with all the 
influencing factors and all its own peculiarities. A few 
special cases will be considered. 

Insects with Complete Metamorphosis. Imagine as 
an example a population of insects with complete 
metamorphosis, the phases of which cannot co-exist. 
Let l(t) be the number of larvae, n(t) the number of 
nymphs, i(t) the number of imagos, and o(t) the number 
of eggs at a time t. Subdivide the complete annual 
cycle into four periods : 

larvae -> nymphs — > imagos -> eggs 

Begin with the larvae. Their original number l(t ) 

93 
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is proportional to the number of eggs laid by the 
previous generation. The co-efficient of multiplication 
of the larvae is simply equal to their mortality; the 
limiting co-efficient may exist in several cases and one 
may put 

which gives the modified logistic law 

m = *(4M . 

at the instant of the change to the next phase the 
number of nymphs is 

the co-efficient k x < i expresses the proportion of 
full-grown larvae transformed into nymphs. 

The nymph phase lasts for a time (t v t 2 ) during which 
the number of nymphs diminishes according to an 
exponential law, 

At the end of this phase the number of imagos is equal 
to 

*■&) = -My « «>-»*-<>'; (2) 

the co-efficient A 2 < i expresses the proportion of full- 
grown nymphs transformed into imagos. 

During the imago phase things go on as in the larval 
phase, that is to say the decrease of population is 
governed by the modified logistic law 

m « i&K (3 ) 

where m 3 and h s are the co-efficient of mortality 
and the limiting co-efficient belonging to this phase. 
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Suppose that eggs are laid at about the moment t 3 . 
The population at this time is equal to i{t z ) and the 
number of eggs laid is proportional to i(t 3 ), 

o{h) = Vfe) (4) 

The number of eggs diminishes during the tune (t z , tf 4 ) 
that elapses between laying and hatching according 
to the exponential law 

o[t) = o{t^)e~ m ' {i -^ 
which gives at a time t t a number of larvae equal to 

« = vw = v(gr»*- f j (5) 

where A 4 'is the- proportion of eggs that hatch. 

Replacing o{t s ) by (4), ifa) by (3), etc., we get the 
relation between the number of larvae at the beginning 
of the next cycle and the original number l(t ). 

__ kkkk m i m sKk) e ~ " hik ~ h) ~ " h{k ~ h) ~ nh[h ~ ' J ~ " hlh ~ h) 

H- htfnJit^e- mAk ~ '•> - »•('•-*>>. (1 — er '»><'>- '■>) (6) 

At the end of a large number of cycles, Z tends 
towards a limit L 

_ m^m^X^X^e ~ '"^ ~ '•> ~ '"* ~ '■> - '"><'• ~ '»» " m < ( '< - '■> — 1 

as long as 

/I /I Jl /L > e'"'!' 1 " '•> + ">>('» -'1) + wil*.-M + »M'i- '.) (8) 

This condition is equivalent to the inequality 

s = n—m >o 

which is an essential for the existence of a stationary 
state other than zero in the case of a logistic population. 
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The same result will be found for the numbers of the 
imagos at times 

<g, k + T, U u + aT, . • • 

differing by T, the period of the whole cycle. 

If the original value l(t ) is equal to L, there is a 

stationary state in the sense that at the end of a cycle 

l(t + T) = l{t ) = L. 

If l(t Q ) is less than L, l(t Q + wT) increases and tends 
towards L, if l[Q is greater than L, l(t + wT) decreases 
and tends towards L, It may again be remarked that 
the stationary state here is the periodic one with a 
period T so that 

l(t + T) = l(f); i(t + T) = i(f) 

o(t + T) = o(t) ; »(* + T) =«(*). 

The most complex cases, where a certain amount 
of superposition is possible, may be treated in the 
same way. 

Case of Co-existent Phases. Suppose that the cycle 
is not a long one and that the co-existence of phases is 
possible. In this case the number of nymphs becoming 
imagos is proportional to the number of larvae, and 
the number of eggs hatching is proportional to the 
number of imagos. In this way the equations of the 
system may be written 

v { t) = n$ - nM - m) - VW \ fe v 

i'{t) = - wit) + k,l(t) - h^(t) ] w 

assuming, naturally, that the larvae and the nymphs 
belong to different stations and consequently that the 
limiting co-efficients h^ and h n are zero. 
Two cases may be distinguished. In the first the 
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iertility v and the co-efficients of metamorphosis 
k± > k B are large enough to diminish the effect of the 
two mortalities, so that 

equilibrium between the two phases is possible and at 
the limit 

Hm I = L, lim i = I 

where L and I are the positive solutions of the algebraic 
equations 

vl = L(A X + Mi + &u L ) ; A 2 L = 1(^2 + ^a 1 )' 
Passage to the limit occurs without fluctuations. 
In the case where 

the two positive factors mentioned are not sufficient 
and the population vanishes in the limit. 

• 2. Animais with Seasonal Reproduction 

From many examples dealing with animals with 
seasonal reproduction we shall choose S. Sewertzoff s 
observations on grouse and black game in the National 
Park of Bashkiria. These populations are composed 
of distinctly different age groups, especially after 
hatching. Sewertzoff devoted attention to the young 
group in the period following hatching (the Russian 
summer months of June, July, and August). He 
showed that the average number of young per nest 
diminishes according to the exponential law. This 
decrease is due to many causes, but chiefly to the 
activity of many predators. Our equations give a 
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satisfactory explanation of this result. Let x be the 
number of young, y the number of nests, proportional 
to the number of adults, and z the number of predators. 
It must be assumed that at the period under considera- 
tion there are no unhatched eggs, and hence the popu- 
lation x does not increase; it decreases from very 
diverse causes such as the death of the young and the 
complete disappearance of nests owing either to the 
death of the parents, the severity of the weather, or 
to predators. Each of these factors has its own place 
in the equations. As far as predators are concerned, 
they are in no way restricted by the presence of the 
adult grouse. In these conditions the equations may 
be written — 

%' = — (fi + m)x — [Xi + X 2 )xz — hxy (10) 

y' = — my — hy % — l^yz (n) 

* / = ~fe_fe a (12) 

These equations apply only to the peculiar conditions 

of the summer months. Let u — - be the average 

y 

number of young per nest. Equations (10) and (n) 
give 

u' ~ — u(fi -\- XJ). (13) 

Assuming that the decrease of the predators is much 
slower than that of the young birds, we get an 
approximate solution of (13) 

u = u{o)e~^ + i ^ (14) 

which reaffirms the statement that the average number 
of young per nest decreases according to the exponential 
law, 
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Using this formula a very close agreement is found 
with Sewertzoff s figures. For 



Date 


June 
7 21 


July 
1 13 23 


Aug. 
2 18 


Average (-observed 
per nest (.calculated 


57 6 
6-2 5-4 


5 4 3-9 
4-9 4'3 3-9 


3'9 2-9 
3-6 3-Q 


and for 


Date 


19 June 


2 Auo. 


17 Aug. 


Average /observed 
per nest (.calculated 


5-4 
5S 


4.6 
4-3 


3-3 
3-9 



These examples show that for animals with repro- 
ductive seasons the dates on which a count is made 
must correspond to the natural cycle ; one can usefully 
compare only homogeneous statistical facts, that is to 
say, those relating to dates separated by a complete 
number of cycles. These facts must be completed by 
a study of the movements of the population during 
a cycle. 



3. A Generalization of the Logistic Law 

A generalization of the logistic law has been men- 
tioned above 

p> = F(# = m - f{p) 

where the function <f>{p) represents ah the favourable 
factors and y>(£) all the unfavourable ones affecting 
the growth of a population. We shall consider a par- 
ticular case which to some extent includes the attempts 
to improve the logistic law. The law is often criticized 
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on the grounds that the decrease of fertility and increase 
of mortality are often more rapid than the logistic law 
implies. Put therefore 

<f> = npe~' p yi = mpe^ (16) 

which gives 

p' =^p[ne- rp - me" p ] (17) 

Two very different cases exist. In the first n <m and 
whatever the original value of p , p' is always negative 
and p tends towards zero. 

In the second case, where n > m, the equation 

ne~' p — me" p = 

has a root co. 

t n 

log- (18) 



v + n m 
Suppose p < co ; the solution of equation (17) 



t-t Q = 



rP ds 



„ sine "~me^ 



(19) 



gives the time t expressed as a simple integral. The 
population p tends towards a limit co. 
Its asymptotic expression when t is large enough 

p ~ co — Ae~ H (20) 

with rj ~ me 1 "* {v + /n) 

is analogous to the logistic law. Moreover, when p 
is small, equation (17) may be replaced at the beginning 
of the process by the Malthusian equation. 

p' = (» — m)p, p = p e("- m)t (21) 

The essential peculiarities of the logistic law thus 
appear again, as we have already shown, save for ease 
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of calculation. This latter argument might be disre- 
garded if it could be demonstrated that equation (17) 
represents the true state of affairs better than the 
logistic law. We have shown that the logistic law is 
well supported by the facts. Moreover, with formula 
(17) the difficulties of calculation become insuperable 
for mixed populations, and in reality homogeneous 
populations do not exist. Hence it is better to make use 
of the logistic law. 

4. Senescence of a Population 

It is often stated that a population, like an organism, 
is subject to old age expressed by a slowing down of all 
vital functions.- Introduce into equation (15) a dimin- 
ishing factor a{t) representing senescence 

f = F(p)o(t) (22) 

This equation is easily reduced to (15) by a change 
in the independent variable. Put 

-( 
u = a{t)it (23) 

J '0 
We get 

dp __dp du _ dp 
dt du dt du 
and equation (22) becomes 

Another form of the equation of senescence may be 

p' = F(p) + F 1 (P)e-" t (24) 

Admit that 

Fo(#) > °- 
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as long as p is less than co, this value being the smallest 
root of the equation 

Fo(0 = o. 

Admit further that F' [co) < o, F^w) > o and that the 
functions F and F x may he expanded by Taylor's 
theorem in the neighbourhood of the point co. As 
t -> x , the population p tends towards co. Two cases 
may appear : 

(i) A = -F' H<(t; 

in which case we have the following asymptotic form 





p ~ m — ce~ u . 


(ii) 


1 = -F' a (a>) xr; 


in which case 






r a—X 



This case does not therefore differ qualitatively from 
the logistic law. A population obeying law (24) always 
has a maximum w which it reaches at the end of a long 
enough time. 

Finally a special case of equation (22) may be 
considered : 

P' ~epe- at . ■ (25) 

The solution is 



;('■ 



-<rti 



P=Pf" (26) 

which, for t-> x, tends towards the limit 

e 

limp = p Q e" 
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Consequently , introduction of the factor of senescence 
in this case gives the same result as the limiting term in 
the logistic equation, but with a very important difference: 
the final value of p is here exactly proportional to its 
original value, while in the logistic case it is quite inde- 
pendent of it. 

This independence appears to be well established by 
observation. Hence, almost without exception, the 
introduction of the factor of senescence into the theory 
of the growth of a population is not justified. On the 
other hand this factor may be of great value in the 
theory of the growth of an organism. 

5. Dynamic Theory or Population-growth 
The chief feature of all these theories of population- 
growth is their kinetic nature. We speak everywhere 
of actions and of effects but nowhere do we introduce 
the idea of force which is so essential both in physics 
and astronomy. The analogy has suggested to several 
writers the idea of constructing a differential equation 
of population-growth as one writes equations for the 
movement of a body or a medium. In such equations 
all the forces acting on the system are balanced up: 
(1) external forces ; (2) internal forces and resistances ; 
(3) inertia which balances all the preceding forces. 
From the mathematical point of view, this 'balance- 
sheet ' is expressed, in the case of a single variable, by 
a differential equation of the second order. 

3-* + '-* 

E denotes the whole of the external forces, I the 
internal forces, and R the resistances. 
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In Quetelet's Physique sociale there is a very inter- 
esting discussion of this problem: "The theory of 
population may be reduced to the two following prin- 
ciples, which I regard henceforward as the fundamental 
principles in the analysis of population-growth and the 
causes that influence it. Population tends to increase 
according to geometrical progression. The resistance, 
or the sum of the obstacles to its growth, varies, other 
things being equal, as the square of the rapidity with 
which the population increases. Thus, when a popula- 
tion can increase freely, without hindrance, it grows 
according to a geometrical progression ; if the growth 
takes place among obstacles of all kinds, which tend 
to check it and which act in a uniform manner, that 
is to say, if the state of society is not changed, the 
population does not grow indefinitely, but tends more 
and more to become stationary." 

Inspired by these ideas of Quetelet, Delevski writes 
■the equation of inertia in the- absence of all forces as 
follows, 

P" = ep', (2 9 ) 

which evidently corresponds to the Malthusian law 
p'=*P<fi*. 

Adding to this terms expressing different positive and 
negative actions, he writes the differential equation 

p"= aa —bp + (e — a')p' - a"p'K (30) 

In this equation a—bp represents an external force, 
and the terms a'p', a"p' 2 correspond to friction or to 
resistance of the medium in equations of mechanical 
movement. 

I do not wish to criticize or condemn these attempts. 
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One knows the opposition encountered by Newton's 
ideas in the scientific world of the seventeenth and 
eighteenth centuries. As far as the idea of force is 
concerned, we return to pre-Newtonian times, and in 
my opinion it is a wise course to take, but one must 
not forget all the good this idea has done for two 
centuries. However, it seems to me that equation (29) 
may be given a preferable form, 

P" = fy. (3i) 

that is to say that the inertia of population-growth is 
proportional to p. Further, the general equation of 
type (30) should be verified by the logistic law but is 
not, and above all there is lacking a resistance simul- 
taneously proportional to the population and to the 
rapidity of its growth. With these corrections the 
equation may be written 

P" = «oo + ( e2 - «io)£ ~ %#' ~ a nPP' ~ <W>"- (32) 

This form is more complete than equation (30), 

more logical, and, in addition, it is verified by the 

logistic law for the following special values of the 

co-efficients 

Another advantage of this form is that it can easily 
be integrated when a m = 0, a 10 — e 2 . 
In this case equation {32) becomes 

P" + Oaf + <*nPP' + <W>' 2 = 
or 

(P" + <W>' V"* + HxP'e aoP + ari>'P'"" = ° 
or again, after integration, 

p> = He-«*-^-^j> + ?± (33) 



io6 

and hence 
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t = 



rp 



ds 






(34) 



*02 



*02 



The function j!> given by this equation is qualitatively 
analogous to the logistic law but is more supple. 

Equation (32) may be completed by integral terms 
expressing the accumulation of residual resistances. 



CHAPTER VIII 

RELATIONS BETWEEN SPECIES 

i. General Considerations 

the problem of the growth of a mixed population 
recalls that of the development of a homogeneous one ; 
it differs, however, to a considerable extent because 
the relations between different species naturally take 
a more complex form than in the case of a homogeneous 
population. This disconcerting complexity seems to 
have evaded all analysis and in particular mathematical 
analysis. However, in the majority of cases these 
relations may be represented by differential or integro- 
differential equations as simple as those studied in the 
preceding chapters. 

It is not a matter of relations between species on 
the whole of the earth's surface, but between groups 
in a station that is more or less isolated, or capable of 
being isolated by a familiar mental process. In this 
complicated microcosm there exist supplies of food, 
shelters of very many kinds, and climatic conditions 
subject to seasonal variation. This assemblage can 
often be separated into microcosms of a lower order, 
with almost no visible connection between them. 

The species that occur in a station may be found 
in a large number of others in similar but not identical 
circumstances. These changes in conditions are ex- 
pressed in variation of the vital co-efficients. But the 
future of a mixed population depends on these co- 

107 
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efficients, so that the probability of this or that result 
of the struggle for existence may be estimated. Thus 
we can to a certain extent foresee a priori the future of 
one or of several species in a set of stations of a certain 
type. 

In addition to the permanent factors we shall con- 
sider intermittent factors, the effects of which appear 
when one or more groups exceed a maximum or fall 
below a certain level, as we saw in the case of small 
rodents subject to epidemics. 

In all these problems two opposed tendencies may 
appear : 

i. The group of organisms occupying a station may 
tend towards a state of equilibrium. 

2. Equilibrium may be impossible. 

In the first case a final stable state is reached in 
practice at the end of a longer or shorter time and is 
maintained thereafter. This state is independent of 
the initial state when equations 25 of Chapter II permit 
only a single set of solutions in the positive region of 
biological space; in other cases the system tends 
towards one of the states of stable equilibrium that 
depend indirectly on the initial state : finally there are 
certain cases where the final state depends directly 
on the initial state. 

The stability of a biological state is not permanent. 
The vital co-efficients vary because of imperfect isola- 
tion, the appearance of mutations, the intervention 
of a new factor, and other causes. In human society 
technical development, improvements in the circum- 
stances of life, a sudden or a continuous variation of 
the social organization, all have a repercussion on the 
vital co-efficients. 
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In the case of instability of a biological system 
separate mention should be made of the more or less 
periodic oscillations about a mean state that is neither 
reached nor approached. It may happen that these 
oscillations, like Lissajous' figures, may entirely fill a 
region or an enclosed area in biological space. This 
dynamic equilibrium has no secular stability. 

It is to be observed that the relations between species 
are not limited to competition and that besides the 
struggle for existence there is co-operation in many 
different ways — unilateral (parasitism), or bilateral 
(symbiosis), sometimes loose or fortuitous, sometimes 
inseparable. This phenomenon is of great importance, 
for upon it depend problems of relationship among 
multicellular organisms, 

The vital co-efficients are extremely unstable and 
the following example gives some idea of their varia- 
tion. According to Kalabuchow, small rodents like 
mice, which in nature follow a plainly seasonal cycle, 
are quite capable in better circumstances of reproduc- 
tion at all times of the year. A little calculation enables 
us to evaluate the variation in the co-efficient of multi- 
plication e. In natural circumstances this co-efficient 
is of the order of one per annum. Seasonal reproduction 
gives an annual birth-rate of about twenty per annum, 
while optimum reproduction gives n ~ 200 per annum. 
This means that in optimum conditions s is of the 
order of 180 per annum. Consequently it is quite 
impossible to look upon the vital co-efficients as an 
unalterable characteristic of the species whatever may 
be the circumstances of its life. 

To understand the construction of the vital co-effi- 
cients, consider a relatively simple case. Suppose that 
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it is concerned with two species in perpetual struggle 
for food, and each eating the other on chance encounters. 
Let p t and p s be their respective numbers. It is always 
possible to make a numerical estimate of the probability 
of the victory of one or the other adversary, and also 
that of an undecided or drawn battle. For example, 
let 77j_ be the probability of the victory of one (fij), 
7T 2 the probability of the victory of the other, (ft 2 ) 
and 77- 12 the probability of an indefinite result. Then 
evidently 

^i + ^2 H - Tus — I - 

Let hpip 2 be the number of encounters in unit time. 
The direct damage from the struggle may be expressed 
by the term — h^Pifiz introduced into the life-equation 
of the group j> x and by an analogous term — htr^p^ 
into the equation of group p 2 . This is not all. En- 
counters that have no result have an effect on the 
fertility and mortality of the two groups. Further, the 
fact that an adversary has been eaten is not without 
influence on the vital co-efficients. This influence is 
proportional to the number of eaters and the numbers 
of the eaten, and also depends on the physiological 
utilization of the food. If the organisms in the two 
groups are physiologically similar, the terms + ^ipip% 
and -f ktr^fz may respectively be added to the life- 
equations of the two groups. Finally, the result of 
neutral encounters is expressed by the terms — hr ls p J j> i 
and — iTTtipip?.. 

In consequence the differential equations of the two 
groups assume the form 

Pi = hti - IhiPx - {hir z - hn x + hrubprfa ) , . 
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or Pi = £ih — h ilPl 2 ~ h l2pifa ] ( 2 \ 

P'z = %Pz — hiPiPz — KiPi J 

where A 12 = Jvrr 2 — ktr x •+• Itt X2 
and /% = hrr x — kn 2 -j- ^tt 12 . 

Further the quadratic terms h^p-? and /; 2 2^>2 2 n °t onr y 
express the reciprocal limiting action of the individuals 
of the same species, but all the external and internal 
limiting factors. 

Thus this particular and relatively simple case of 
two predatory species that are physiologically similar 
gives co-efficients of inter-action h 12 and h 2X which are 
not identical, and are positive or negative according to 
circumstances. In more complex cases complementary 
terms corresponding to every activity that occurs in 
a microcosm find places in the biological equations. 
According to what has been seen in the preceding 
chapters all these terms are of the order o, i, or 2, in 
relation to the numerical effectives of the groups. 
Terms of a higher order can, however, be foreseen : we 
may complete the system of two predatory species by 
a third group p s which feeds solely on the remains of 
the combats between j> x and p 2 . Then the advantage 
that it receives from this is proportional to the product 
of the three effectives p x ps,p 3 . 



2. Whelks, Mussels, Barnacles 

We shall make a study of an extremely interesting 
cycle observed by E. Fischer-Piette. It concerns the 
relations between shell-fish (whelks and mussels) and 
sessile crustaceans (barnacles) observed during a period 
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of ten years on the mussel-beds of the Pointe du 
Decolle\ The study of these relations enabled Fischer- 
Piette to perceive the interactions of the three groups 
after a disturbance of the faunal equilibrium, which 
had the result of progressively restoring the first 
equilibrium state. Moreover, the phenomenon of 
periodicity in mussel-beds is well known. Thus, Fauvel 
found that in the neighbourhood of St Vaast-la-Hougue 
the mussels from the Bay of Isigny from time to time 
rapidly colonize the neighbourhood of St Vaast and 
then disappear before the attacks of their enemies, 
man and particularly starfish. This process is made up 
of the following phases- — 

1. Rapid development of the mussels. 

2. Starfish arrive from all round and, finding an 
abundance of food, multiply with astonishing speed. 

3. Their voracity and the destructive action of man 
quickly obtain the upper hand of the mussels. 

4. As prey fails the starfish tend in their turn to 
disappear, leaving the held clear for a fresh invasion 
of mussels. 

5. A new cycle begins. 

In the case observed by Fischer-Piette the phenomena 
are more complicated. Whelks can feed on mussels or 
on barnacles, but have a marked preference for bar- 
nacles and eat mussels only in the absence of other food. 
This predilection is the more surprising because a diet 
of mussels seems to be much better for them. We shall 
show that this preference can create and maintain 
periodic oscillations ; but first we shall rapidly set out 
Fischer-Piette's results. 

1. The cycle begins by the reproduction of the 
barnacles, which cover the rocks with a dense popula- 
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tion. This provides food for whelks, which ignore the 
barnacles below a certain size. The whelks eat barnacles 
fixed to mussels, not the mussels themselves. 

2. The arrival of an abundant hatching of mussels 
causes the death of the barnacles over a wide area 
where the latter are swamped by the mussels. 

3. The whelks which had previously eaten these 
barnacles are now deprived of them. They change their 
diet and begin to pierce the mussels. They multiply 
enormously ; although they only adopt their new diet 
under pressure of absolute necessity it is much better 
for them. This increase of whelks results in the rapid 
diminishing of the mussel bed. 

4. When the time comes for the barnacle larvae to 
settle, surfaces freed from mussels by the activities of 
the whelks are covered afresh with barnacles. 

5. As soon as these barnacles have grown big enough, 
the whelks return to their original diet and start to 
eat them. The mussels, left alone by the whelks, are 
henceforward safe from destruction. The cycle begins 
again. 

Thus, if the mussels happen to increase abnormally, 
and so have disturbed the equilibrium, one direct 
consequence of this disturbance (the suppression of the 
barnacles) is to let loose an antagonistic factor — the 
whelks begin to destroy them. This destruction will 
go as far as complete extermination if it does not bring 
into play another factor to curb the action of the first : 
the return of the barnacles which attract the whelks 
from then- prey. 

Let us now try to express these relations in differ- 
ential equations. Let ft x be the number of barnacles, 
p s the number of mussels, and p a the number of whelks. 
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Wiist Stags. The first stage may be summarized by 
the following differential equations : 

P'i = hPi ~ KPiPz — hsPiPz - hiPi 2 \ 

P'l = %p2 + KPlPs - hiPi (3) 

P's — B 3ps + htPiPz — hsp3* ) 

These equations refer to the initial period, character- 
ized by an abundance of barnacles. The absence of 
the variable p 2 from the third equation implies that 
the whelks prefer to feed on barnacles (the term h 3 p x p^. 
The term + KiPiPz i n the second of equations (3) shows 
that the mussels' smothering of the barnacles has begun. 
Equilibrium is impossible in this phase if the fertility 
of the barnacles is not great. It is very probable that 
the inequality 

e^^f+e^f (4) 

"22 "33 

is satisfied at this stage. In consequence, population p x 
diminishes, and when it falls below the level p x < tt x 
relations are radically altered. 

Second Stage. The barnacle population now exists 
only sporadically, and the whelks begin to eat the 
mussels. Then : 

P'2 :=z: %Psi~'K! i3 p z p s ~B-Mpr? (5) 

P's = E3P3 + H 32 M)- H 33 jV J 
During this time the multiplication of the whelks is 

extremely rapid and it is probable that E 3 satisfies the 

inequality : 

E 3 > e 2 g?. 
■n-23 

This implies that equilibrium is impossible during 
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this stage also and that population ft 2 diminishes, 
liberating rock surfaces for the return of the barnacles, 
which little by little re-establish themselves. At this 
point the third stage begins. 

Third Stage. The mussels almost disappear and the 
differential equations assume the form 

P'i = hPi— VuPifa- VuPi 2 ] 

p' 2 = E 2 p 2 — r] Z2 p 2 2 (7) 



P'z = £ 3^3 + VsxPiPs — %sP; 



2 



3 



The barnacles are eaten by the whelks, but the 
second attack of the mussels is only to be seen towards 
the end of this stage. Thus the relative equilibrium 
establishing itself during this phase is upset, and the 
cycle begins again. 

Fischer-Piette's memoir gives no figures except for 
the duration of each phase. From 1925 to 1929 mussels 
were almost absent : thus the third stage covered four 
years. In 1929 equilibrium was upset by the arrival 
of a large hatching of mussels. In 1932 the whelks 
began to perforate mussels, so that the first stage 
occupied three years. The shrinking of the mussel-bed 
has been observed since then. 



3. Residual Action 

In Chapter V we gave a complete numerical treat- 
ment of a special case of residual action; it was con- 
cerned with the poisoning of a population by accumula- 
tion of its own metabolic products. Other residual 
actions, accumulating and altering the vital co-efficients, 
can easily be suggested. For example, theoretically 
one can foresee the accumulation of hereditary effects 
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produced by the action of the environment on the 
organisms, and thus express Lamarck's theory in the 
form of equations. 

We are here limited to the study of reciprocal toxic 
action between two groups, p x and p 2 , composing a 
mixed population. Replace the co-efficients of multi- 
plication by the following equations : 



E 1 = e 1 - c n J k n (t — s)p 1 (s)ds — c X2 

kjt - s)p 2 {s)ds 

k zl (t—s)p 1 (s)ds — 

t 
K4$ — s)p 2 (s)ds 



Ej = £ 2 — <, 21 I rt 2lV (. — o )Fl ^)uo — i- 22 



J 



(8) 





The differential equations may be changed into 
integro-differentials 

P'x = A( E i - KPx ~ KPi) 
p\ = p 2 (E z — h^pi — h 22 p 2 ). 

Common sense insists that a residual action should 
diminish with time, and the functions of residual action 
k ih {ii) are decreasing functions. Further, the presence 
of the co-efficients of residual action c u , c 12 , c 2l , c 22 , 
allows us to write, in general, 

k ih {o) = I. 

The form of these functions can be theoretically 
deduced from experiments ; actually it is not possible 
to quote any example in which this research leads to 
a satisfactory result. These functions can then be 
given the simplest mathematical form. It has been 
seen that the solution of these integro-differential 
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equations is a very difficult matter, even in the case 
of a single variable. Accordingly we here confine 
ourselves to the study of a few simple examples in 
which two first approximations are sufficient and these 
approximations can he satisfactorily calculated. Put 
throughout 

k n (u) = kjaHu) = k zl (u) = k 22 (u) = 1. (10) 

Example 1. Imagine a case of two species with zero 
co-efficients of internal competition, that is to say a 
mixed population of low density. Equations (9) 
become : 

P'^P^-KPi) i (II) 

Expand j> x and fi 2 by Taylor's theorem in respect of 
the parameters c n , c 12 , c 21 , c 22 , . . . 

Pi — P10 + c nfiu + c nPiz + CmPiz + c wPu 4- • 



(12) 

Pi = P20 + C llAl + C 12p22 + CoPm + C 22pU + • • • / 

Substituting these expressions in equations (11) gives 
as a first approximation 

P'io = Pio{ei — fhsPso) ] ( „n 

P'M = p2o(e2—KPio) J V ; 

which only differs in notation from equation (18), 

Chapter II. We then get the following equations for 

the approximations of the higher orders, 

^i = qj p^- h ^Pio^-Pi<Ai 

/; > (i4) 

P20 

where q x and q z denote a couple of corresponding 
functions in expansions (12), and A ± and A 2 are func- 
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tions depending solely on approximations of a lower 
order. These equations can be easily integrated. 

Multiply the first by -^-, the second by -~^-> and 

AoAo A0A0 

take the difference of the two results, we get 



^-^■-»b(s) + *s(s;) 



Ai^'ao Ag^'j 



10 



Integrating this we have 



y 20 



P: 



'10 



qd>'io— ViP'zo = PiaP: 



20 



Ai^'20 A 2 ^' 10 " 



ds = B(t) (15) 



L Ao Ao J 

in which B(tf) denotes the second member. In the first 
of equations (14) put the function q 2 determined by 
equation (15) , replacing p' w by its value from equation 
(13) and we get 

V'lP'w— 1iP"w Ao A i W>io- 



A 2 



■ 10 AlO 

putting, for simplification, 



\' 2 
' 10 



B 



B(t) = p 10 p w 



'Ajp'n A2P 



^20 



Ao J 



ds. 



Integrating this equation we get 

J L P »(*) J 

and an analogous expression for q % 



&o(*) 
Aio(*) 



dz (16) 



= -A: 



20 



"a,W-^ bw 



Aofc) 



dz (17) 
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These formulae allow us to calculate the second 
approximations. Put 



Pio = 



pjo(s)ds 



p 20 (s)ds 



Bu(«) = Ad(*)Ak>(*) 
Bsi(x) = Ao(#»(*) 

B 2l{ Z ) = Plo(z)p2 (z) 



A 



10 



"lO* 20 



y 20 



* P_20^io 
Ao 

*20.P 21 



<fe, 



ds. 



<fc, 



<fe. 



.0 Ao 
We then get the following expressions 

Plo( Z ) 13 /„\J„ x Al^'lO B 2l(0 



^"^'-lo^-W 



BgiW'te Ai = - 



# 



20 



P\ 



20 



(18) 



p 22 — h n p 20 
A3 = ^Wio 

i 

Pu = kiaP'w 



&o(*) p ,., j, _ JW* M 

' Ao(«) 



A 



B #(x 9 ) 



20 



£ 



20 



*v(*) 

Ao(*) 



Bu(«)& ^23 = 



A^'20 B u (0 



^'10 
iw'20 



A: 



10 



(20) 



7— ¥w 

P 10 f 10 



iV(*)' 

The remaining approximations can be calculated in 
the same way. Formulae (18) to (21) are not very 
convenient, but they allow several qualitative conclu- 
sions to be drawn in any case of a concrete biological 
problem in which auto-concurrence can be neglected. 
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Example 2. Consider another case in which the 
later approximations follow from the possibility of 
integration in the first approximation. Let two species 
be fairly close, so that 

h n h 22 —h n h 21 =o. ( 2Z ) 

Admitting the possibility of expansion (12) we get 
in the first approximation the equations — 

P'w = Pioih— hiPio~ hzPw) 
P'w = Aota ~ K\Pm — htPm) 

We easily get the relation 

P10 —Pw _<«.*., -«.*..H 



(23) 



«■/'" a 2 h " 



(24) 



where a x and a 2 denote the initial values of p x and p 2 . 

Let jj and q 2 be any pair of corresponding functions 
in expansion (12) ; we get 



it 
l'i = "V^ - IhxPxoli — hiPMz — AA 

?' 2 = ~br^ ~ Vrfi — KPrnQi ~ Ao A 2 



(25) 



in which the functions A x and A 2 only depend on the 
preceding approximations. From these equations we 
get 

^s(£)- *»£(£) -M.-V*. 

which gives after integration 

A^s (26) 




Ao Ao 



AnflS ~~ ftgil 
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This equation, together with (25), enables us to obtain 
q x and § 2 in the following form : 



ft = An(0 [fto(*) c (*)*2a - A 1 (^)]F(i, «)<?* 
J 

fc = - A»W J o [Ao(*)C(^ 21 + A 2 (*)]F(*, *)<fe 



(27) 



in which 






Aorfs — 



A t ds 



Apply these formulae to the search for the second 
approximation. We get 



Ai = KPw ?(*, *)&»(«)& *V S 

Jo Jo 

P10 W21 ,, 



Ai = 



^ai L P\ 



20 



P in ds 



(28) 



J 
Ao 

21 



JTonrfS 



•20' 



P. 



P BO tfe 



^21 L i>20 

Aa = M10 F(*. *)Pm{*)&* *V S 
Jo Jo 



Pw = 



' KP. 
L Ac 



*n 



(29) 



(30) 
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Pu = hlPlO F(*. #2o(*)^ P 20^ 

Jo Jo 



•4, -tn 

Fa — 7, 



'K1P1, 



V L £ 



A, 



■10 



KJS 



} 



(31) 



These formulae show that when /%, h i2 , h ilt h 2Z , 
e ± and £ 2 are positive, the functions p n , p lz , p 2 a, Pa are 
always negative and the functions fi ilt p 2i , p^, pn 
always positive. From this we immediately get the 
foEowing inequalities. 

\ 



Ai< 



A 3 < 



KP, 



"ii 



KP11 
K 



20 ! Pj^S 



Pm^S 



Aa < 

Pli< 



hiP: 



ao 



hiPw 



Pods 



h 



21 



(3a) 



P w ds. 



J 



These results can be expressed in ordinary language. 
In the expansion (12) the term p 10 gives the size of the 
population p x in the > absence of toxic actions ; the 
terms c u p n and c 12 p 12 express respectively auto-intoxica- 
tion and intoxication of group p x by the metabolic 
products of group p 2 ; hence it is natural that these 
terms are negative. On the other hand, the terms 
c ziPis an d <W>u express the advantage which the group 
derives from the auto-intoxication and intoxication of 
its neighbours; it is thus natural that these terms 
should be positive. Similar reasoning applies to group 
p 2 . The inequalities (32) give a rough expression of 
these advantages and disadvantages. A mixed popu- 
lation of n groups can be treated in the same way. 

Mixed Bacterial Populations. As material for ex- 
periment there are the figures published by J. Regnier, 
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and the additional information that he has been good 
enough to send me, concerning a population composed 
of Bacillus coli and Staphylococcus. Regnier's data 
concerning pure cultures have already been used, and 
we supplement these results by what can be learnt 
from mixed cultures. From a qualitative point of view 
Regnier's figures confirm the theory. 

Let p x be the number of B. coli and p z the number of 
Staphylococcus in the mixed culture. Using the values 
for the vital co-efficients 

% = 1-14 h lx = -0003 c u = -023 

e 2 = 1-05 h 12 = -ooo2 c aa = -0265 

calculated from pure cultures we can calculate the co- 
efficients of inter-action, K& ~ -00005 and h zl ~ -00014, 
but the distribution about the average is very wide. 
These values show that in all cases Staphylococcus 
suffers more from living in company than does B. coli. 
The much more rapid decrease in the population p 2 
shows that the toxic action of B. coli on Staphylococcus 
is greater than the reverse action, but for the moment 
it is impossible to calculate the co-efficients of residual 
action with any certainty. 



CHAPTER IX 
RELATIONS BETWEEN SPECIES {continued) 
i. Two Species Competing for the Same Food 
Consider a population composed of two species fi x and 
p z in competition for the same, sufficiently abundant, 
food. Each meeting between the competitors may be 
followed by a more or less definite struggle or may 
produce only a mutual inconvenience. Repeating the 
reasoning of the last chapter, we get the equations : 

P'x = e xPx ~ hxPx % — K-lPxPi \ /js 

P'z = Hpz ~ hiPiPz — h 2ZP* i 

The co-efficients e x and e 2 are positive, for the two 
species can very well do without each other. The 
co-efficients h n , h^, h 21 , and h iZ are positive, for it may 
be supposed that in every case competition between 
individuals is a net loss to the community. 

These equations have been studied by V. Volterra 
and A. J. Lotka. Volterra assumes that the two com- 
peting groups are not widely separated biologically 
and that the limiting co-efficients are related by the 
equation 

MM-*u*n.= o (2) 

This hypothesis implies that, from the point of view 
of competition, the two partial populations j> x and p z 
may be replaced by the imaginary population ri x p x 
+ %Pz- 

124 
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In these circumstances equations (1) become 

P'i = hPx ~ hMViPi + %p2) ] m 

P'a = %h ~ hMVifi + V2P2) j ' 

This hypothesis has been discussed in the last chapter. 
Equations (3) can be integrated : 

^. = ^_ e (.A-M.W (4) 

Pt of* 

where a x and a 2 are the original numbers of the two 
groups. Two possibilities exist. If 

the first species has a marked advantage over the 
second, and at the limit 

T £ 1 I- 

lim p x = - — ■, hm fiz — 

t -> DO 'hVt i -> DO 

In the alternative 

lim^ = o, lim^> 2 = ^- 

Finally, in the most improbable case in which 
£^2 =■ ejt^, the co-existence of the two groups is pos- 
sible. Thus, according to Volterra's hypothesis, equili- 
brium between the two groups is in general impossible, 
and the smallest disadvantage leads to the disappear- 
ance of the corresponding group. 

Lotka studied equations (1) in a general way, without 
imposing condition (2) on the limiting co-efficients. 
Equations (1) cannot be integrated as easily as equa- 
tions (3) but it is possible to point out the general 



136 MATHEMATICAL BIOLOGY 

character of their solutions. The following cases may 
occur. 

In the first case, the vital co-efficients satisfy the 
inequalities 

hht > e As> hhi < e Ai (5) 

which are equally favourable or unfavourable for the 
two species; the two groups then tend towards a 
limiting stable state 

llrv, -A S Jh2— £ 2^12 1,- L £ 2^1 1 ~ £ 1"21 ffS 

f-H>- DO "U«22 — 'h.l n 2\ f ->- 30 /2 U"22 "12''21 

In the second case the inequalities 

£^22 > E z h iz , ejifr > e % h n (7) 

are wholly favourable to the first species, which alone 
survives : 

11:01^ = — lining = (8) 

t -> 30 "U I -*- 30 

In the third case, on the other hand, 

ejAjjjj < ea^a, eAi < e zhi (9) 

the second species is favoured and is the survivor : 

lim& = o, limp 2 = -~. (10) 

< — >- DO ' -> DO "22 

In the fourth case, the inequalities 

£^22 < £A a , 6^21 > £ 2^1 (ll) 

are the reverse of (5). Here things are more compli- 
cated; the process always ends in the disappearance 
of one of the two competing groups, the other tends 
towards a limiting stable state determined either by 
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equations (8) or equations (ro). The inclination of the 
system towards one result or the other depends on the 
original values of a x and a a . It may be said, broadly 
speaking, that a much more numerous group alone 
finally survives. 

These results are very interesting from a biological 
point of view. Equilibrium between two species com- 
peting for the same food is possible, but only in a 
quarter of all the cases which may occur. It follows 
from equation (5) that in this case we have 

"11"22 > "12"21' 

and hence that the struggle between the two groups 
is less violent than their internal competition. The 
contrary would be more probable in the case of two 
allied species. 

We must now apply our remark on vital co-efficients 
and environment. Two species can live together in a 
large number of stations in very varied conditions. 
Vital co-efficients may thus be looked upon as numbers 
chosen by chance, when the probability of equilibrium 
between the two groups is only -25, while that of the 
disappearance of one of them is -75. The presence of 
other groups further decreases the probability of equili- 
brium. Moreover, the vital co-efficients are not alto- 
gether independent, and an increased mortality is 
almost always accompanied by an increase in the 
limiting co-efficients. Consequently the chance of co- 
existence for related groups is very small and in general 
the influence of natural selection is towards an exagger- 
ation of the differences between the- groups. We shall 
see later that this result remains true for groups in 
which living hybrids may occur. 



128 



MATHEMATICAL BIOLOGY 



2. Three Species, Two of which Feed on 
the Third 

In the preceding paragraph it has been assumed that 
the food supply is ample and is independent of the 
presence of the consumers. If this is not the case, a 
third equation must be added. Let p r and ft z be the 
active numbers or biomasses of two predatory species 
and n the biomass of the prey species. The differential 
equations take the form : 

P\ = hi— e i + s i n - KiPi - KPzl ] 

P'2 = p2{-s 2 + s s n~h 21 p l — h 22 p 2 ) (12) 

n' = n(e — o^ — a 2 ft 2 — hn) j 

Note that here the co-efficients e v e 2 are no longer 
the same as those of equations (11). They represent a 
negative multiplication of the predators in the absence 
of their prey. 

For brevity the following notation will be used : 

A = hQ^h^— A]Ai) + s z {h n a z — \ z a^ 

E i = ~ £1(^22 + <V2) 4- e 2 (M 12 + ffasj 
-f e(— A 12 s a + AjjSi) 



E 2 = e^hhn + OiSg) — £ 2 (%i + ^1) 



(13) 



— e{— hjjs^ + h^Sj) 
E = — e^h^— hvfj) + fi 2 (/ r %o' 2 — }^a^ 
+ #iA a — W 

These expressions satisfy the following equations 

E i(/^ii + CiSj + E^M^ + a^J = fa— Ae x )A 
E i(^2i + CT i s a) + E a (M 22 -f <7 2 s 2 ) = (es a — fo 2 )A 
E^e + ff l£l ) -f E 2 (A 12 e + 0-363) = (es 1 — h&JE 
E i(^2ie + oi^a) + E 2 (A 22 e -f <r 2 e 2 ) = (ss 2 — Ae^E 



(14) 
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Five cases of stable equilibrium can be distinguished, 
that is to say, five possible orientations of the system. 

First Case. Suppose that 

E x > 0, E 2 > 



es 1 —he 1 > 0, es 2 — he 2 > o, h^h^ — h iz h 21 > 0} ' 
and therefore that A > o, E > o. 

These suppositions imply that the situation favours 
the two predatory species, but that neither of them 
holds a special advantage over the other. In this case 
stable equilibrium between the three components is 
possible and we get the following limiting values : 

E E 2 E 

lim fa = -!, lim f % = — , lim n = -r (16) 

*— >3o A f -*- do A <-*-:» ^ 

The limit has been reached without any fluctuation. 

Second Case. Suppose on the contrary that the 
situation favours neither species, that is to say that 

£S a — he-^ < o, es 2 — hs s < 0. (17) 

In this case the predators disappear and the prey multi- 
plies to saturation point. 

£ 

lim p 1 — 0, lim fi% = 0, lim n =— . (18) 

t -=► DO t — =► DO I -> DO 



n 



Third Case. Suppose that 

ss 1 — e-Ji > 0, E 2 < 

es 2 — eJi < 0, E x > o 



(19) 



In this case the first species is definitely favoured and 
finally survives, in equilibrium with the prey, while 
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the second predator disappears; we get, finally, the 
following limiting values : 

hmp 1 = —I—i-, hm p z = o, hm n = 1 ,*~[ (20) 

Fourth Case. In this case the second predator is 
favoured and we have 

es z -e 2 h>o Ex < o j , , 

es x — s x h < E 2 > o j v ' 

which in the limit gives 

hm p x == 0, hm ^> 2 = -7-^ — , hm n = a a *~ 22 (22) 

( -si- 30 * -> DO '*"2S ~I~ °2 5 2 ( ->• DO ""22 T a i s i 

Fifth Case. In this case we have the inequalities : 

E x < o, E 2 < 0, es t — fe x > 0, es 2 — Ae a > o (23) 

which are partly the reverse of (15). Everything de- 
pends on the initial values p^o), fi%[d), n(o). The 
system may move either towards the values of equa- 
tions (20) or towards those of equations (22). 

In comparing these results with those of the preceding 
paragraph it is seen that variation in the nature of the 
food effectively lessens the chance of equilibrium 
between the two species. The equations of the preced- 
ing section apply to optimum conditions. Equations 
(12), on the other hand, take into account the eventual 
changes in the food-supply. Hence the probability 
of equilibrium expressed by formula (16) is smaller 
than that expressed by formula (6) of the earlier 
section. On the other hand, according to inequalities 
(17) the probability of extinction of the predators 
reaches -25, while it is zero in the earlier section. 
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3. Two Species, One of which Eats the Other. 
Volterra's Biological Laws 

Consider two species of which the first, p v has access 
to abundant food, which is independent of the presence 
of the second species, p 2 ; while the latter is a predator, 
feeding solely on the species p x . Differential equations 
have already been obtained : 

P'x = Pi( £ i ~ KPx — hzPz) \ <„,s 

expressing the relation between these species. These 
equations have been considered from the point of 
view of the predatory species in order to determine 
the influence of food on the rate of increase of the 
consumers. Further, in the special case of the absence 
of quadratic terms, equations (24) express the relation 
between the animal and vegetable parts of the organic 
world. These equations will be reconsidered in order 
to study the influence of rapid variation in the vital 
co-efficients on the numerical development of the two 
components. 

Suppose, first, that there is no internal competition, 
either for prey or predator, that is to say that 

h n = b, h Z2 = 0. 

This is the cyclic case of Volterra 

P'l^Pii^i — ^izPi) \ hc * 

corresponding to a low density of both populations 
p 1 and p z . First write the integral : 

p l "p^e~ h '^~ h " p ' = H (26) 
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which has been mentioned more than once. Then : 
KP'i = E i J hiPi — ^AiPifiz 

thzP'z = ~ hhzPz + A iAvAA 
p\ 

— e 2 "7 — £ 1 £ 2 H~ e 2^12r2 

A 

p% 

Adding these relations we get 

KP\ + W« -^j-4 a = o 

Pi Pi. 

the immediate integration of which gives formula (26). 
H is an integration constant 

H ^fl^V- *""•-*'••' (27) 

where <h — pi{o), H = PM- 

In the graph of (p v p z ) curves (26) are closed and are 
without intersections (Fig. 4). Hence the functions 
p-iif), p 2 {t) are periodic, with a common period a>: 

M + »)=&(*). A (< + »)= M)- 

This periodicity is the subject of the first biological 
law enunciated by Volterra under the name of the law 
of the periodic cycle. It has been seen elsewhere that 
the co-ordinates of the common centre of the closed 
curves (26) are equal to the averages of the numbers of 
individuals of the two species during a period of time 
co: 



CO 



t + m />t + m 

tx*t = £,I*= T -\ P*dt=£- (28) 
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This fact is the subject of Volterra's second biological 
law, the law of conservation of averages. 

Having recalled these facts we pass on to Volterra's 
third biological law. Suppose that for some reason the 
mortalities of the two species increase. Equations (25) 
become 

P'x = A(«i — 0i — h 12 p 2 ) \ /- n v 

4>' 8 = &(-<*- a + Mi) J K9) 

in which 8 1 < e x and 6 2 are the complementary mor- 
talities. Let a*! be the corresponding period of these 
modified equations. Equations (28) assume the form 

ft+io x 



1 
w 1 

I 
co 9 



Pjdt = Ez + 



, 6 } (30) 

J.M £ 1~ d l 

t "21 



Volterra gives these relations the name of the law 
of disturbance of averages. It is seen that the comple- 
mentary mortality of the prey results in a decrease of 
the average number of predators. Similarly, comple- 
mentary mortality of the predators results in an increase 
in the numbers of the prey. 

Experimental verification of these laws is extremely 
difficult because periodic fluctuations corresponding 
to type (25) have never been observed. In most 
periodic examples may be seen the periodicity of relaxa- 
tion, the mechanism of which is very different from 
that of equation (25) and the laws of Volterra are not 
applicable to the periodicity of relaxation. Every 
direct confirmation is therefore extremely valuable. 
Statistical studies of fish have shown exactly the 
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indirect repercussion on the prey of an increased 
mortality of the predators, which confirms Volterra's 
third law. Gause's experiments on protozoa and on 
ticks only rarely gave closed curves comparable to (26) . 
In the case of a population composed of Paramecium 
bursaria (predator) and Saccharomyces exiguus (prey) 
it is possible, by artificially thinning the population, 
to obtain curves resembling the cycles of Volterra. 
In this way Gause has succeeded in obtaining three 
consecutive cycles. 

One of the causes of the rarity of truly periodic cases 
is the presence of quadratic terms. Limiting actions 
expressed by these terms always exist, especially 
among predators, and it is equation (24) that most 
accurately expresses the development of the system 
ifv Pz)- Two very different cases, which we have 
already examined, exist. In the^first case the inequality 

is equally favourable to both species ; for the prey it 
means that multiplication e x is rapid enough to coun- 
teract the destructive power of the predators ; for the 
predatory species it means that the decrease in num- 
bers, e 2 , is relatively small compared with the products 
of hunting. Thus as a result there is an equilibrium 
between the two components 

limfc ^ £ ^ + *A* , lim^ = _ e Ai-*Ai , . 

_ «-»■ DO «11«22 + %*u *-* »" M22 + Khz 
Finally we get 

]im &>%-,■ lim p 2< lL, ( 32 ) 

'-> DO "21 *->0O "12 
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that is to say the action of the limiting factors is 
advantageous to the prey but not to the predators. 

When the co-efficients h n and # 22 are small enough, 
the functions p x and p 2 tend towards the limits of 
equations (31), passing through decreasing oscillations, 
exactly analogous to those of mechanical systems with 
internal resistance. This analogy between limiting 
factors and internal friction is very interesting. In 
the language of dynamics which we have used in 
Chapter VII, the term — a n pp' in equation (32) of 
that chapter precisely expresses an internal resistance 
proportional to the population and to the speed of its 
growth. 

On the other hand, when h n and k w are large, the 
functions p x and p^ tend towards the limits of equation 
(31) without oscillation. 

In the second case we have 

which is unfavourable to the predators, so that they 
disappear while the first species tends towards satura- 
tion point : 

lim p ± = r— ■> lim p % = 0. 

1 _^. ^ /z u « -i» DO 

These simple results seem not to be in agreement 
with experiment. The final stable state (31) is but 
rarely observed. Among the causes of discrepancy may 
be mentioned the varying age of the prey and the 
preference of the predators for one or more of the 
age-groups. 

Effect of Age of Prey. Suppose that it is the adults 
which are chiefly attacked by the predators, while the 
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young- are better protected, either by their smaller size 
or by their taste or by living in a different station. This 
case can easily be expressed in equations. Let % x be 
the number of adult prey, # 2 the number of young prey, 
and p the number of predators. The following equa- 
tions may be written : 

x\ — — mx x + a.x 2 — h^Xjp \ 

x\ = tix 1 — (m -f a)x 2 (33) 

p' = — ep+ hgjfix! J 

instead of equations (25). The term ax 2 represents the 
number of young growing into adults, the term nx x the 
birth-rate proportional to the number of adults. 
A stable stationary state 

e ns a,n~ m z —■ ma, 
x i = 7—' %•> = 1 rr~ > P = — ; — ; ; — (34) 

will establish itself as long as the prey is saved from 
destruction by its fertility n and hy a sufficient survival 
of its young «, a condition expressible by the inequality 

an > m(m -f- a). 

On the other hand, if the mortality of the predators 
is great enough, 

{m -\- «) {cm — W 2 — ma) 



e > 



« 2 w 2 



the stationary state is a focus, or in other words the 
system passes through an infinite series of decreasing 
oscillations, very similar to the cycles of Volterra. 
Finally, if the mortality of the prey is too great 

m(m + a) > cm 

the two species ultimately disappear. 
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Hence, equations (33), which here replace equations 
(25) of the cyclic case, are less favourable to the future 
of the system and the influence of the varying age of 
the prey is expressed in a reduction of the chances of 
survival for either species. 

Heterogeneous Predators. Another cause may com- 
plicate the simple cyclic curve or the final equilibrium, 
and that is complexity in the predatory species. This 
species may consist of two forms of different sizes. If 
the food p 1 is plentiful the two races co-exist ; if food 
becomes insufficient the dwarf form is favoured. It 
requires less food and it is better able to catch the prey 
in its hiding-places. The prey being less severely 
attacked multiplies and so reproduces conditions fav- 
ourable to the giant form. This produces oscillations 
of a very different kind from the fluctuations of Vol- 
terra: it is a case that occurred in one of Gause's 
experiments. It goes without saying that the two races 
may differ in other characters. 

Didinium nasutum, Paramecium caudattim. In 
this experiment of Gause, the population is composed 
of two species of protozoa, one of which, Didinium 
nasutum, feeds on the other, Paramecium caudatum. 
The ferocity of attack is remarkable and in all Gause's 
experiments the result was the same : Didinium quickly 
exterminated Paramecium and then died of starvation. 
This result is not foreshadowed either by equations (25) 
or by the more general equations (24). In this latter 
case we had two types of equilibrium, neither of which 
corresponds to reality. How can this discrepancy be 
explained? Return to equations (24). The terms 
^hiPiPz an d ^ziPiPi. express the interaction between the 
two species. As long as the number of encounters is 
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relatively small, all is well, but if this number becomes 
large enough the interaction between the species can 
no longer be thus expressed. A simple example will 
show the correction to be applied to equations (24), 
If it is a question of meetings between males and females 
of the same species, the number of such meetings will 
be proportional to p 2 and the fertility should be ex- 
pressed by a quadratic term ; now it is in fact expressed 
by a term of the first order and the number of births 
in unit time is proportional to the population, not to 
the square of the population. Hence, if meetings 
become too frequent, statistical method introduces a 
term of the first order instead of a quadratic term. 
Here the extraordinary energy with which Didinium 
seeks and attacks Paramecium should lead to the same 
result, and equations (24) should be replaced by 

P\ = *iA - wh - hipi - KPii>* ] , -x 

P'z = — £ 2pS + NPz + KPlPz — KlPz i K ' 

as long as the prey exists. In the cyclic case the same 
correction would give 

P'i = E iPi ~ V1P2 ~ KaPiPa \ ,„ fi x 

P\ = - Ata - /*») + KPiPz J (3 ' 

Consequently, when e 2 > /- l z> we g& a diminished 
limit for the prey, 

t -$- 30 "21 "21 

and the influence of discontinuity in variation would 
lead to the disappearance of the prey, which the 
biological equations do not provide for. 
Gause has noticed in Didinium a certain dimorphism 
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leading to unexpected results. If the density of popu- 
lation of the prey falls below a certain limit, difference 
in size becomes an advantage for the small Didinium. 
The large predators disappear more quickly than the 
small ones, and this circumstance may sometimes save 
small Didinium from destruction. Under experimental 
conditions, in the absence of proper hiding-places the 
last Paramecium are quickly exterminated by the last 
of their small enemies and the whole population finally 
disappears. This does not occur in nature, where 
hiding-places are not lacking, and the result may be 
different. Gause tried to reproduce natural conditions, 
providing the Paramecium with a refuge inaccessible 
to the predators, but he did not obtain the expected 
result because the Paramecium did not leave their 
refuges and the predators died of starvation. The only 
way to obtain a relatively stable population was to 
supplement the two groups at regular intervals by 
suitably measured additions, and Gause thus obtained 
the desired result. 

To complete these results Gause undertook an inter- 
esting statistical experiment. He placed five Para- 
mecium and three Didinium in each of thirty identical 
tubes. Two days later he found the predators dead in 
four tubes and in twenty-six tubes a mixed population 
containing a number of Paramecium varying from two 
to thirty-eight. In another similar experiment he 
found, on the sixth day, mixed populations in seventeen 
tubes out of twenty-five. This result shows that the 
phenomenon is very complex, and that a distribution 
that agrees with theory in one of these microcosms 
cannot be used as evidence in favour of or against any 
theory whatever. 
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Bursaria truncatella: Paramecium bursaria. 
In another experiment of Gause's, Bursaria truncatella 
(predator) devoured Paramecium bursaria (prey). The 
latter were fed on yeast, regularly supplied. Because 
of the difference in size, the predators did not succeed 

in exterminating all 
the Paramecium and 
perished rapidly as soon 
as their food began to 
fail'. Then, in the 
absence of predators, 
Paramecium multiplied 
in the usual way (Fig. 9). 
This process is well 
suited to mathematical 
interpretation, provided 
the vital co-efficients 
are made to vary in con- 
formity with reality, 
and the result of discontinuity in biological variation 
is observed. Suppose that from the starting-point A 
(Fig. 9) to the point B the process is regulated by 
equations (24) or (25). From the point B the scarcity 
of food becomes appreciable and the few victims that 
some of the predators happen to catch are not sufficient 
for the whole group : from this moment the co-efficients 
h 21 and h 12 rapidly diminish, and at about the point C 
both co-efficients disappear, or in other words the last 
surviving Paramecium are no longer threatened by 
predators. From this moment equations (24) or (25) 
must be replaced by 



Fig. 



g. — Relative Growth of 
a Mixed Population 
(B. truncatella, P. bursaria) after 
Gause. 
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At about the point D there are no more predators. 
They practically disappear as soon as their numbers fall 
below a certain limit. The development of Paramecium 
then continues normally. 

Cheyletas eruditus; Aleuroglyphus agilis. Con- 
sider further the case of the mites Cheyletus eruditus (pre- 
dator) and Aleuroglyphus agilis (prey) studied by Gause 
and Smaragdova. A new factor appears here, the effect 
of metabolic products. This factor is the cause of an 
apparent paradox. Cheyletus placed in a population 
of its prey of low density can accustom itself to it and 
even establish a sort of cyclic process. The same mites 
placed in a dense population of their prey perish rapidly 
though it would seem that the abundance of food 
ought to be very favourable. In this case it is again 
possible to express the phenomenon in equations. 
According to Gause the smell of the prey is unbearable 
to its predators. At all events this individual odour 
of Aleuroglyphus does not prevent a Cheyletus from 
eating its prey if the surrounding atmosphere is fit 
to breathe, but robs it of the desire to attack when the 
personal atmosphere of the prey fills the whole of the 
space in which they live. It appears that in this case 
it is not a question of the accumulation of poisonous 
products but the absence of pockets of respirable air. 
Let p x be the density of the prey, p 2 that of the pre- 
dators, II 1; the value of p x marking the dividing line 
between a respirable and irrespirable condition for the 
predators. The process can be represented by the 
following equations : 

P\= A( e i - hipi — W> 2 ) ij. ^u)\ (*i) 
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The co-efficient E 2 represents the increased mortality 
of the predators in irrespirable air. The value p 1 = E^ 
cuts the graph (p lt p^ into two portions, and the result 
of the struggle for existence between these species 
depends on the original numbers of the two populations. 
Equations (37) and (38) give, broadly speaking, an 
explanation of Gause and Smaragdova's observations. 
It is seen that in many cases Volterra's equations 
are sufficient to explain the phenomena, but it is not 
possible to represent a natural process by a single 
formula which does not take account of sudden or 
gradual variations in the vital co-efficients. 

4. "Wolf — Goat — Cabbage" 

Populations of several species in chain-like relation, 
where each link is simultaneously predator and prey, 
are not uncommon. One case of this kind of population 
is that of the Island of Komodo in Malaya, inhabited 
by giant carnivorous reptiles, and by mammals — their 
food — which themselves feed on the rich vegetation of 
this island. We may write the differential equations 
of this biological system, assuming that competition 
exists only among the plants and that the presence of 
the reptiles has no direct influence on the vegetation. 
Let p x be the biomass of the mammals, p 2 the number 
of reptiles, and p 3 the biomass of plants. We get the 
following equations : 

P'x = Pi(- H - K2P2 + hzPa) ) 

P\ = Pi(~ H + V>i) (39) 

fi'z = Ps( e s — KxPx ~ haaPa) ) 
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Three different cases exist. 
First Case. Suppose that 

s zhJhi— £ As^3i— £ Ai^33 > (40) 

Stable equilibrium is possible with the following limiting 
values : 



lim^ 

t—> DO 

lim^V 
*-> 30 

lim^ 

f — > DO 



^21 
£ 3 «1 3 «2 



e %'hs"S. 



fi M 



&,& 



; 3"21 ' 



f 2"31 



/2n,&. 



21«33 



(41) 



Procedure to the limit may occur either with fluctua- 
tions or asymptotically. We do not stress the condi- 
tions that make it possible to distinguish the two cases, 
but note that fluctuations occur when the co-efficient 
h 31 is sufficiently small. 

Second Case. Suppose that 



£ 3%3"21 tyhsrSL' 



°1 



(42) 



In this case the food absorbed by plants is sufficient 
to ensure the existence of a reduced number of mam- 
mals, which do not, however, suffice to support the 
reptiles. Thus we get the following limiting values : 

e 3^i3 e i"3a 



lim p x = 

(-> DO 

lim pi== 
t— > DO 

\imp z = ^~ 

1 ~> do r ha 



faJi. 



13"31 



(43) 
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or, in other words, the reptiles disappear and the 
residue, (p v p a ) tends towards a stable limiting state. 
Third Case. In this case we have 

e s h 13 —E 1 h 33 <o; (44) 

the food supplied by plants is not enough to feed the 
mammals, so that 

lim p x = 0, lim p z = o, lim p z = — ■ (45) 

i -> CO * -> 30 t -> 30 H 33 

Thus all animals disappear and the remaining plant 
population develops alone to its maximum. 

What is the probability of these results? To judge 
from the sign of expression (40), the probability of 
equilibrium between the three groups is equal to '5, 
and that of the two others together is equal to >5. 
These figures may explain the relative frequency of 
three-sided equilibria in nature, but they do not give a 
reason for the disappearance of these reptiles from all 
parts save Komodo. This may mean that the estima- 
tion of probability in this particular case was not 
accurate. In fact, in this argument we have implicitly 
assumed the reciprocal independence of the vital 
co-efficients, which is not always true. On the other 
hand, the appearance of new factors is almost always 
unfavourable to the equilibrium of a system. 



CHAPTER X 

SYMBIOSIS AND PARASITISM 

i. Symbiosis 

The phenomena of symbiosis and parasitism, interest- 
ing and important as they are in themselves, allow us 
to obtain a better understanding of the mechanism of 
the evolution of species and the development of 
organisms. Without sharing the views of several 
eminent biologists (Famintzine, Merejkovski, Portier) 
who consider all living matter to be an instance of 
cellular symbiosis, it is difficult to deny, for example, 
the symbiotic origin of multicellular organisms. An 
association of two species, fortuitous at first, may 
become closer and closer. All possible degrees are 
known, from a very loose association easily broken 
to a practically indestructible union in which the 
partners form an integral whole that cannot be resolved 
into a mere juxtaposition of two components. This is 
true of lichens. Admittedly dissociation is still a possi- 
bility under laboratory conditions, but these conditions 
have nothing in common with natural circumstances. 

Problems of this kind may be put into the form of 
equations, not with a view to obtaining numerical 
results, but in the hope of qualitative conclusions 
which may have a certain value. Two fundamental 
types of such equations will be considered. 

Groups of Simple Structure. Take first the case of 
two species of a simple structure, with vital co-efficients 
almost independent of age. Let x v x 2 be the numbers 
of the free individuals of the two ultimately symbiotic 
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species, % the number of symbiotic couples, n x and n 2 
the fertilities in the free state, v x and v 2 the fertilities 
in the symbiotic state, m x and m 2 the mortalities in the 
free state and jjl x and /x % the mortalities in the symbiotic 
state. It may be supposed that a certain number of 
couples, vx, are born already associated, and that a 
complementary mortality p. of symbiotic couples must 
be added to the specific mortalities. Moreover, in 
certain cases the death of one partner does not involve 
the death of the other. Besides these factors, which 
give terms of the first order, meetings between free 
individuals must be considered, producing a.x t % 2 pairs 
in unit time. Finally, limiting factors must be con- 
sidered. In these circumstances, the following differ- 
ential equations are obtained : 

x' x — e x x x + [v x — v + f$ x )x — h xx x x 2 — h m x x x — ux x x 2 \ 
tf'a — e 2 # 2 + (v s — v + @ 2 )x — h w x<?— h^x — ctx x x 2 I 
x' = {v — n x — p 2 — ft)x — h^x* — h 31 x x x — #32*2* I ^ I 

In these equations each term has an exact biological 
meaning. The co-efficients are constant but do not 
possess secular stability and vary with time, sometimes 
slowly, sometimes rapidly. Before discussion of equa- 
tions (1) the history of a symbiotic union must be 
passed in review a history which, in the present case, 
finds expression in a variation of the vital co-efficients. 
The story of an association begins with the close 
approach of two species, originally quite independent 
of each other. 

At the outset, equations (r) have a simpler form .' 

x x = e x x x — • h xl x x hjgXjXz 



X 2 - — £'2^2 ' 'l%i%\%p '^22^2 



to) 
■* 9 — fc a-*a "asc'a 
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Of course the co-efficients of these equations are 
different from those of system (1). The two species are 
never connected and their inter-actions are virtually 
non-existent : 

h.% ~ °> Ki ~ °- 

In any case, for these co-efficients negative values 
are more probable than positive values. With this 
simplification equations (2) become: 

■ 2 === ^2^2 "22 

Suppose that these living objects inhabiting the same 
environment can fortuitously enter into very loose 
associations, no more, perhaps, than a mere approach 
that is advantageous to both. The chances of such 
meetings are not zero, and any event of which the 
probability is not zero is bound to occur.. Hence by 
the side of free individuals will appear a certain number 
of couples, the chance partners in which find an advan- 
tage in their association. Equations (3) will become : 

■ % 1 - £1 Ai -\- fJ\% '11 1 ' 1 3*1 GCiv^A'R 1 

^2~ a 2,%2 ~r P2* "22^2 "23*2''* OW^j > (4) 

x' = — If^ + /h)x~ h^— h^x% x ~ h 32 xx 2 + oa-^J 

The terms $ x x and fi z x correspond to the eventual 
separation of the couples. The term xx x x 2 expresses the 
formation of couples; the limitative co-efficients hj 3 , 
A 23 , h 31 , and h 3Z express the inter-actions between couples 
and free individuals. The association has advantages, 
and hence the limitative co-efficients satisfy the follow- 
ing inequalities : 

#•33 < hi <Kx<hz \ (5) 

"33 *■■- '^32 *■■- ^22 ^ ""23 
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In the initial stage (3) there is a state of equilibrium : 

e l e 2 

In stage (4) the equilibrium satisfies the following 
inequalities : 



% "22 



h < t^-' % < — ' x i <■ 



This implies that the formation of the third group 
(the couples) results in a decrease of the two free groups. 
On the other hand, too great a decrease in the numbers 
of free individuals is not an advantage to the symbionts. 
However, since the association is admittedly advan- 
tageous, the working of competition and of elimination 
may lead to the disappearance of free individuals and 
in consequence to the disappearance of symbiotic pairs. 
In these circumstances two solutions are possible. 
First, the symbiotic condition may become hereditary. 
The formation of a symbiotic pair becomes possible in 
what may be called the common domain of a couple. 
Some of the gametes of one of the associated species 
may already contain the gametes or even the adults 
of the other. In the case in which this phenomenon is 
in partial existence the process is regulated by equa- 
tions (1) which follow chronologically upon the ex- 
pressions (5). The next stage is the total disappearance 
of the free groups. System (1) is then reduced to the 
logistic equation : 

x' = (v — p)x ~ hx z (6) 

This is the case with regard to lichens. The wood- 
eating Termites which contain symbiotic flagellate 
Protozoa have not reached this stage, but they have 
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lost the ability to live independently. Other examples 
of this kind of evolution might be indicated. 

Groups of Heterogeneous Structure. We shall now 
consider the case of two different species, choosing 
among many examples the association of Eupagurus 
and Adamsia. The larvae of these animals lead a free 
life in different circumstances from the adults, and 
develop in different oecological zones. Let x be the 
number of symbiotic pairs, x n and x 12 the respective 
numbers of larvae and adults of Eupagurus, x iX and x n 
the corresponding numbers of Adamsia. It must be 
supposed that in the free state the larvae are not in 
competition with the adults and that there is no inter- 
action between the two species. 

First Stage. The first stage is expressed by the 
following equations : 

(?) 



% 11 — w l%2 ~ 
X 12 = nj^n 


T\ % 12 — gl X U 


x 21 = 7^22 
% 22 == ""$11 


r 2^22 £s#22 



(8) 

Equations (7) refer to Pagurus and equations (8) 
express the development of the Adamsia population 
before the formation of the symbiotic couples. A state 
of equilibrium between the two groups, each composed 
of one species, is theoretically possible but is never 
realized. Pagurus is a mobile animal, ill-protected 
against its many dangerous foes. It provides itself 
with a movable shelter in the shape of the shell of a 
gasteropod mollusc. Adamsia is an animal unable to 
move by its own power, feeding on prey brought to it 
by the current and provided with acontia as a means of 
defence. By chance encounters, the same shell may 



(9) 
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protect the Pagurus and carry the Adamsia. It must 
be supposed that at this initial stage the acontia are 
no danger to Pagurus, or otherwise the formation of 
the association would be impossible. The two partners 
in the association have an advantage over their free 
brethren. Their mortalities and their limitative co- 
efficients are certainly diminished, and equations (7) 
and (8) by the addition of compensating terms and a 
fifth equation taking into account the formation of the 
couples would then assume the form — ■ 

X 12 ~ "l^ll " r %12 Sl X 12 "13^12-*' a;v 12-*'22 

% 01 - — • ^0^22 — ~ //po^voi — '^2 21 2 21 "" I ^2 

X 22 ~ "2#21 T 2, X 22 §2*22 ^23^22* a 12' t '22 

X ' = — (/"l + ^2 + fA X ~ h^ — h\V 
»32^22' 1( ' ~T V-%\$"V2.- 

These equations represent a state of affairs created 
and maintained by purely fortuitous meetings during 
the search for shelter by Pagurus and for resting places 
by Adamsia larvae. There has been no pre-adaptation 
of the two species to a common life. A few individuals 
have merely been favoured by chance, and nothing 
more. 

Second Stage. This state might have persisted 
indefinitely but for the working of selection. In some 
instances selection eliminates fresh mutations, in others 
it removes the less favoured ones. The latter eventuality 
seems to be the more probable. Hence the working 
of selection increases the fortuitous advantage of the 
commensal pairs. This process may be illustrated by 
several examples of Pagurus in more or less close 
symbiosis with anemones. For example, the case of 
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Eupagurus prideauxii and Adamsia palliata represents 
the highest point in this evolution. The reciprocal 
adaptation is very close. Adamsia, very sensitive to 
external irritation, nevertheless offers no resistance to 
its capture by a Eupagurus prideauxii; it is unable 
to live if separated from the hermit-crab, and the 
latter, if deprived of its companion, does all it can to 
restore the symbiotic condition. 

To understand the working of the evolutionary 
mechanism, look closely at equations (9). Suppose 
that a state of equilibrium is nearly reached. The 
values of the variables satisfy the following inequalities : 

"■33 

It follows that each decrease in ar 12 % results in a 
decrease in x. The struggle for existence among the 
hermit-crabs is very keen, and it would seem that any 
advantage acquired by the symbiotic pairs must result 
in the total destruction of the species. This is where the 
second solution, of which we have spoken, comes into 
consideration. The existence of the larval phase of 
the two species explains everything. The appearance 
of mutants with tropisms resulting in mutual attraction 
increases the co-efficient a and consequently increases 
the proportion of couples in this mixed population. 
At the same time, the average age at which association 
occurs falls in both species. Free adults become rarer, 
especially adults of mature age, and the terms %%, 
« a # 22 , expressing the number of larvae proceeding from 
free groups, become negligible in comparison with the 
terms n^x, n^x, which express the numbers of larvae 
of symbiotic origin. On the other hand, the number of 
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couples formed in unit time becomes more and more 
proportional to the product x n x 12 , for young animals, 
after having emerged from the larval state, either find 
symbiotic shelter or very rapidly disappear. In this 
way equations (9) are simplified : 

x\ x = n x x — m 1 x n — h^x ix — h^x-^ 

x' 21 = n 2 x — m 2 %. — & 2 % ~ K x zi } ( I0 ) 

x 22 ~ b% X 21 — T 2, x 22, — El x Vl "iiX^X (XX-^X^i 

X ' = - (/"l + N + M) X — k 33 x2 + «*11#21 

If the co-efficients of fertility n v n z are large enough, 
equations (10) lead to a stable limiting state. Among 
other species of hermit-crab the association with 
anemones is less intimate than that of Adamsia -palliaia 
and Enpagurus prideauxii. Thus the anemone Sagartia 
associates with a hermit-crab but is also able to live 
without it. 

2. Simple Parasitism 

Parasitism may, in general, be treated as a special 
case of symbiosis by modifying the hypotheses con- 
cerning the vital co-efficients in such a way as especially 
to favour one species while maintaining the life of the 
other. But parasitism profoundly changes organisms 
and so deserves special study. 

Example 1 : Entomophagous parasites. Following a 
series of notes by W. R. Thompson on the action of 
entomophagous parasites, A. J. Lotka completed that 
writer's analysis with a system of differential equations 
of the first order. 

Admit, with W. R. Thompson, that two species be 
found together: the host, a noxious insect, and the 
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parasite, introduced to destroy the insect. A parasite 
lays on the host a certain number of eggs of which, on 
an average, k hatch and produce larvae, which in 
turn are transformed into free-living adults. Each 
oviposition of the parasite leads ultimately to the 
destruction of the infested host. Let p x be the number 
of uninfested hosts and ft z the number of free adult 
parasites seeking an insect host on which to deposit 
their eggs. The development of the two species may 
be expressed by two equations of the type — 

P'x = hPi ~ h nPi - VA ] (JJ) 

P'i^kKP\Pi~-VhPi-KP£ J 
The squared terms represent, as always, the total 
action of all the limitative factors. The term kh 12 pip z 
expresses the total number of parasite's eggs hatched 
out in unit time. Equations (n) permit three stationary 
states. 

1. Destruction of the two species. 

pi = o, p 2 = ; 

this state is stable only where e x is negative. 

2. Repulse of the parasites. 





Pi~T-' Pa = > 

"11 




this state 


is stable only when 

ke-Ji 12 < m z h n . 




3. Go~existence of the two species. 




Pi- 






this state 


is stable only when 

kejiw > m 2 h u . 
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We shall discuss these results from the practical 
point of view. Theoretically, the introduction of para- 
sites is not sufficient for complete destruction of the 
host. In the most favourable case the system {pift 2 ) 
tends towards an equilibrium value which is not zero. 
However, by choosing a sufficiently large original 
number of parasites the curve can be made to approach 
as closely as is wished to the p s axis and so to obtain 
total destruction of the hosts. 

hi practice, it is essential to take into accimnt the 
interdependence between the numbers of the hosts and the 
fertility of the parasites. In the experiments of Smirnov 
and Kusin on Mormoniella vitripcnnis, a parasite on 
the fly Calliphora erythrocepliala, a number of flies 
varying fi-om one to a hundred were put at the disposal 
of a pair of parasites. Each experiment lasted for an 
average of twelve days: the results are given in the 
following table : 



Number or 

NYMWIS 


I 


2 


3 


5 


10 


25 


50 


0+ 


98 


No. of parasites 
pernyoipli . 

Total no. ol 
parasites 


22-G 
33-6 


S3'8 
47'5 


67-2 


22'3 

1H'4 


i8'8 
187'5 


i.yi 
352-5 


9'0 
45r 


107 
08s 


7-8 
745 


Similar results were obtained with the common house- 
fly: 


Number of 
Nymphs 


I 


2 


3 


5 


10 


25 

5-7 

r+2 


50 


8.t 


So 


No. of parasites 
perilymph . 

Total no. of 
parasites 


12'+ 

IS'4 


8-3 
i/-5 


7-5 

22*6 


7-3 
3b-5 


6-4 
6>4 


5-4 
368 


3'5 
212 


4 
354 



It may be asked to what extent this multiple infesta- 
tion contributes to the destruction of the flies. The 
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result depends on the host. The proportion of surviving 
flies, which is very small for Calliphora, is considerable 
for Musca domestica. Hence total destruction of flies 
does not occur and equations (11) of Lotka need serious 
revision. 

In the following examples we shall study the action 
of parasites on the vital co-efficients and on the fate 
of their hosts. 

Example 2. Consider a case in which parasitism is 
maintained by the co-existence of healthy and infested 
hosts. Suppose, first, that multiple infestation does 
not occur and that the parasite does not multiply 
within the host; secondly, that infestation is due to 
free larvae coming from without; thirdly, that the 
number of parasitic larvae is proportional to the 
number of infected hosts ; fourthly, that the parasites 
do not influence the vital co-efficients of the hosts; 
and fifthly, that the hosts are born healthy. This gives 
very favourable conditions for the hosts. 

Let x be the number of free hosts and y the number 
of infected ones. Write the equations : 

%' = ny + ex — hx 2 — hxy — <xxy \ , , 
y' — ~~ my — hxy — hy 2 + axy j ^ ' 

The term cuey expresses the number of healthy hosts 
infected in a unit of time. Note that the whole popu- 
lation obeys the logistic law : 

*' + / = s(x + y) - h(x + y)K (13) 

Two possibilities exist for the distribution of this popu- 
lation between the two groups. In the first 

, m e 

a. >h, r < - 

a—- h h 
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which implies that the production of parasites is large 
enough to ensure the maintenance of the group of 
parasitized hosts by the side of the healthy group. 
Then at the limit : 

n ,. eoc— iih 

hm x =— . limy = 



* — > 30 K < — > 30 



ah 



In the second case 



hm 
h < ■ — » 

E 



which implies that the reproduction of the parasites 
is not sufficient to keep the infected group from ex- 
tinction. In the limit : 

lim x = y> lim y — 0. 

1 -* DO ' l 1 ->- DO 

Hence, in the first case the two groups co-exist, in the 
second the population rids itself of the parasites. 

Example 3 : Parasitic Castration. Suppose that the 
hypotheses of Example 2 remain valid save that con- 
tamination renders the host incapable of reproduction. 
Equations (12) are then to be replaced by others : 

x' = ex — hx* — hxy — axy ) , , 

y' = — my — hy z — hxy -f axy j * *" 

which differ by the absence of the term ny representing 
the birth-rate among infected hosts. This hypothesis 
does not alter the final result : the two possibilities of 
the preceding paragraph remain. In the first case 

as — hn > 0, 

which means that the reproduction of parasites is 
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sufficient to ensure the maintenance of the infected 
group by the side of the healthy group. At the limit 

m 1m , . eoi — hn 
hm x = — I — j> hm y — — 

* — > OO a K I -»- 30 x 

In the second case 

ae — hn < 

which means that the multiplication of the parasites 
is insufficient and that the process ends with the total 
disappearance of contaminated hosts : 

lim x = -> lim y = 0. 

t —> do " t -» 30 

Example 4 : Parasitic Action Increases the Mortality 
of the Host. Suppose again that the hypotheses of 
Example 2 remain valid, save that the mortality of 
infected hosts increases. The equations then assume 
the form : 

%' = ex + ny — hx* — hxy — a.xy ) , . 

y' = — [m -f- fj)y — hy 2 — hxy -f- a.xy ) ^ ' 

The qualitative result of the discussion does not differ 
from Examples 2 and 3. In the case in which 

e(a— h) > (m -f- f.i)h 

parasitic reproduction is sufficient to ensure the main- 
tenance of the parasitized group. In the case in which 

e(a — h) < (m + ju)h 

the parasitized group disappears, and at the limit 

lim x — -> lim y — o. 

t->2Q " *-> x 
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Example 5. Parasitic Action Increases the Limita- 
tive Co-efficients without Aflecting the Intrinsic Mor- 
tality and Fertility. In this case the equations are 

%' == ex -f ny — J^x 2 — h^xy — a.xy | , ,. 
y = — my — h zl xy — h w y z + axy j 

These equations are more complicated than the pre- 
ceding ones, but the result is the same. There are 
always two results. In the case in which 

the co-existence of the two groups is possible, and the 
limiting values of x and y are not zero. On the other 
hand, when 

e(« — h zl ) < mh u 

we have in the limit 

lim x = ~-> lim y = 0. 

t —. > do 'hi t—> 30 

In these four examples we have considered the possible 
influence of parasitism on all the vital co-efficients. The 
result is the same in every case. Parasitic activity is not 
sufficient to exterminate the hosts directly, but it weakens 
them both physically and numerically. This weakening 
may disturb the equilibrium which was possible between 
the given species and its environment, and may lead to 
its disappearance if its numbers fall below a certain 
critical value. 

3. Multiple Parasitism 
Multiple parasitism may arise in several ways — 
either by multiple contamination or by multiplication 
o'f the parasite within the host. The two cases have a 



SYMBIOSIS AND PARASITISM 159 

different mathematical aspect, but it is not always easy 
to distinguish them, and the simple criteria that are 
often recommended are not very certain. 

Multiple Contamination. Suppose, first, that multiple 
parasitism is due to multiple contamination. Let x 
be the total number of hosts, x k the total number of 

DO 

hosts with k parasites, p = J> hx k the total number 

k = 1 

DO 

of parasites, q — / x k the number of infected hosts. 
^— — i 

Suppose that the adult parasites are unable to lead a 
free life, but that their larvae live for a comparatively 
short time in the environment and penetrate the hosts 
as they chance to meet them. Choose the unit of time 
so that no host can be infected more than once per unit. 
In these circumstances the differential equations may 
be given the following form : 

DO CO 

p n s x s — m x — c px — x > h oa x, 

s = o * «=* o 

x' k = — m k x k — c k px k + c k - x px k - x — x k y h ks x s 



X n — 



(17) 



The contamination co-efficients c k are not the same 
for all the groups. When immunity follows contamina- 
tion the co-efficients form a diminishing series. In the 
case of parasitic enfeeblement of the hosts the co-effi- 
cients form an increasing series. The mortality of the 
parasites within the bodies of the hosts is neglected. 

Suppose that the biological system that is being 
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studied has reached a stationary state. Eliminating 
x' k we have, in general, 

7*- CJ ^~ < l8 > 



c k p + m k + J> h ks x s 



Suppose that the presence of the parasites weakens 
the hosts and diminishes their resistance to later 
infections. In this case 

and in these circumstances equation (18) gives 

which implies that the series 

x x-l * 2 ... (19) 

is diminishing, or, in other words, that hosts with 
k parasites are rarer than hosts with k — 1 parasites. 

Suppose, on the other hand, that the presence of 
parasites forms a bar to the entry of other larval para- 
sites, that is to say that 

C k-\ > C k- 

In this case it is qr^te possible that for certain values 
of k 

x k > X k~l 

and the series (19) loses its diminishing character. It 
is even possible that in some instances it may be 
increasing. Hence it is incorrect to state that the 
absence of reproduction of parasites within the host 
makes series (19) a diminishing one. 

Multiplication of Parasites within the Host. Con- 
sider now another case in which there is no multiple 



(20) 
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contamination and in which reproduction of the para- 
site is the sole cause of multiple parasitism. The 
differential equations now assume a slightly different 
form 

DO DO 

x' = y> n s x $ — m x — cpx — x \ li 0a x 3 

J n= S " 

DO 

x\ = — mi x x — Pfa + cpx a — % y h l3 x, 

s ^ o 

DO 

x' k = — m k x k — fi k x k + /W^-i ■- x k \h ka x s 

s ^ o 

The terms (3 k x k represent the reproduction of the 
parasite in the host. Here also two different results 
may be obtained. It is possible that the presence of 
many parasites in one host does not prevent them from, 
increasing (Malthusian multiplication). In this case 

h > /W 
Then in the stationary state 

*- *S — en 



x 



*-l 



Wfc -h h + jy h>x, 



so that series (19) is decreasing or in other words 
hosts with k parasites are rarer than hosts with k—x 
parasites. 

Suppose on the other hand that within the interior 
of the host the parasites are subject to the same 
limiting factors as affect any population. In this case 
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and it may happen that for certain values of k 

so that series (19) loses its diminishing character. It 
would therefore be incorrect to state that the absence 
of multiple contamination does not allow series (19) 
to diminish. 

These simple considerations show how much care is 
necessary when we wish to draw conclusions as to the 
nature of multiple parasitism from the numerical aspect 
of the distribution of parasitized hosts. If series (19) 
is diminishing a choice may be made between multiple 
contamination facilitated by an enfeeblement of the 
hosts and unlimited multiplication of the parasites 
within the host. If the series is not diminishing a choice 
exists between increasing immunity of the host and 
limited multiplication of the parasites within. It may 
be pointed out that simultaneous action of two factors 
may further complicate the position and make the 
choice impossible. 

We shall take a few examples which will allow a 
better appreciation of the influence of different factors. 

Example I. Internal multiplication does not 
occur. The parasites do not affect the vital co- 
efficients. 

Equations (17) assume the form : 

x' — nx — mx — cpx a — hxx (22) 

%\ = — mx k - cpx k 4- cpx k ^ ~ hxx h (23) 
(It = 1, 2, . . . ) 

The sum of these equations is equal to 

x' = (n — m)x — hx 2 (24) 
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which shows that all the groups together satisfy the 
logistic equation. On adding the same equations 
multiplied by 1, 2, 3, . . . , we get 

p' = p[— m + (c - h)x] (25) 

The solution of equation (24) is 

a(n~m)e { "~ m)l 



n — m — ha -|- hae in ~ m)t 
where a denotes the original value of x. We also get 



(26) 



P = P(°) e 



,— mt 



ha ha 
_____ j pi" - »•)' 



h~ (27) 



n — m n — m 
The total number of hosts tends towards the limit 

n— m 



lim x — 

As far as the distribution of the parasites is concerned, 
two possibilities have to be considered. In the first 
case we have 

c[n — m) > nh 

and the total number of parasites increases indefinitely ; 
it is clear that this state of affairs is biologically im- 
possible, for if parasites become more and more num- 
erous they end by upsetting their hosts, and the vital 
co-efficients will change in such a way that the problem 
will be completely different. In the second case we 
have 

c(» — m) < nh; 

the number of the parasites approaches zero, as does 
all x k except x , which approaches — 7 — , and as a 
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result the hosts rid themselves of their parasites. The 
inequalities 

c(n— m) ^ nh 

lend themselves to many biological interpretations, 
which may be expressed in any case in accordance with 
the peculiarities of the species under consideration. 

In the following examples we shall study the effect 
of the disadvantage due to the parasites, in the absence 
of any internal multiplication. 

Example 2. The parasites modify the co- 
efficients of mortality and of meeting the host. 

Suppose that parasitic action modifies the mortality 
and meeting co-efficients in the following way : 

m = m, m k — m -\- 8, (k = 1, 2, . . .) 

h 0s = X, h k , = X + ft, (It = 1, 2, . . . , s = o, 1, 2, . . . 

The co-efficients of fertility and contamination remain 
unchanged. In these circumstances equations (17) 
become 

x' = nx — m% — cpx — Xx Q x (28) 

x\ = — (m + 8)x k — cpx k + cpx k - 1 -~ [X + fi)xx k 

(k = 1, 2, . . .) (29) 

From these we get 

p' = p[- (m + 0) + (c - A - n)x] (30) 

x' = (n — m—6)x—(X + fi)x 2 + dx -f- /j,xx (31) 

The three equations (28), (30), and (31) enable us 
to evaluate the three functions x , x, and p. We here 
limit ourselves to the equilibrium state. Two possi- 
bilities are to be considered. In the first case we have 
the inequality 

c— X — n < 0; 
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the function p decreases and approaches zero, which 
means that the variables x v x 2 , . . . , also tend towards 
zero ; the host rids itself of parasites, and finally 

,. n—m 

hm x — lim x a — — - — 

* — > 30 * — > do *■ 

In the second case we have 

c — X— fi > \ 

(c — A — fi) {n — m) > Mm -f 0) (32) 

c(m -\- 6)- >n(c— A— /i) J 

Equations (28-29) give the following equilibrium 
values 

m + 



X 



X 



Xjfc = 



c — A — ju 

{m + Q)[(n - m) (c — 1 — p) — Urn + 0)] 
(c — X — /()[<:(m + 0) — «(c — I — ,u)] 
(w + 0)[c{ m + 0) — n{c— A — //)] 

(^T^^w + 0) + 0(c- A- ^)] 
c*P*X P*X 



[ m + + c p + (A + /iJX]* (P + X)* 



Thus the values of X k form a geometrically decreasing 
series, which, in certain cases corresponds closely with 
observation. The equilibrium is stable and the process 
tends towards it when the inequalities (32) are satisfied. 
In this case co-existence of healthy and infected hosts 
is possible. 

The following examples serve to show that in the 
case of a progressive increase in the co-efficients of 
mortality and the limiting co-efficients the only possible 
results are either the co-existence of the groups % 
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x v x z , . , '. , or the complete extinction of the parasitized 
groups. 

Example 3. The parasites progressively change 
the limiting co-efficients. 

To fix our ideas, suppose that the parasitic action is 
expressed by an arithmetical progression of the limita- 
tive co-efhcients : 

h ks = X + P&. 

Equations (17) become 

x' = nx — mx a — cpx — Xx x (33) 

x'k = — ««* — opx k + cpx k - x — (A + fyt)*** 

(* = r, 2, . . .) (34) 

These equations give first 

x' — (n — m)x — Xx % — fifix (35) 

Two stationary states are possible. In the first case 
we have 

[iP -+- AX = n •— m 

X, cP 

< 1 



X M w + cP + (A + A»»)X 

X "X 

m + cP 4- AX 

which enables X k to be expressed in terms of P 

r n(n — m~ pP) A*c*P* 
x k = - k 



[In + AP(c- p)] n \hx + A{c- n)P 

(A = 1, a, . . •) (36) 
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Replacing X 7i by these expressions in the equations 



DO 



p — / kxj. 

we get, to determine P, the equation 

30 
t— P) "ST -1 



30 

n(n —m—V) 

hi + X{c ■ 



_ (37) 

II [Xn -]- l(c — fi)V + sju[h — m — /^P)] 

The values of P with which we are concerned must 

0j ) for otherwise x k 

would have negative values, which is impossible. The 
right-hand side of equation (37) disappears when 

p — or when p — -*- — . Hence it passes through a 

maximum in the interval ( 0, ). Equation (37) 

has certainly one root other than zero in this interval 
if the inequality 

In < (c — fj) (n~m) (38) 

is satisfied. This inequality may be interpreted in 
several ways; it may be said, for instance, that the 
contamination of the hosts expressed by the co-efficient 
c is a factor which is opposed to the destructive ten- 
dencies represented by the limitative co-efficients A 
and fx. Suppose that the inequality (38) is satisfied 
and let P be the first root, other than zero, of equation 
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(37). The first stable stationary state of the system is 
expressed by the values 

n(n — m— /*P) 
^ ==P ' * 0= A»+ AP(c-^T) 

and by formula (36). The series x k is decreasing. 

Suppose on the contrary the inequality (38) is not 
satisfied. Then the system 

n— m 
p = o, x — x = — — , x x = x 2 — . . . = o 

represents the second stable stationary state in which 
the factors favourable to the parasites are not sufficient 
to ensure the existence of parasitized groups. 

Example 4. Progressive increase in the co~ 
efficient of mortality. 

Suppose that the co-efficient of mortality takes the 
form 

m k = m + Qk 

and that the other vital co-efficients are unaffected. 
Equations (17) become 

x' = nx — mx — cpx — hx x (39) 

x' h = — (m + kQ)x k — cpx lc + cpx^ — hxx k 

k = 1, 2, . . .) (40) 

The equilibrium state of this system is given by the 
equations 

dp = x(n — m — hx), (41) 

nx _ c k p k x 

m -f cp + hx 



v o = „, , „j. , ,,„ ' %k — h 



II (cp + hx + m + sO) 
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Substituting these expressions in the formula 

DO 

s = r 

we get 



P ~ m + cp + hx^ * 

ft = 1 11 (6- 



(42) 



(cjf) + /j# -\-m-\- sO) 



Equations (41) and (42) enable p and * to be deter- 
mined. In the case in which 

h(0 + n) > c(n — m) 

there is the stationary stable state 

x = x = — - — , ^ = 0, * fc — 0, {k = i, 2, . . .) 

In the case in which 

h{6 -\- n) < c{n — m) 

there is a stable state of equilibrium with p equal to 
the first root, other than zero, of equation (42). 

Thus, continued elimination of contaminated hosts is 
the only means of protecting a population against an 
indefinitely increasing contamination. In every case an 
infected population can rid itself of parasites by a suitable 
modification of the vital co-efficients. 

Example S. Pagurus and Chlorogaster, We shall 
end this chapter by a study of the only case of para- 
sitism for which we have exact figures. It is the case 
of the Pagurus parasitized by Chlorogaster. These are 
rhizocephaloid crustaceans, generally presenting several 
external visceral sacs on each infested Pagurus. Each 
sac has its own system of apparently independent roots. 
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C. Perez examined more than two thousand Pagurus 
and worked out a formula for the distribution of the 
visceral sacs on infested Pagurus. According to his 
anatomical studies the hypothesis of multiple con- 
■tamination may be disregarded, and there remains 
the hypothesis of reproduction of the parasite within 
the host. Here there is a choice between two possi- 
bilities, limited multiplication or unlimited multipli- 
cation of the parasites. Only the former can give a 
non-decreasing distribution of parasites. In the cases 
observed by Perez the distribution is not decreasing 
and hence the multiplication of visceral sacs is limited. 
The internal multiplication of parasites may occur in 
several ways : 

1. It may be supposed that the original infection was 
due to a single larva, but that the parasitic mass, after 
becoming internal, was subdivided precociously into a 
number of independent fragments, which have then 
developed simultaneously. 

2. It may be supposed that the parasitic mass 
produced by inoculation with a single larva remains 
itself undivided and gives rise to a continuous root 
system, a long-lived stock from which successive groups 
of external sacs are budded off by discontinuous growth. 
As a result of his most recent anatomical studies of 
contaminated Pagurus, Perez has adopted this hypo- 
thesis of budding. 

Admitting the hypothesis of budding, suppose that 
there is some relation between the number of sacs and 
the age of the host. Let « be the age of a parasite. 
This age is always less than that of the host, Thus, if 
the number of sacs increases with the age of the 
•parasite, the Pagurus which carry the greatest number 
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of sacs are on the average the oldest. Let N„ be the 
total number of contaminated hosts and N„^(«)rf« 
the number of parasites of age (u, u + du). This 
number decreases in accordance with the mortalities of 
Pagurus and of the parasites. Denote the combined 
co-efficient of mortality by x(u) ; then evidently 

*M=-Trv #«)=#>)*-/>* 

Now bring in the frequency function of the parasites 

N e ^(w)y(w, p) du dp 

in respect of age u and number of sacs p. This function, 
which gives us the number of parasites with (p,p + dp) 
visceral sacs having an age («, u -\- du) is wholly 
analogous to the functions used in demographic 
statistics. Let the function y{u, p) take the form 



*fc«-?5o v Uj} 



Denote by N(p)dp the number of Pagurus carrying 
(p, p + dp) sacs. Then 

N c = f^(p)dp, 1 = fv(p)dp. 

Jo Jo 

From this the integral equation easily follows : 



m \m, 

which makes it possible to determine the frequency 
function if all the other functions are known. In fact, 
we get 



V(4) = -2- 



p-"Y{z)dz 
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putting 

/•CO 

N{ft)ft" dp 

F(*) = ji~ ■ 

<j>{u)[p{u))" du 

J 

The functions /?(«) and %{u) are unknown. Give 
them an arbitrary but possible form 

x(u) ^ B?]^- \ p[u)=ku° 

On the other hand the empirical function 

vw^ m\P/ 



P / m 

represents very well the observed distribution of the 
parasitic sacs, P being the number of sacs to which 
the maximum of N(£) corresponds. As the following 
table shows, the approximation is sufficiently exact. 
The parameters taken are : 

Males n — 7 m = 175 P = 5-0 N(P) — 29 

Females n = 7 m = 175 P = 4-5 N(P) = 31 

Taking n as equal to - and m as there is no diffi- 

£ e 

culty in calculating the functions V(p) and y(u, ft). We 

get 

Q3(«)]- - • exp {- —^ |j8(«)]"^- 
Wm, *) = const 
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p 

N0>) 


N(«<J 


NW? 


N(p) Total 


Obs. 


Calc. 


Obs. 


Calo. 


Obs. 


Calc. 


I 
2 
3 
4 
5 
6 

7 
S 

9 

10 
II 


4 

7 
9 

24 

29 

25 

11 

10 

6 

3 
2 


2-3 

7-8 
157 

24-5 
29-0 
24-2 

15-5 
8-8 

5-i 
3.0 
1-9 


2 
10 
17 
3° 

29 

17 
7 
5 

2 


3-0 

I0-O 
20-O 

29-3 
29-3 

I9'5 
IO-8 

5'7 
3-2 
1-9 


6 

17 
26 

54 
58 
42 
18 

IS 
6 

5 


5-3 
17-8 

357 
53-8 
58-3 

43-7 
26-3 

14-5 
8-3 
4.9 



This is one of Pearson's frequency curves; quali- 
tatively it corresponds well to the probable distribution 
of Pagurus in respect of u and fi, if it be admitted that 
the number of sacs increases with age up to a certain 
limit. From this the distribution according to Perez 
follows immediately. Consequently the hypothesis of 
single contamination supplemented by a hypothesis of 
limited budding is well able to explain the facts, 



CHAPTER XI 

GROWTH OF ORGANISMS 

i. General Considerations 

By a natural and simple generalization an organism 
may readily be compared to an animal or vegetable 
population, such as has been studied in the preceding 
chapters ; the cells, and sometimes the organs, are 
considered as almost independent elements. It is 
indeed possible to make up a series connecting the two 
extremes, and it is very difficult to point to the exact 
place of separation between the two phenomena. 
However, this impossibility only means that the two 
extremes are blended and that the same laws govern 
the numerical increase of a population without social 
organization, composed of wholly independent indivi- 
duals, as govern that of a higher organism with many 
complex organizations. It may happen that, in a first 
approximation, the two may be represented by similar 
equations, but it is certain that neither a second 
approximation nor the interpretation of the vital 
co-efficients will be the same. There is not, in fact, a 
demographic analogy even between populations with 
different social organizations. A primitive society, 
almost anarchical and almost homogeneous, may be 
compared to a rarified gas. A modern society, capitalist 
or socialist, more closely resembles a very dense, 
heterogeneous, viscous fluid, with well defined currents, 
obeying a set of laws and forces unknown to the other, 
but the most important difference consists in the 

174 
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existence of a division of labour, and of special organs 
for social functions. In an animal population, each 
individual defends and provides for himself according 
to his needs, a cell in an organism shares in the life 
of the whole, receives food, and is relieved of its meta- 
bolic products. The food necessary for a population is 
proportional to the number of individuals; in an 
organism the food needed is proportional neither to 
the number of cells nor to the weight of the body. 
The same is true of all relations with the external 
environment; for an animal or vegetable population 
the exchanges are proportional to the number of indivi- 
duals, while for a warm-blooded animal they are 
proportional to its surface area. The capacity for 
reproduction in an animal or vegetable population 
belongs to all individuals,* but this is not true of the 
cells of an organism. In short, the difference between 
the two types is not a quantitative but a qualitative 
one, and the dialectic passage from quantity to quality 
is real, palpable, and undeniable. 

The Development of Multicellular from Unicellular 
Organisms. The successive stages in the evolution of 
multicellular from unicellular organisms can be recon- 
structed. The beginning of the process must have been 
symbiosis resulting from imperfect separation after 
conjugation. Incomplete separation reduces the total 
area of the two bodies and diminishes their exchanges 
with the environment. In some circumstances this 
reduction is an advantage, particularly among Protozoa 
with a specialized surface. The formation of chains, 
resulting from division unaccompanied by separation, 

* I exclude animal societies such as the beehive and the white 
ants' nest, which require separate treatment. 
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may sometimes be seen among the Ciliata. By a 
modification of the physico-chemical conditions the 
complete separation of a multicellular organism at the 
beginning of its formation may be artificially produced, 
Among a variety of environmental conditions it is 
certain that a favourable situation will be found in 
some more or less isolated spots. The work of selection 
has done the rest. In multicellular organisms thus 
formed specialization of cells is inevitable, and the first 
specialization is a differentiation between external and 
internal cells, or sometimes between upper and lower 
cells. At the same time retardation and localization 
of certain functions can be seen. Thus a Stylonychia 
undergoes about five divisions a day, but the multipli- 
cation of the cells of an organism is subject to a slower 
rhythm. Further, the capacity for growth is not found 
in the same degree in all parts of an organism ; among 
plants it is strictly localized, among animals localization 
is less pronounced but nevertheless it still exists. 

Function of External Surface. In the life of a popu- 
lation the frontier plays a minor part; it is in the 
inside that nourishment takes place, and in constructing 
equations one only takes account of the movements of 
immigration and emigrations, which give several com- 
plementary terms, the analytical nature of which is 
analogous to that of the fundamental terms. In an 
organism, on the contrary, immigration and emigration 
of cells can scarcely be conceived, while currents of 
nutritive matter and of energy are of first importance. 
This fact suggests that the fundamental expenditure 
in warm-blooded animals, that is to say the energy 
necessary to maintain the body at a constant tempera- 
ture is proportional to the area. This law of surface 
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area has been satisfactorily verified by experiment. 
But functional expenditure is not proportional to the 
area. Again, the utilization of, food in an organism is 
governed by special laws, at present imperfectly known, 
but which in no case result in a proportional relation 
between the food and the energy produced by its 
consumption. 

Constructive Expenditure. Besides the basal ex- 
penditure and the functional expenditure, the expendi- 
ture necessary for construction must be considered, 
that is to say the difference between the number of 
cells reconstituted and destroyed in unit time. The 
expenditure of energy necessary to account for this 
difference is very small in comparison with the other 
two, but its results are important. The constructive 
potentiality of an organism depends, first, on glandular 
activity, secondly, on the accumulation of metabolic 
products in the body, thirdly, on the local retarding 
action, or in other words, the density of the cell- 
population, and fourthly, on the localization of cell- 
division. All these factors vary with the age of the 
organism. The organs of internal secretion do not 
come into activity simultaneously. Food varies with 
age, both as regards quality and quantity : in mammals 
the period of lactation is very different from the rest 
of their lives. The intervention of the sex glands 
profoundly modifies the vital co-efficients of an organ- 
ism. The accumulation of metabolic products with a 
consequent degeneration of tissues leads first to a 
retardation of growth and then to the dissolution of 
the organism. Finally, in addition to these accelerating 
or retarding factors the external form of the organism 
must be taken into account. The development of form 
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in its turn obeys laws which are doubtless simple but 
which are hitherto unknown. A few statements of fact, 
chiefly of a negative character, are of interest. First, 
the growth of an organism does not occur simultane- 
ously in all directions. The old observation that a 
man grows first in height, then in breadth, and finally 
in "bulk" applies mutatis mutandis to all organisms. 
It follows therefore that at different times in life growth 
is not proportional either to weight or to area or to 
stature. Is it proportional to some power of the weight 
w m 7 If so the exponent m depends on the age and 
perhaps also on the weight, 

m ~m(t, w). 

The same is true of all other biological characters 
and factors. 

Growth Problems expressed as Equations. The final 
result of this discussion may be summarized as follows : 

1. The development of an organism or even the 
growth of any bodily character cannot be represented 
by a single equation or a single formula. 

2. The process must be divided into several parts (for 
instance, for man; embryonic stage, breast-feeding, in- 
fancy, adolescence, sexual maturity, old age) and each 
part must be represented by the appropriate equation. 

3. It may be supposed that numerical characteristics 
of any part of the body whatever are linked with the 
weight w by a relation of the type 

4. Hence each stage of increase of weight can be 
represented by the integro-differential equation 

w' — e„w l '' — h i w t >— c s w m> j w"'(u)du (t s _ 1 <t<t,) (1) 
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with a suitable choice for each stage of the co-efficients 
e„ h„ c 3 , and the indices k c , l s , m a , n ss and also of the 
initial and final moments. 

5. In this equation the term e s w hi expresses the whole 
of the constructive activities, the term h,w 1 ' expresses 
the whole of the retarding activities and the integral 
term expresses the residual activities. 

This form of the equation is not very precise, never- 
theless it expresses the antagonistic influences which 
act upon an organism during its growth. 

Insufficiency of Biological Statistics. A few more 
words are necessary on the numerical data of biology. 
In all biological statistics serious methodological defects 
are to be found. Thus in data concerning the growth 
of animals, each average is deduced from a different 
number of animals, which makes some publications, 
full 'of figures and very expensive, altogether useless. 

A deplorable error is to be found concerning feeding 
and metabolism in animals. Express the weight, the 
metabolism, and the food consumed by an animal each 
day as a function of its age :. 

Weight —f(t) ; metabolism = </>(t) ; food = y>(t). 

Eliminating t from any two of these we get equations 
of the type 

Metabolism = ¥[w) ; Food per day — <b(w). 

Again, from observations of animals of the same age 
but different weights, relations between the same 
biological variables can be deduced. 

Metabolism = F 1 (ze') ; Food per day = ^(w) . 

but this does not mean that 

F{w) ss F t (w). <b[w) ss Ojiw). 
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In reality these functions are completely different. 
This circumstance is often forgotten. The true relation 
between weight and metabolism of an animal is given 
by the function F x and not by F. 

Growth-curves of animals possess a truly remarkable 
peculiarity. Save for man, these curves never pass 
through a maximum. This means that only the ascend- 
ing portion of the life of an animal has been observed, as 
if there were no aged animals in existence. It seems that 
this is due to the extreme vulnerability of animals in 
natural circumstances. Hence, contrary to the general 
opinion, natural conditions are wholly unfavourable 
for animals, and it is to civilization that man owes the 
descending portion of his growth curve. 



2. Utilization of Food by an Organism 

A simple relation can be deduced between the food 
absorbed by an animal and its utilization. This relation 
gives a plausible explanation of the empirical law 
deduced by S. Brody and R. C. Procter. It concerns 
the relation between the net energy % taken in by an 
organism and the gross energy X absorbed as food. 
According to Brody and Procter this relation is ex- 
pressed by the following formula : 

* = A(i-e-" ,x ) (3) 

I shall show that this formula can be replaced by 
another equally simple and corresponding equally well 
with the observations, but based on a few almost 
evident premisses which permit a wider range of 
application. 

Suppose that the utilization of large and small 
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amounts of alimentary energy do not obey the same 
laws and that the two gross amounts of energy A and 
B which separately produce the net energies a and b 
give, acting in conjunction, the net energy 

a 4- b 

<a + b. 



i + k 2 ab 



Denote by the sign (+) the operation of physio- 
logical synthesis so that 

If A and B are very small, we get approximately 

A(+)B ~ a + b (4) 

and the net energy a derived from the gross energy A 
is proportional to A 

a = oA (5) 

In the same way the symbolical operations of sub- 
traction, multiplication, and division can be interpreted. 
Multiplication is the only one that concerns us here. 
We find easily 

A(+)A = 2(x)A 2a 



i + kW 

23 

+ a 






2 



i + k 2 a 

put " <x)A = Sw (6) 
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It is easy to see that the polynomials P„(«) and Q„(«) 
satisfy the following relations 

V» + *(a)=-p n (a) + aQ n (a) 

and also that 

P n («) = - P,(- a), Q„(«) = Q„(- a). 

Equations (7) give 

kP n (a) + Q» = (1 + ka)» 

which, taken together with (8) permits P n (a) and 
Q„(a) to be evaluated 

(1 + ka) n - (1 - ka) n 



(7) 
(8)' 
(9) 



P„(«) = 

Q» = 



2^ 
(1 + te)« 4 (1 - ka) n 



and we get finally 
«(X)A = 



(1 + £a) n — (1 — A«) n 



(10) 



(11) 



k[(z + ka) n + (1 - *«)"] 
Suppose that A and a are small and n is large. Replace 

a by crA, A by — . w( x ) A by x. The preceding formula 

becomes 



= - lim 

K 
»->■ DO 



v^r- 



foxy 



( 



. kaX \ n . 1 
i + — j +[i. 



*aX\" 



or 



1 6 

x — - 



kaX 



-haX j 

zm = - k th(kaX) 



1 JkaX~J^TaX — J >>">\™^) (12) 

When k is small enough equation (5) reappears. 
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The maximum possible value of the net energy is 

-■ The form of curve (12) is the same as that of curve 
n 

(2) and curve (12) represents the numerical data of 

Brody as well as curve (2) . 

One important point remains to be cleared up — the 

dependence of the co-efficients k and a on the weight 

of the animal. It may be expected a prion that the 

digestibility of food will not be the same for animals 

of different weights. Brody's figures and equation (12) 

suggest a solution of this problem. The following figures 

relating to the cattle of Illinois (111. Steer B in Brody's 

notation) exist : 

Average weight = 613 kg., - = 38,500, ak = -0000126 
which give 

The figures for another breed (Steer 36 PA) are : 

Average weight = 433 kg., - = 24,000, ak = -0000195 
whence 

° = ^ h = 56 - 

This result is most interesting. The co-efficient a is 
a physical co-efficient which is constant for the different 
breeds of a species, and its constancy is an encouraging 
confirmation of our theory. It is further seen that the 
maximum assimilated energy is proportional to the 
weight of the animal, and that therefore : 

I = aw (13) 

which in a general sense confirms the observations. 
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In these conditions, formula (12) becomes 



x = awth ( — (14) 

\ aw I 

Two extreme cases may appear. In the first — 

is so small that the hyperbolic tangent can be replaced 
by its argument, which gives the relation between x 
and X 

x~<rX; (5') 

and in this case the weight of the body does not ex- 
plicitly appear, and the net energy depends on w only 

by the intervention of X. In the second case — is 

a.w 

so great that the hyperbolic tangent may be replaced 
by- 1, which gives 

x~«.w (15) 

that is to say that the net energy obtained by an 
animal from a very abundant food-supply is propor- 
tional to the weight of its body. The use of these 
formulae in the theory of the growth of animals will 
appear later. 

3. Growth of Tissues 
B. Ephrussi and G. Teissier have carried out experi- 
ments on the growth of tissue cultures, concerning the 
residual growth of fragments of tissue in a purely 
protective medium. According to Carrel the residual 
activity of tissues arises from the small reserve which 
every fragment of tissue taken from an organism 
carries with it. Among the factors concerned in the 
growth of these cultures are two to which the authors 
attribute a decisive part. 
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1. The gradual exhaustion of certain substances. 

2. The peripheral localization of the majority of the 
mitoses. 

Let s be the area of the culture at the moment t 
and S the final area when growth is finished. The 
differential equation of the problem takes the form 

s' = ^(S - s) (16) 

which merely implies that growth is proportional to 
the periphery s* and to the available reserves of nutritive 
material. This reserve can be calculated. At the 
beginning the fragment of surface area s includes a 
quantity of reserve material as„. At the end of time t 
the area becomes s and this growth is accompanied by 
a consumption of the reserves proportional to the 
increase of area s — % The amount of reserve e which 
remains available at the time t will then be equal to 

e = «s — /S(s - s ) = (a + /?)s — /3s. 
Hence 

S = s^±i 
Integration of equation (16) gives 

"(Vs + V^)e" tVS - (Vs- VT y 



s = S 



~ Vs ) J 



.(Vs + Vs )fi" (V§ + (Vs 

To simplify this, put 

S = s A a , fiV$ = b = fj,X-Vs Q 

and 



(17) 
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Equation (17) becomes 



S k Sn 



(A-i)- 



,(A + i)fi d + (A-i)J 



This formula conveniently expresses the experimental 
data. In particular it shows that for two tissues 
characterized by the same co-efficient X the final areas 
are proportional to the initial areas. This proportion- 
ality does not exist during growth because of the 
presence of the index e. 



CHAPTER XII 

EMBRYONIC GROWTH 

i. The Free Embryo 

Consider a free embryo. It may be represented 
schematically as a whole containing the organizer, of 
weight w, and its nutritive reserve, of weight n. Sup- 
pose, first, that all the changes within this system and 
between this system and the atmosphere are always 
proportional to these weights; secondly, that the 
absorption of food in unit time is equal to aw -f- bwn ; 
thirdly that the constructive utilization of the food in 
unit time is equal, according to formula (5) of Chapter 
XI, to r(aw + bwn) , if it be unnecessary to take account 
of the set of limitative factors or of intoxication by 
the products of metabolism ; fourthly, that the limita- 
tive term is equal to hw 2 ; and fifthly, that the toxic 
action is cwW, denoting by W the integral 



-f. 



W= I wdt. 




Thus, the equations of the problem assume the form : 

w' — r(aw -f- bwn) — hw 2 - — cwW \ 



n' = — aw — bwn 



(1) 



These equations can easily be integrated. The second 
of equations (1) gives 
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Replace n by this expression in the first of equations 
(i) , expressing w' and w in terms of W, 

W" + AW' 2 = rbW'tn + ~) e~ bw - cW'W 

Multiply this equation by e™ 

(W" + hW*)e m = t(« 6 + a) We (h ~ S)W - cW'W«s* w 

and on integrating we get 

!» = W = F(W) (3) 

which gives 

F(W) = T y + a) Us- ™-e-™) + w a e~™ 
h — o 

c . w c cW , x 
A 2 ^A 2 h V4; 

By integrating equation (3) we get the time t as a 
function of W 

as 

(5) 



* = 



F(s) 



Thus we have the complete solution of the problem 
in parametric form. We shall study this solution. 
Write the equation for the relative growth of the 
organizer 



' VhrinJ) 4- a) c , "I 



br{n i + a) m c 
h-b h 



The theoretical curve has three peculiar points. 
First, the function F(W) vanishes for a certain value 
of W = W w , which corresponds to w = ; secondly, 
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the function w' vanishes for a certain value of W 
= W,„ < W M , which corresponds to the maximum 
value of w — w m ; thirdly the function n vanishes for 
a certain value of W 

W = W = |log(i+^) ( 6 ) 

The fate of the system depends wholly on the food- 
supply. 

First Case. It may happen, that the initial reserve 
of food is insufficient so that 

t(« Z> + a) < w h, w'(o) < 0; 

the weight of the organizer decreases, and the embryo 
dies unless metamorphosis occurs before the disappear- 
ance of the reserves. 
Second Case. The initial reserve is sufficient, so that 

t(w 6 + a) > w h, w' > o ; 

the weight of the embryo begins to increase, reaches 
a maximum and must then tend towards zero. Meta- 
morphosis must be supposed to occur in the neighbour- 
hood of the maximum and to coincide approximately 
with the moment oi disappearance of nutritive reserves. 
The hypothesis of the simultaneous occurrence of 
these two phenomena is expressed by the equation 
W, — W„, which gives the original weight of the 
organizer 

„^(^ + -)( I+ „„i)! 



h-b 

/tr(%& + a) 



+ 



IN 



h-b 
as a function of the reserve of nutritive material n 
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In this argument it has been assumed that the 
utilization of the food is complete. Suppose that the 
reverse is true, and apply formula 15, Chapter XI. 
This gives a slightly different system : 

w' = sw — hw 2 — cwW, n' — — aw — bwn (8) 

The first equation is familiar : it also gives a curve 
with a maximum. Further, equation (6) gives the 
the value of W at the moment of total consumption 
of food. If these two moments are almost coincident, 
metamorphosis occurs at about this critical time. Thus 
the two possibilities, equation (5) and equation (15) 
of Chapter XI, give the same qualitative result. 

2. Embryo in the Mother's Body 

An embryo that develops in the mother's body lives 
in totally different circumstances from those of the 
preceding section. It receives nourishment regularly, 
and the products of its metabolism are more or less 
eliminated. In these exchanges surface area does not 
play so important a part as in the adult condition, 
and it may be supposed that in the absence of limitative 
factors the vital balance of the embryo should become 
positive and proportional to its weight. Food is, in 
fact, proportional to the weight w. Hence, replacing 
X by Xw in equation (14), Chapter XI, we see that the 
utilization of the food is proportional to the weight: 



awth 



® 



Basal metabolism and metabolism of function are 
zero and all this energy is used for growth. Thus at the 
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beginning of embryonic development we return to the 
Malthusian equation 

w' = ew. 

This equation must be completed by limitative terms. 
First, there is a kind of mechanical pressure which 
increases with the growth of the embryo ; further, a 
certain competition between the cells of the embryo 
can be foreseen; lastly, it must be supposed that in 
spite of elimination there are accumulated products of 
metabolism coming from the embryo itself and from 
the increasingly disordered maternal tissues. It is 
extremely difficult to estimate even approximately how 
great these opposing actions will be. A more or less 
accurate qualitative representation can be obtained 
by the equation 

w' = ew — hw 2 — cz^W (9) 

denoting as usual by W the integral 

W = f wdt 

We summarize the peculiarities of the solution of 
equation (9) in the three cases of practical interest 
to us. 

1. In the absence of toxic action, this is the logistic 

law. When w(o) < - the function w increases and 
approac.es i; " 

2. In the absence of internal competition when w(o) 
is very small we have, approximately 

_ 4e%(o)e e( 
W== [2e 2 + w{o)ce el ] z ' (I0) 
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the function w at first increases, reaches a maximum 
equal to 

w - = £ (ii) 

at the moment 

4» = - log — t-\ (12) 

e cw[o) 

and then decreases to zero. 

3. In the general case the solution of equation (9) 
may be represented by the following approximate 
formulae 

w(a)w„ 
w = 




2 

denoting the maximum value of w by w m and the time 
at which w reaches its maximum by t m . When h is 
sufficiently small this time is expressed well enough by 
formula (12). The maximum value w m is equal to 

~c -\- eh— hhv(o)' 



»--'h-tM 



(14) 



Experimental verification is always very difficult, 
for the statistical data are always defective and will 
remain so until a means is found of weighing an 
embryo without killing it. In averages in which each 
individual appears once only the dispersion is too great 
to permit more than qualitative conclusions. We shall 
examine a few numerical examples. 
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Example 1. White Rats. 
following statistics exist : 
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For these animals the 



Ace 
(Days) 


Weight 


(Grams) 


Age 
(Days) 


Weight 


(Grams) 


Observed 


Calculated 


Observed 


Calculated 


13 


■040 


■063 


19 


i-58 


1-64 


14 


■112 


•112 


20 


2'63 


2-54 


15 


■168 


■196 


21 


3.98 


3'68 


16 


■310 


•342 


22 


4-63 


473 


17 


•548 


■59 


Birth 


4-68 




18 


I -000 


I'OO 


— 


— 


— 



Before the thirteenth day the embryos are too small 
to appear usefully in the table, and for the same reason 
the figures for the 13th, 14th, and 15th days must be 
considered doubtful. Some useful conclusions may, 
however, be drawn. At the beginning, limitative and 
toxic actions are slight and w is expressed by the 
Malthusian law 



w = w(o)e ei 



or, in logarithmic form 
log w = 



log w(o) + st 



(15) 



Consequently, with semi-logarithmic co-ordinates, the 
beginning of growth may be represented by the straight 
line (15). If, on the contrary, growth were proportional 
to the surface area of the body we should have, in 
the same conditions : 



w 



ew + 



Wc=L [e(l— k)t -\- W Q X 
or, if w Q is very small 

w ~ [e(l — k)f\~> 



k~- 
3 
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and in logarithmic form, 

log[e(i — A)] log* 

log w a — =-^ — z 1' 



i- ft 



k 



(16) 



hence, with logarithmic co-ordinates log w, log t the 
beginning of growth would be represented by a straight 




Fig. io. — Embryonic Growth of the White Rat 

line. For white rats the possibility of equation (15) 
is realized from the 13th to the 21st day (Fig. 10). 
The conclusion is that in these animals pure embryonic 
growth is proportional to the weight of the body and 
not to its area. The parameter e and the initial value 
of w are equal to 



e = -53; 



w = 7-10 3 gram 



Consider now the final stage of embryonic growth. 
Application of the logistic law shows that the hypothesis 
of the absence of toxic action is inadmissible. In fact 
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the limitative action only becomes noticeable from the 
21st day, but in a most emphatic way, which is in- 
consistent with the logistic law. Equation (9) may be 
improved by supposing that internal competition is 
proportional to a higher power of w. 

w' — ew — hw 1 + "' (m > 1) (17) 



The solution of this equation is 
w = w(o)e ei 



w^he emi + e - w m h 



(18) 



This formula, applied to the case of the white rat, 
will represent its embryonic growth when m is given a 
value of at least 2. Possibly justification may be found 
for putting m > 2, but I do not see it. On the contrary, 
formula (10), which expresses pure toxic action fairly 
closely, represents the observed growth-curve. If 
c = -0295 a satisfactory curve is obtained, reaching a 
maximum at t m = 23-7 days. The calculated values 
of w are found in the third column of the table above. 
Birth must occur a little before the critical point. 

Example 8. Rabbits (after Chaine). For the rabbit 
the following figures exist : 



t (days) 


1+ 


16 


18 


so 


22 


24 


26 


28 


30 


w (obs.) . 
w (calc.) . 


•16 
•17 


•37 
■40 


■94 
■94 


i'6i 

2'17 


5-52 
4-96 


II-O 
II-O 


25-0 
25'4 


40'0 
41-6 


60-0 
59-7 



As in the first example, log w is very well represented 
up to the 26th day by formula (15) with the following 
values for the parameters 

w(o) = 5-10" + gram ; e = -415. 
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Curve (10) reaches its maximum at 



t m = 31-6, 



c == -0014. 



Hence, as in the case of the white rat, birth takes 
place a little before the maximum point. The agree- 
ment between observed and calculated values of w is 
close up to the 18th and after the 22nd days. In the 
interval, iath~2ist, the deviation is considerable. 

Example 3. Guinea Pig (Ibsen and Draper). To 
show the small degree of certainty that statistical data 
on growth present we compare the curves of embryonic 
growth of the guinea pig, which we owe to Ibsen (I) 
and Draper (D) : 



Age (days) 


30 


30 


40 


5° 


60 


64 


67 


"Weight (grams) I . 
Weight (grams) D . 


•062 
•080 


1-32 
1-50 


9'4 

IO-I 


36-3 
29-9 


62-5 
6i-6 


84 


96-6 



The semi-logarithmic straight lines (15) represent 
the observations better than the logarithmic straight 
lines (16), and hence pure embryonic growth is propor- 
tional to the weight of the body. 

For guinea pigs the limitative actions become notice- 
able much earlier than for rabbits or white rats (Fig. 11). 
Also it cannot be affirmed that in this case one or the 
other of these factors predominates. In both series 
the co-efficient of multiplication e is equal to -n. In 
the same way 

h = -ooii, w = -22, c = -00007 

but the calculated curve agrees very imperfectly, with 
the observations, especially for small values of t. 

The other examples more or less confirm our con- 
clusions, which may be stated thus : 
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i. Purely embryonic growth in mammals is propor- 
tional to the weight of the body and not to its area. 
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Fig. ii. — Embryonic Growth of 

the Guinea Pig 

after Ibsen and Draper. 

2. Equation (9) and its two special cases give good 
agreement with the observations. 

3. In some instances residual action seems to be 
more important than the action of internal competition. 

In this respect it may be remarked that the nature 
of residual action is not necessarily chemical and that 
it is quite possible to conceive the accumulation of 
residual effects of a mechanical or physical nature. 



CHAPTER XIII 

POST-EMBRYONIC GROWTH 

i. Human Growth 

In the study of the post-embryonic growth of mammals, 
metabolism and nutrition must he taken into account 
far more than in the study of embryonic development. 
Further, the process must be subdivided into a number 
of stages, each characterized by the constancy of its 
vital co-efficients, with a sharp change in the co- 
efficients on proceeding from one stage to the next. 
We shall consider first the growth of human beings, 
and the better to show the existence of the critical 
points we shall employ the following method, Sup- 
pose that the differential equation of growth takes 
the form 

w' — F{w) (i) 

explicitly independent of the time t. Making use of 
the numerical table of human growth, calculate the 
values of w' and construct the curve F(w) (Fig. 12). 

Consider this curve. During the first eighteen months 
growth is very rapid ; then it becomes slower, doubtless 
because of the change in diet, and remains at almost 
the same value until the ninth year. At this point a 
revival of growth occurs, and increases till the 15th or 
16th year: this is the maximum value of w', later it 
decreases rapidly and disappears at about the thirtieth 
year. After thirty years w' is negative. The form of 
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this curve shows that human growth should be divided 
into several characteristic periods of very different 
duration. 
The first stage, that of breast-feeding, is characterized 

w' 

6- 




20 years 



Fig. 12. — Function w' = F(w) for Man 
after Quetelet. 



by a very high value of w' , and then by a rapid decrease 
to the level of the next period. From the second to 
the ninth year the annual increase in weight maintains 
the same level, with perhaps a slight tendency to 
increase. Suppose that during this time w' is propor- 
tional to some power of w 

w' — ew k (1-5 < t < 9) (2) 
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Human Growth, according to Quetelet 



Age 
(Years) 


Height 

(0 
(Meters) 


Weight 
(Kilos) 


w' 


V 


»«1 


T 





•5° 


3-5 






1-518 


3-04 


I 


•70 


9-4 


3-9 


•145 


2-IIO 


3-ol 


2 


•79 


n-3 


I> 5 


•080 


2-244 


2-84 


3 


.8b 


12-4 


1-45 


■070 


2-315 


2-69 


4 


•93 


14-2 


1-65 


■070 


2-422 


2-60 


5 


I-OO 


157 


i'5 


■°55 


2-504 


2-50 


6 


1-04 


17-2 


17 


•050 


2-581 


2-48 


7 


I-IO 


IQ-I 


175 


•060 


2-673 


2-43 


8 


1-16 


20-7 


175 


•060 


2-746 


2-37 


9 


1-22 


22-6 


1-9 


■055 


2-827 


2-32 


10 


1-2'/ 


24-5 


2'25 


■°55 


2-904 


2-29 


II 


1-33 


2yl 


2-65 


•°55 


3-004 


2-26 


12 


1-38 


29-8 


3.6 


•055 


3-i°3 


2-25 


13 


1-44 


34-3 


4'45 


■055 


3-248 


2-26 


14 


1-49 


387 


4'65 


•°55 


3-383 


2-27 


15 


i-55 


43'6 


5-45 


•055 


3-520 


2-27 


16 


l'59 


49-6 


4.6 


•050 


3-674 


2-31 


17 


1-63 


52'8 


4'i 


•040 


3752 


2-3° 


18 


1-66 


57-8 


2-6 


•°35 


3-866 


2-33 


20 


1-67 


6o-o 


1 '93 


•023 


3-9I5 


2-35 


25 


1-08 


62 -9 


•37 


■008 


3-977 


2-37 


3° 


1-68 


63-66 


•18 


■001 


3-994 


2-38 


40 


i-68 


63-6 


— -OI 


. 


3'99i 


2-38 


5° 


1-67 


D 3"4 


— -09 




3-99 


2'39 


60 


1 '64 


61-9 


— -20 




3-96 


2-41 


70 


1-62 


59-5 


— '21 




3-90 


2-41 


So 


i-6i 


57'8 






3-87 


2-4° 


90 


l-6i 


57'8 






3-87 


2-40 



and that the limitative factors are negligible: it is 
then found, that the index k is positive and does not 
exceed one-third : 



< k < ■ 



(3) 



How can this value of k be explained? It may be 
noted that in this respect man does not differ from 
other mammals. The same phenomenon is found both 
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in rats and sheep. For the explanation we return to 
equation (14) of Chapter XI, which expresses the 
utilization of gross energy X by an animal of weight w. 

x = a.wth I — ■ J (4) 

\a.w / 

A part of the net energy is used by the animal in 
growth. Let x be this constructive energy. Suppose 
that it is related to % by the equation 

x c = xf(w) (5) 

Replacing x by its value from equation (4), we get 

x e — <x.wf(w)th ( — 
\ aw 

Further, w' must be proportional to x e : 

w' = EW k = yx = ycnwf(w)th ( — ) 
This equation enables us to determine the function 

m = - 



o-X 



yath , 

\a.w / 



Two extreme cases may be considered. When 



a.w 



is small, equation (6) gives 



£161* 



and when — is large, we get 

EW k ~ I 
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Suppose that the gross energy X is proportional to 
w m , where m lies between § and i : equation (7) can 
be rewritten as : 

i{w) = 'Ew k - m (7') 

and there are the inequalities 

1 2 

< k < -1 - < m < 1. 

3 3 

Consequently the two equations (7) and (8) can be 
replaced by one 

/(«>) = Hot-" (9) 

in which n satisfies the inequality 

- <n <i. (10) 

3 

Thus, supposing the constructive portion of x to be 
inversely proportional to a power of w varying from 
% to 1 (linear growth — superficial growth — cubical 
growth), we get an explanation of the value of k within 
the limits of expression (3) . 

The following explanation may also be suggested. 
According to the experiments of Ephrussi and Teissier 
on the residual growth of tissues, the hypothesis of 
peripheral growth is in good agreement with the 
observations. It is admitted that this hypothesis may ' 
also be valid for non-residual growth. In this case the 
relation between the constructive energy x and the 
net energy x may be expressed by the equation : 

— — f(w) — Xw~ i (11) 

X 

or, replacing x by its value in equation (4) : 

x = fawHhl — )• 
\ aw J 
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When — is very small this formula becomes 

C/.W 



(- <m < 1 ) 



and when — is large, 

At the age of 9-10, growth, which was beginning to 
slacken, accelerates. If the principal part of the growth 
is still to be represented by equation (2), k must be 
given a value at least equal to unity. Further, the 
form of the curves w and w' clearly shows the action 
of the limitative factors. In fact, the equation 

w' = ew — hw' 2 — cwW, 

with the parameters 

t = 9 yrs., e = -09, h = -000023, c = -000007, 

expresses the observations very well. 

How can the value k = i be physiologically inter- 
preted? Return to equation, putting h ~ 1, and we 
get 

M—jar- 

VK.W — 

\cat) 
in the two extreme cases, we have 

f( w ) = — % I — sma11 ) 

K ' yaX \a.w J 

f(w) = — ( — large 

ya. \a.w 

It appears therefore that after 9-10 years the con- 
structive energy is proportional to the net energy 
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obtained by the organism from its food. Hence at this 
age there is an abrupt passage from one stage to the 
next. This increased growth accompanied by a more 
complete constructive utilization necessarily provokes 
resistances which were absent in the earlier stage. 

It might perhaps be possible to attribute this change 
to the beginning of the activity of the sex glands or 
perhaps to the combined action of these glands and 
others. On the whole it is at the age of 9-10 that the 
change takes place, and at 15-16 that growth reaches 
its maximum intensity. Without this hypothetical 
factor the organism would quickly die. Can there be 
seen in this factor, and in analogous factors appearing 
at critical moments, traces of the evolution of the 
species? It is quite possible that each species has 
passed through critical points and that the introduc- 
tion of a new factor, developed and strengthened by 
selection, has enabled the species to survive. 

A few words may be written on the linear growth of 

the body. The fifth column of the Table shows that 

the increase of stature, V, is very different from the 

function w'. The function I' decreases rapidly from 

t = to t = 4 years, and also after 15 years ; in the 

interval from 4 to 15 it remains almost constant. 

w 
The last column of the Table gives the quotient — . 

v 

If the growth of the body were nomothetic, the ratio 
-j would remain constant. But this ratio first de- 
creases, passes through a minimum about the 12th year, 
and then increases, which means that until the 13th 
year the increase in stature is more rapid and after the 
12th year less rapid than that of weight. 
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This confirms the suggestion that a profound differ- 
ence exists between increase in weight and increase in 
height, and explains the differences observed between 
the basal metabolism and the area of the body calcu- 
lated from the formula 

Surface ~ [mfi. 

2. Mammals 
Problems of growth are not the same for warm- 
blooded and for cold-blooded animals. In the former, 
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Pig, 13.— Function m' = F(w) for the Rat 
after Donaldson. 



as in man, there may be observed critical moments in 
the development, angular point's in the growth-curve ; 
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and the hope of expressing the growth of an animal 
by a simple formula has to be abandoned. This dis- 
continuity of the growth-curve was first established 
by Lotka, and we owe the detailed study of the 
phenomenon to J. S. Huxley and G. Teissier. 

Rats. Lotka, who treated the problem of growth as 
a statistician, noted that the first three months of a 
rat's life could be expressed by a logistic curve, but 
that later this curve gave values much lower than 
those observed in practice. As a matter of fact the 
first three months of the life of a rat are not very closely 
represented by the logistic curve. To confirm this, 
look at the empirical curve (i) for the rat (Fig. 13). 
Just as in man, it shows critical points. The first of 
these occurs at the age of twenty days and separates 
a phase of growth in which w' is almost constant from 
the following phase in which w' rapidly increases. This 
lasts until the 75th day, after which w' decreases fairly 
rapidly, almost in a straight line, and reaches zero on 
about the 290th day. Curve (1) shows clearly that the 
development of the rat consists of at least three different 
phases. It is curious that the differential equation 

w' = £ze>* — hw f ' (12) 

leads to a growth law of the type 

w = A[thX{t + t)]« (13) 

which well represents the weight, w, of the rat during 
the whole of its growth, but curve (12) does not corre- 
spond at all to that in Fig. 13. 

Guinea Pig, according to Minot. The growth curve 
of the guinea pig presented by Minot is a very interest- 
ing one. The curve (Fig. 14) passes through two 
maxima separated by a minimum. The first maximum 
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occurs on about the 15th day, the minimum about 
the 100th day, and the second maximum about the 
180th. It looks as if the growth energy, weakened after 
the first maximum, experienced about the 100th day 

w'| 
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Fig. 14. — Function w' = ~F(w) for the Guinea Pig 
after Minot. 

a recrudescence under the influence of a factor which 
had remained hidden until this moment. The pheno- 
menon might be compared to a rejuvenescence of the 
organism. 

The growth of the sheep is remarkably well repre- 
sented by the equations 

w' — a — bw, w = tei~ (w — w )e~ ai , (14) 

and this is true for several breeds. 

In large cattle there is a critical point with the 
transfer from a milk diet to pasture. Similarly there 
are seasonal critical points corresponding to changes 
from summer to winter regime and vice versa. It 
seems that these changes not only affect the co-efficients, 
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but the exponents as well. Hence these seasonal 
variations must always be taken into account in the 
study of yearly statistics. 



3. Fish 
The growth of fishes also presents some interesting 
peculiarities. We have chosen an example which shows 
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Fig. 15. — Growth of thb Bream, after WasnetzoS 



A. Southorn Bream : 

B. Northern Bream : 



. Volga- — Caspian. 2. Azov. 
:. L. Xuusula. 3. L. Pslrov. 



the importance of the geographical factor in splitting 
a species into varieties. Wasnetzoff gives very com- 
plete statistics of the growth of the bream in the 
fresh waters of the European parts of the Soviet Union. 
This vast region contains the basins of rivers flowing 
into four seas, and may be roughly subdivided into 
two parts, one of which feeds the Baltic and the Arctic 
Ocean, the other the Black Sea and the Caspian Sea. It 
seems that the growth of the bream in the northern half 
is governed by different laws from those in the south. 
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Fig. 15, A and B, gives a comparison between the 
growth-curves of the bream in these two regions. The 
abscissae represent the logarithm of the age and the 
ordinates the logarithm of the length of the body. It 
is seen that : (1) the original length of the northern 
bream is less than that of the southern ; (2) the growth 
curve of the northern bream is almost rectilinear, 
while that of the southern curves towards the abscissa. 
The difference in original length may be explained by 
the action of climatic factors. The delayed action of 
the limitative factors in southern bream is more difficult 
to explain. Many hypotheses may be suggested of 
which the most plausible is that we are witnessing 
the appearance, perhaps already completed, of two 
varieties under the influence of different climatic factors. 

The growth of the bream, in its initial stage, may be 
expressed by the formula 

I = I/. 

The average values of l x and i for the two varieties 
are: 



Bream 


h 


i 


Northern 
Southern 


5'i 
8-i 


■87 
i'o6 



The values of * do not differ much from unity. The 
whole curve for the southern variety can be calculated 
from the formula 

hi 

I — L— (L — l )a 
which corresponds to the differential equations 



(15) 



3 



i), 



IS) a= BW* — 



(16) 
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The considerable deviation of the index i about the 
average i = i may be explained by the fact that the 
length is not proportional to the cube root of the weight. 

4. Growth and Ecdysis 

In many animals the relation between moulting and 
growing takes an interesting mathematical form. In 
arthropods linear growth between the moults is pre- 
vented by the rigid chitinous skeleton, while increase 
in weight continues to occur. At ecdysis the animal 
rids itself of its tight skin, another more suited to its 
weight replaces it, and the animal soon adapts itself 
to the new regime. Przibram and Megusar supposed 
that the weight of an arthropod should be doubled 
between one moult and the next, which corresponds 
to a doubling of every cell in the body. In consequence 
linear dimensions should increase from one moult to 
the next in the proportion of V2 = i"26. 

In spite of the existence of many aberrant cases this 
rule has been confirmed by observation. In some 
aberrant cases the linear dimensions are multiplied 
after a moult by 1-59 = i-26 8 or by 2 = I-26 3 . Thus 
Przibram's rule corresponds to an objective reality, 
but there are many factors which may make growth 
diverge from it. On the one hand, growth does not 
always take place in a nomothetic way : on the other, 
the body of an animal has a complex internal and 
external structure, and the organs and parts of the 
body do not grow proportionately to its weight. 

Finally, indefinite growth is impossible, for limiting 
factors intervene in the invertebrates just as in the 
vertebrates. These factors may act in two ways, either 
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by increasing the interval of time between two moults 
or by reducing the growth between one moult and the 
next. Lengthening of the 
time-interval is seen in 
several growth - curves, 
for example in Dixippus 
morosus and Tenebrio 
molitor. The curve for in- 
crease of weight between 
two moults in Tenebrio 
molitor resembles the 
growth-curves in the case 
of residual actions. The 
weight at first increases 
rapidly, then the increase 
is retarded, the weight 
reaches a maximum 
some time before the next moult, and finally decreases, 
but the moult intervenes and the next cycle begins. 

It would be very interesting to determine the nature 
of the factor that produces this decrease of weight 
before a moult. Dehydration of the tissues may be 
suggested, as well as chemical action. It is a factor 
which may produce actions of the "hereditary" type, 
in the mathematical sense of the word, as the experi- 
ments of Ilyin have shown. 



150 days 

Fig. i6, — Individual Growth of 
Tknbbrio molitor, after Teissier. 

A, weight in function of time. B, spread of 

tbc antennae in function of time, m, m, 

minima corresponding to the moulds. M 

maximum shortly preceding a moult. 



5. Growth of Organs and Parts of the Body 

Denote by q the weight of an organ in an animal of 
total weight w. There is an empirical law connecting 
these two quantities 

q = Kw* (17) 
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or, in logarithmic form 

log q = log K -f t log w. (18) 

This law is based upon a great number of observa- 
tions, and within fairly wide limits of variation of 
w and q the growth-curve of q in relation to w is ex- 
pressed in logarithmic co-ordinates by a straight line 
or by a broken line composed of straight elements. A 
question may be put here. What is the meaning of 
equation (17) and what are the values of i? If the 
growth of the organs were proportional to the weight 
of the body, then i would be equal to unity. 

But i is almost always different from 1, and truth 
to tell there is no reason why it should equal 1. In 
many cases organs may be considered to be superficial 
formations. It is then not a matter of the outer surface 
of the body but the surface of the embryonic tissues, 
the foldings, invaginations and proliferations of which 
produce the development and structure of the body. 
It may be supposed that a superficial excrescence is 
nourished through the portion of the surface that it 
occupies, and if this portion is not too large the weight 
of the excrescence will be proportionate to this area 
and not to the weight of the body. I do not mean that 
things always occur in this way, but the hypothesis 
that i is different from 1 is in every case much more 
probable than that i should be equal to 1. G. Teissier 
calls laws of type (17) when the exponent i is different 
from 1, laiys of disharmony. It would appear rather 
that there is not harmony when i ~ 1, and that there 
is not disharmony when i is different from 1. 

To confirm this point of view, two numerical exam- 
ples may be considered. In the particular case of the 
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rat's heart, i = 75, a comparison with another em- 
pirical law affords an explanation of this value of i. 
The length of the rat's body also follows a law of 
type (17) 

I = ~Lw' 3 (19) 

and the surface area of its body, s, can be represented 
in the form 

s = Sro' 7 (20) 

Hence the weight of the rat's heart is proportional 
to the area of its body. 

The other case concerns the orthopteron Dixipfius 
morosus, of which G. Teissier has measured the total 
length, the length of the abdomen, the length of the 
prothorax, the length of the head and the diameter of 
the eye at different periods of its life. Let w and I be 
the weight and length of the body, w ± and l x the weight 
and length of the head, and w 2 and l 2 the weight and 
diameter of the eye. Suppose that 

I = Lw*, ly = ~L l w l i , l 2 — LjjWa*- 

Suppose, now, that the head is a superficial excres- 
cence of the body and the eye is a superficial excrescence 
of the head. In such circumstances we get 

Wx = wy - Q X P, w z = w 2 v = Q/ 

and in consequence 

l x = ijfi, i z = jy». (21) 

Suppose that the length of the body increases pro- 
portionately to a 

I' = al. 
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Equations (21) show that the corresponding values 
of ^ and l 2 become 

and hence 



V 

7 BSfl ' 



l'x 



h. 



(22) 



/ 

t, "a 

The following table gives the values observed by 
Teissier and those calculated from formulae (22) : 





15 


19 


28 


37 


A 6 


64 


85 


'1 


1-15 


i-35 


1-65 


2-1 


2-6 


3"i 


3-8 


'2 


•38 


•43 


•47 


•57 


-65 


•75 


•85 


V 

a = J 


1 


1-27 


1-87 


2-47 


3-07 


4-27 


5-C7 


a § = ~ obs, 


1 


1-09 


i-43 


1-83 


2-26 


270 


3-30 


ji calc. 


1 


1-17 


I'52 


1-83 


2'12 


2-62 


3-18 


at = -p obs. 


1 


1-13 


1-24 


1-50 


I-7I 


i-97 


2-24 


-p calc. 


1 


I'll 


1-32 


1-49 


1-64 


1-90 


2-l6 



If the uncertainty to which these measurements are 
subject be taken into account, it cannot be denied 
that the agreement is remarkable. The hypotheses 
used are very approximate and need many important 
corrections, but in this particular case they completely 
explain the observed phenomena. 



CHAPTER XIV 
FORMS OF LIVING OBJECTS 
In this chapter we confine ourselves to a few remarks 
on the role of symmetry and polarity in the economy 
of an organism, and on the role of movement in the 
exterior shape of living objects. 

i. Symmetry; Polarity 

In most cases the forms of living objects are sym- 
metrical either about a point, an axis, or a plane. 
Further, there are often poles or points of activity, 
either internal or peripheral, in situations that are 
related to the type of symmetry of the organism. 
This phenomenon may be explained by the principle 
of economy of the organism, a maximum result at a 
minimum cost, and by considerations of stability. 
From the mathematical point of view these two prin- 
ciples are united in the extreme problems, in which 
are sought the conditions causing an expression to 
assume an extreme value, either maximum or minimum. 
It is not possible to assert that all the forms of living 
animals are very satisfactory in this respect, but a 
general tendency is undeniable. 

Sphere. Begin with the spherical form. Of all bodies 
of the same surface area the sphere has the largest 
volume, and of all bodies of the same volume the 
sphere has the smallest surface. These two extreme 
principles are of importance in the exchanges between 
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a living object and its environment. When these 
geometrical considerations are supplemented by a study 
of physical and chemical forces, the sphere becomes 
the figure of equilibrium. Thus a convex mass of 
liquid, either homogeneous or heterogeneous and strati- 
fied, assumes the form of a sphere under the action of 
hydrostatic pressure, osmotic pressure, capillary attrac- 
tion, and viscosity, and it is the only form of equilibrium 
possible in the absence of other forces or of internal 
movement. If an irregular mass is divided into several 
parts it is very probable that these will take on a 
spheroidal form. If in a spheroidal mass the osmotic 
equilibrium with the surroundings is disturbed, the 
mass changes its volume but retains its spheroidal 
form. 

Spheroidal arid Ellipsoidal Forms. Radial symmetry 
can be preserved when the living mass contains a focus 
of activity nearly coincident with its geometrical centre, 
but this is relatively uncommon ; hence the appearance 
of a focus of activity produces a change from the 
spheroidal to the elongated or flattened ellipsoidal 
forms, The intervention of the force of gravity produces 
anon-concentric distribution of the materials composing 
a cell. An internal complexity due to one of these 
causes or to their combined actions changes a spheroidal 
body by creating an axial in place of a radial symmetry. 
A specialization of the surface accompanies this, with 
permeability, capillary tension and osmotic pressure 
varying from one pole to the other. 

Movement, Suppose that an organism thus consti- 
tuted comes to propel itself by ciliary action, pulsations, 
or any other method. In the case of vertical motion, 
axial symmetry may be retained ; in horizontal motion 
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a deformation is produced, which increases as the 
normal speed of horizontal motion increases ; and the 
body loses its axial symmetry, replacing it by symmetry 
about a plane. This is the commonest case. 

Divisions. Consider now the disturbances of equili- 
brium utilized by organisms. The approach of other 
bodies may produce, directly or by induction, local 
changes of superficial permeability and the appearance 
of pscudopodia. In the same way the disappearance 
of the pseudopodia when the superficial stimulus ceases 
to act may be explained. It is curious that this pheno- 
menon which appears to be so disadvantageous to the 
organism, should come at length to be used so profitably. 

Cell-division is another form of equilibrium disturb- 
ance. In another connection Henri Poincare, discussing 
the formation of the heavenly bodies, studied the 
problem of equilibrium-forms of a homogeneous fluid 
mass rotating under the Newtonian attraction of its 
particles. This problem, although very different from 
the one we are at present considering, nevertheless 
offers very useful analogies and suggestions. It is 
known that the spheroidal form of a heavenly body is 
stable; further it is known that the Earth-Moon 
system, composed of two spheroidal bodies, is stable. 
It must, however, be admitted that this system was 
formed by the break-up of a single gaseous mass, 
owing to the loss of stability. Consequently in com- 
paring the present state of the Earth-Moon system 
with its hypothetical distant past, a single gaseous 
mass, we find a series of stable forms at the beginning 
and at the end, but losing their stability at some 
unknown intermediate moment. Similarly, a cell or 
a protozoon is relatively stable before division and the 
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two daughter cells become stable some time after 
division, but these two ends of a chain are separated 
by a link in which stability is disturbed. Henri Poin- 
care\ discussing the series of equilibrium forms, has 
demonstrated the existence of critical points at which 
the series ceases to be stable, and stability then charac- 
terizes a new series of equilibrium forms. In this way 
the division of the gaseous mass into two separate 
masses is foreseen. 

A similar phenomenon probably occurs with cells. 
The spheroidal form may be stable, while the forces 
in equilibrium satisfy certain inequalities, for example 
as long as the mass of the cell does not exceed a certain 
limit, depending on all the factors present. Near this 
limit stability vanishes and the smallest external action 
or the least internal asymmetry may produce a change 
of form or the division of the mass into two or more 
portions. 

Division of a Fertile Egg. In the particular case of 
an egg there are four interesting circumstances : (i) an 
unfertilized egg possesses a really remarkable stability ; 
(2) the first breakdown of stability is due to a foreign 
body penetrating as far as the nucleus, and the plane 
of the first division passes through the two poles of the 
egg and the path of the spermatozoon ; (3) the divisions 
that follow have no need of external stimulus ; (4) the 
first two divisions still keep a certain symmetry and 
produce four blastomeres to some extent similar to the 
original egg, while later divisions produce a deeper and 
deeper differentiation. 

This last peculiarity is explained by the plane of the 
first division passing through the axis of symmetry. 
The four blastomeres resulting from the second division 
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have an elongated shape and are least stable in the 
equatorial region, which explains the orientation of 
the third division. Further the ahomothetic stratifica- 
tion of an egg explains the difference between the four 
micromeres and the four macromeres resulting from 
the third division. We shall not proceed to the later 
stages of the formation of a multicellular organism, 
although fairly simple geometrical and physical con- 
siderations seem able to explain some phases of the 
process. 

What is less clear is the automatic mechanism which 
performs all these divisions after the first loss of 
stability. These uninterrupted divisions cannot be 
attributed to the presence of a foreign body, for 
parthenogenetic development obeys the same laws. 
Everything proceeds as if an essentially unstable 
system were held, pending the action of a releasing 
agent. This is very different from the division of 
unicellular organisms, producing in general a series of 
identical generations, separated by periods of stability. 

2. Environment, Movement, and Form 

The shape of a living organism depends in general 
on many and varied factors, but in certain cases the 
action of a small group of these factors dominates and 
wholly explains the external form. This is true of 
fishes, which lead an active mobile life in water, and 
sometimes out of it. By the phenomenon of convergence 
similar forms are found in other animals, very different 
from fish in origin and internal structure, but living in 
the same environment. 

Consider first several experimental and theoretical 
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results from aero- and hydrodynamics. A body moving 
through a fluid needs a motor, a propeller, a method 
of changing speed and direction, and a stabilizer. The 
technical problem may have several more or less 
economical solutions. 

First, horizontal motion must be distinguished from 
vertical motion. In horizontal movement in a liquid 
medium, without being supported on the bottom, there 
is a necessity for a propeller driven by a motor to 
produce, maintain, and modify the speed, and of a 
suspensor to preserve a certain level. Consider pro- 
pulsion first. A propeller is opposed by the resistance 
of the medium and by friction. Resistance is propor- 
tional to the square of the velocity and to the density 
of the medium, 

R = k P v 2 
where R is the resistance, p the density, v the velocity, 
and A is a constant depending on the dimensions and 
shape of the body. Aerodynamic experiments show 
the sometimes enormous advantage of fish-shaped 
bodies over spheres, cylinders, and other forms. A 
fish-shaped body can move in a fluid with less effort 
than a much smaller sphere. The reason may perhaps 
be attributed to the formation in the surrounding 
fluid of discontinuities, separating regions of different 
conditions. Friction between the body and the medium 
is often negligible, but viscosity between two fluid 
regions of different conditions contributes to the 
formation of rotary movement, of whirlpools behind, 
before, and sometimes at the sides of the moving body, 
and these whirlpools increase the resistance to a 
remarkable extent. Sometimes they accompany the 
solid body as it moves when it is as if the propeller 
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were forced to drive the body itself plus the disturbed 
regions behind and before it; sometimes they detach 
themselves and their sucking action produces an 
increased resistance. But the whole system formed by 
these regions and the body itself presents a kind of 
fish-shaped figure. Suppose this system were solidified. 
The current admirably adapts itself to this fish-shaped 
model and the discontinuities due to the movements 
of such a body are overcome without the formation of 
dead water, whirlpools, etc. Experiment shows that 
for a given velocity there is an optimum relation be- 
tween the dimensions of a fish-shaped body. If the 
body is too long, whirlpools appear at certain critical 
places on the lateral surface, if it is too short, whirlpools 
are formed in front and behind. Fred6ric Houssay, to 
whom belongs the merit of having perceived this 
relation between the form and the movement of fish, 
also foresaw the existence of these critical spots, which 
coincide with the fins ; he is mistaken only in the nature 
and exact role of the whirlpools, and hence his hydro- 
dynamic analogies appear to be doubtful. 

We pass on to vertical movement. The first thing 
to be done is reduce the difference between the weight 
of the body and the weight of liquid occupying the 
same volume. This little technical problem is very 
easy to solve in the case of water and the density of 
fishes differs very slightly from that of water. With air 
it is more difficult, and the specific gravity of birds and 
insects, though much reduced by a special arrangement 
of the tissues, still remains much higher than that of air. 

Aeronautical technique admits four methods of 
vertical locomotion : 

I. A special vertical propeller, as in helicopters. 
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2. Difference of specific gravity, as in balloons. 

3. A horizontal propeller in conjunction with inclined 
planes. 

4. Vertical currents in the atmosphere. 

All these methods are used by living objects, and 
there are others which our technique has not succeeded 
in reproducing. The swim-bladder serves to change 
the specific gravity of fish, and is not intended for 
stabilization. The fins also seem to play the part of 
inclined planes. In some cases the whole body of the 
fish is transformed into an inclined wing, as in certain 
Trygonids, and these fish closely resemble aeroplanes 
or hydroplanes. This adaptation of form to medium 
makes all other modes of vertical displacement un- 
necessary, and very well illustrates the principle of 
economy. 

There exists a fish able to propel itself in two media 
as different as air and water — the flying fish, whose 
extremely aerodynamic shape, with its well-developed 
pectoral fins, is altogether comparable to that of an 
aeroplane. In water the caudal fin acts as a propeller, 
with intermittent action, like the scull at the stern of 
a boat. This propeller gives the initial soaring move- 
ment which allows the fish to leave the water and 
accomplish a fairly long volplane. According to the 
experiments of Schouleikine with a model of this fish 
placed in an air current in a wind tunnel, the support 
given by these fin-wings is quite enough to hold 
the fish up during its flight. It appears furthermore 
that the flying fish can make use of ascending air- 
currents. 

We shall not linger over the form of birds, as this 
is a slightly more complicated example of the same 
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principle with a different propeller. It may merely be 
remarked that birds make great use of ascending air- 
currents, which enables them to reduce their expendi- 
ture of energy to a minimum. 



CHAPTER XV 

EVOLUTION 

i. Life and Entropy 

Before attacking the problem of evolution, it is 
necessary to understand the exact meaning of the word. 
All existing definitions vary between two extremes. 
Some apply the word to any natural process whatever, 
contrasting evolution with immobility; others look 
for an expression of progress, describing as evolution 
all that is progressive, ascendant, and opposed to 
levelling. The former try to preserve an objective 
attitude, the latter paint the universe rose-coloured or 
black, according to their inclinations. 

It seems to me that the word may be given a quanti- 
tative definition by applying it to every process that 
leads from a probable state to another less probable, 
or which is opposed to the increase of entropy. The 
objection to this definition may be that in applying it 
to partial processes or to ascending portions of cyclic 
processes, there is a risk of losing oneself in the intan- 
gible. But it is not necessary to do so. The distinction 
is very difficult to determine, as difficult as that in 
the paradox of the heap of corn. However, the heap 
exists and is very different from a grain. Evolution 
exists, and is very different from a particular, even an 
ascending, process. 

Evolution local, Entropy universal. At this point a 
very important distinction must be established — 
between local tendencies and universal tendencies. 
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The latter will, without any doubt, move towards a 
more probable state, if this state is not already reached. 
Local tendencies may be ascendant or descendant and 
there is no doubt that the probability of the former is 
far less than that of the latter, but it is not zero. The 
realization of a slightly probable state is not impossible, 
and the non-realization of such a state is very slightly 
probable. Hence evolution must be a very rare pheno- 
menon, highly localized both in time and space. Com- 
binations that are very rare but which are more or less 
regularly realized may form a kind of separate universe, 
superimposed on chaos. If the appearance of matter 
is due to these chance combinations, each particle is 
a local phenomenon, but their aggregate may have a 
universal character and last a sufficiently long time to 
permit the realization of events much less probable 
but not impossible, such as living matter. Consequently 
life is a link in an ascending chain of very slightly 
probable events, the realization of which is in strict 
agreement with the calculus of probabilities. This 
ascent has not come to an end. Its latest stage is that 
which leads from the appearance of life to the appear- 
ance of human consciousness. 

As a matter of fact this forward march is accom- 
panied by a great wastage, and adaptation is inevitably 
accompanied by all the imperfections of Nature; all 
the species that have not succeeded, all the germs that 
have perished before maturation, before fertilization, 
before hatching. Each living object is in some way a 
winner in the lottery of life. Considered from this 
point of view the vital process ceases to be an exception 
to the law of entropy and becomes its most striking 
confirmation. 
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2. The Struggle for Existence. Selection 

In the process of evolution, which has effected the 
transformation of species, what is the exact role of the 
struggle for existence? In the preceding chapters we 
have examined several cases in which a slight dis- 
advantage has led to the total, disappearance of one or 
other of the groups composing the population of a 
neighbourhood. We have seen that the disadvantage 
may be of a varied nature, but it always amounted to 
additions to the mortality. Hence the mechanism of 
elimination of the least adapted becomes quite clear 
when it concerns animals belonging to already existing 
species, but how can it apply to related groups belong- 
ing to the same species and giving fertile hybrids ? 

An example taken from Gause showed us a popula- 
tion composed of two different races of Didinium 
nasutum feeding on Paramecium. All was well in this 
little world as long as food was sufficient, but a scarcity 
of Paramecium led to a rapid disappearance of the 
giant strain of the predators, while the dwarf strain 
was able to survive for some time. What happens to 
the survivors of the dwarf strain on the reappearance 
of food? The population of Didinium is restored. By 
a mechanism which escapes us, but which is perhaps 
related to the behaviour of the chromosomes in con- 
jugation, a giant mutation may reappear and multiply. 
It is to be supposed that their conjugations with the 
dwarf strain produce hybrids and that Mendel's laws 
determine the size of these hybrids and of their pure 
or mixed descendants. 

In what follows we are not concerned with the 
mechanism of the appearance of mutations, We shall 
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consider a population which is composed in the first 
place of one hybrid and two pure groups. Each of these 
groups has its own vital co-efficients. In the case of one 
homogeneous isolated species the logistic equation 

p 1 = np — mp ~ vp% — iip" 1 (i) 

gives the increase in numbers of the species. The co- 
efficients of mortality m and p do not require special 
treatment. Each group has its own. The co-efficients 
of fertility n and v may vary in the different crossings 
according to which of the three groups includes the 
mating pair. In the simplest case there are nine types 
of crossing if the species is bisexual, and each type may 
have its particular fertility. Let p x and p z be the 
numbers of the two pure strains and p 3 that of the 
hybrids. If the crosses between individuals of these 
three groups are not controlled by any special law the 
number of births resulting from crossings between 
groups p { and p k must, without the intervention of 
limiting factors, be equal to n^pip-klp where 

P=Pl+p2+P 3 - (2) 

Suppose that the action of these factors on the females 
decreases the fertility n ik so that 

°* = v vd>i + h k pz + v slc p 3 

and the analogous action on the males has an effect 

r i = PilPt + PiSpZ + Pi3p3- 

In these circumstances the corrected co-efficient of 
fertility is 

N« = »« ~ er ft - r«. (3) 

Further, the offspring of these crosses is distributed 
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according to Mendelian laws. We finally get the 
differential equations of the problem 

, , NuA* , (N 13 + N 31 )A^ 3 , Nrf 

* 1 = — + Tp + -# — m ^ 

~ Mftih + /W>2 + fhzPa) 
., N 2 ^ 2 (N 23 + N 32 )A &, Nrf 
*■ = — + — ^ + ^- 

- m2p s — puitenpi + Unfa + /W^) 
(N la + N 21 )AA , (N 13 + N 31 )^3 



^3 = T h 



Zp 



(4) 



"^ 2^> 2£ 

— OT3A — PaiHiaPi + HiPz + MaaPa) ■ 

The problem of selection appears in this form when 
it is concerned with two strains giving fertile hybrids. 
If other species inhabiting the same station are in- 
volved, system (4) must be completed by new terms 
and equations. 

Consideration of equations (4) shows that : 

1. In conditions unfavourable to the hybrids the 
more numerous pure strain will survive ; consequently 
a new mutation which must necessarily be restricted 
in numbers has a very small chance of survival. 

2. In conditions that are very favourable to the 
hybrids the three groups can co-exist and tend towards 
a stable equilibrium. 

3. In the case in which the hybrids occupy an inter- 
mediate position between the two pure strains, only 
the pure group with the more favourable vital co- 
efficients survives, 
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We proceed to confirm these conclusions in a few 
relatively simple special cases. 

Special Case 1. Suppose that the coefficients v ik 
and p iS are zero and the only difference between the 
three groups consists of different mortalities. Equa- 
tions (4) are simplified and become : 

Pi = j »fA - h PiP 

P 3 « — J ; »hp a - HzP 

In the case of a low birth-rate less than the death- 
rate : 

n < m v n < m 2 , n < m 3 (6) 

all three groups tend towards zero : 

lim fix — lim fi z = lim p a — o (7) 

( — > do t -»- DO I ->• DO 

If the death-rate of the first pure strain is low 

n„> m v m 3 > m t (8) 

the hybrids and the second pure strain vanish, and we 
get 

hm p t — — r— *, lim & = o, hm ^ 3 = 0. (9) 

t~>DO " < -> DO < -» DO 

If the death-rate of the second pure strain is low 

n >m 3 , m 3 >m 2 (10) 
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the hybrids and the first pure strain disappear : 

lim j> x — o, hm fi 2 = — 7 — -, hm p 3 = o. (n) 

* — > DO * -*• DO " t — > DO 

This implies that the survivor of the pure strains is 
the one whose co-efficient of increase is positive and 
whose death-rate is lower than that of the hybrids. 
What happens if the inequalities (8) and (10) are both 
realized? In this case, as has already been said, 
everything depends on the original values of ^x(o), 
ft 2 {o), and ft 3 (o), and the less numerous of the pure 
strains vanishes with all the hybrids. Finally, in cir- 
cumstances that are very favourable to the hybrids, 

w 3 < m lt m z < m z (12) 

co-existence of the three groups is possible. 

When the value of m 3 lies between % and m 2 , the 
survivor will be the pure strain whose death-rate is 
lowest and whose co-efficient of increase is positive. 

Special Case 2. Take another case in which the 
co-efficients v (k and p ik are zero, the death-rate and the 
limiting factors are the same for all groups, but fertility 
varies according to the crossing. To fix our ideas put 

W,t 4" tlS / 

«U = n i> n 2Z = «B» »« = , (13) 

2 

that is to say, fertility in the hybrids is intermediate 
between those in the pure strains. Ih this way we get 

% + n % 



Was = «xa = n i- 

3% + n» 

n 13 = « 31 s= 

4 

4 



(14) 
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Equations (4) then become 

jp 1 = ~r- i n — — + 
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■ w/> x — hp x p 



8p 



h' ^*^t* I (3 M 2 + w x )^a . K + ^ 2 )A 2 
^ 2 p~. 4p 8p 



P'z 



— mp 2 — hp & p 

(3% + Wa)A& , (3«2 + %)jV>3 

(j% + « 2 )^ 2 (Mi + ra„);p 3 : 



(15) 



- mp s ~ hp a p. 

The three states of equilibrium are in this case incon- 
sistent with one another. In the first case, in which 
fertility is low 

n t <m, n 2 <m (16) 

all groups vanish. 

In the second case the co-efficient of increase of the 
first pure strain is positive and its fertility is greater 
than that of the second pure strain: 

n t >m, % > n a (17) 

which in the limit leads to the survival of the first pure 
strain: 

lim p x — -^-=r — , Una p 2 = 0, Um p a = 0. (18) 

< -» do " 1 -> DO f -> DO 

The third case is the opposite of the second: 

n a > m, w 2 > % \ 

i- , nn—m .. L \ (19) 

lim & = 0, hm^ a = —; — , hm p 3 « 0. 

(-v-30 '->DO "■ <-»DO } 
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Thus in this case the most fertile strain alone sur- 
vives, provided its co-efficient of increase is positive. 

Special Case 8. Suppose that the fertility and 
mortality are the same for all groups and all crossings, 
and that the limiting co-efficients satisfy the following 
relations : 

Pi* = v nc = O) Put — Pv lhn — Pa> fhio — Pa 
so that equations (4) become 

'3 



P 



p\ = m p t - H p x p 



^ 3 = ■- mp 3 - jj, 3 p z f 

Everything depends on the limiting co-efficient ^3 
of the hybrid group. Suppose that n — m > and that 
,« 3 is intermediate between p x and p % . In this case the 
most resistant pure strain will survive, that which has 
the smallest limiting co-efficient. 

Suppose again that 

Pa > Pz> fa > Pv n— m > o. 
In this case the less numerous group disappears very 
quickly and so does the hybrid group. 

Finally, in circumstances that are very favourable 
to the hybrids, with 

Pa < Pv P% < P%> n~m>o 
the three groups can co-exist. 
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Thus the final result of this discussion is in agreement 
with our conclusions. 



3. Indirect Selection 

Suppose that a species p composed of two pure 
strains p x and p 2 and their hybrids p s occupies a station 
in conditions favourable to the hybrids. Equilibrium 
is established and maintained until the appearance of 
another species x, which upsets the equilibrium in 
this little world. The mathematical and biological 
problem which results from this is too complex to be 
attacked here in general terms, and so we confine 
ourselves to the study of a particular case. Suppose, 
to fix our ideas, that the new species x serves as food 
for the species p ; the equations of the system have the 
following form : 



P\- 



™[Pi 



+ 



P s 



mp t — hp x p + l x p x % 



P\ = 



«(A + ty 



mp % — hp s p + l 2 p 2 x 



P'* 



Pi 



MPi+ r i){p* + i; 



p% 



up, 



x' = E*~-Hx s 
also M 



m, 



— hp a p + l 3 p a x 
WhPi + hP% + hPa) 
and e — 



(21) 



■m > 



(22) 
Suppose first that the fertility n is fairly large and 
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the co-efficients of utilization of the food satisfy the 
following inequalities : 

l x > l a , HE > eU x 

These conditions are favourable to the first pure strain, 
and result in the ultimate disappearance of the second 
and of the hybrids. The first pure strain survives in a 
state of equilibrium with the food-species 

E^ + eH 



lira p x = 
lim p & 
lim x 



u? + m 

o, lim p 9 

hE - eMj_ 



tif + Hh 
If the second pure strain is favoured 
l z > l a , hE > eU 2 



when 



lim p x 



lim^ a 


Eh 


+ ftH 


u^ 


+ Rh 


lim A 


= 




lim x 


hE 


— eM% 



us + m 

In the case of insufficient multiplication, 

WE 
S< ~~R' (^ = I, 2, 3) 

the predatory species vanishes : 



(23) 



(24) 



(25) 



E 



lim p x = lim p^ = lim p 3 — o, lim x — — (26) 
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Finally, if l s is large enough and e is positive, the 
co-existence of the three predatory groups and their 
prey is not impossible. Thus the ultimate fate of the 
system depends chiefly on the hybrid group and espec- 
ially on the l 3 , the co-efficient of utilization of food. 
If l s is intermediate between l x and l % , the survivor is 
the pure strain that makes greatest use of its food. 
If the hybrids make less efficient use of the food than 
either of the pure strains, the result depends on the 
initial composition of the predatory species. If the 
hybrids make better use of the food the co-existence of 
the three groups and their prey is possible. It goes 
without saying that supplementary conditions expressed 
in inequalities are in every case essential in order that 
these conclusions may be valid. 

All these phenomena remain possibilities if e is nega- 
tive, that is to say when species x constitutes almost 
the sole diet of species p, with the proviso, however, 
that inequalities (25) are not satisfied. 

No doubts can therefore be entertained as to the 
reality of the working of selection. In the case of two 
closely allied species which do not produce hybrids, 
the probability of the disappearance of the less favoured 
species is greater than that of the co-existence of the 
two. In the case of two strains which produce hybrids, 
the ultimate fate of these strains depends on the vital 
co-efficients of the hybrids. If the hybrids are less- 
favoured than the pure strains, it is most probable 
that the less numerous strain will disappear, and this 
fact determines the fate of all mutations that produce 
feeble hybrids. In contrary circumstances the co- 
existence of the three groups is possible, and in the 
intermediate case the more favoured pure strain alone 
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will survive. This is what happens when all the char- 
acters of the hybrids are either dominated by those of 
one of the pure strains or are intermediate between 
those of the two pure strains. 

The words "more favoured" and "less favoured" 
are very vague. We have tried to replace them by 
inequalities connecting the vital co-efficients and to 
make a separate study of the action of the three chief 
groups of factors, and the conclusions obtained are in 
complete agreement with actually observed results. 

Two points remain obscure — the causes of the appear- 
ance of mutations and the number of varieties that a 
species may include. The two problems are not inde- 
pendent of each other. Everything would seem to 
depend on the initial and the final distribution of the 
genes in the chromosomes. Physical and chemical 
factors may act in this process, and it may be asked 
whether in critical circumstances a sudden change in 
external conditions may not produce abnormal distri- 
butions, some of which may be advantageous and 
permanent. 

It seems further that the number of stable genetic 
arrangements is considerable and that the phenotype 
may show almost continuous variation. We cannot go 
into details of this mechanism in the chromosomes. 
What is important is the possibility of almost unlimited 
variation in a species by the working of two indepen- 
dent mechanisms — selection and the appearance of 
mutations. 
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